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ABSTRACT 


The primary objective of this report was the development 
of an approximate method of analysis to predict the behaviour of 
multi-story shear wall-frame structures with reasonable accuracy 
and study the behaviour by a test of a reinforced concrete shear 


wall-frame structure. 


The analysis was performed on an equivalent lumped model 
under incremental lateral load and fixed vertical load and traces 
the second order elastic-inelastic response of the structure. 
Conditions of equilibrium are formulated on the deformed structure 
to include the secondary, P-A, effects in columns and shear walls. 
The analysis does not consider the axial shortening of the columns 
and shear wall but the width of wall has been taken into consideration. 
Time effects, temperature changes and shear deformations in the members 


and joints are neglected. 


The analysis uses an iterative procedure. The computer 
program developed for the analysis uses moment area and slope 
deflection equations modified to consider the presence of plastic 
hinges in the members. A method has been developed to derive the 
equivalent lumped structure. Using this analysis, good correlation 


was obtained with other analyses. A design example has been presented. 


Rationalized methods were developed to predict the moment- 
iii 
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iv 


thrust-curvature relationships for the cross-sections of the members. 
The computer program developed for the purpose considers rectangu- 
lar reinforced concrete cross-sections and can consider tensile 


stress jn the concrete and unsymmetrically placed reinforcement. 


A four story, one bay shear wall-frame structure was tested 
under incremental lateral load and fixed vertical load applied at 
the top of column and wall. The behaviour shown by thtis test was 
compared to the behaviour predicted by the present analysis and 


a good correlation was obtained. 
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CHAPTER I 


INTRODUCTION 


The social and economic development of urban areas requires the 
construction of high rise commercial and residential buildings. Recent 
years have seen the construction of many tall buildings and their numbers are 
increasing rapidly. Such construction has proved the necessity for more 
and more knowledge of the behaviour of such buildings and, analyses capable 


of accessing, with sufficient accuracy and rapidity, the overall strength of 


the buildings. 


The planning and design of tall buildings requires the provision 
of sufficient resistance against lateral loads in all directions. This 
stiffness can be achieved in many ways depending upon the height and plan 
of the building. The most usual planning concept uses shear walts to 
provide the required rigidity. In the twentieth century large number of 
tall reinforced concrete buildings have been built with shear walls so 
located that they act as columns in addition to resisting the lateral forces. 
In most cases, in addition to carrying vertical and lateral loads, shear 
walls are used for the purpose of enclosing service areas. As a result, 
Shear walls can take a variety of shapes. The object of this thesis is to 
present an approximate method of analysis for shear wall-frame structures 


and justify this analysis by laboratory test and discussion. 


The modern trend of reinforced concrete design has shown increas- 
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ing use of ultimate load theory. This design method tends to be more 
economical and logical than working stress design because it predicts the 
dependable load carrying capacity of the structure and the design can be 

SO arranged that the appropriate load factor can be achieved. Most of the 
presently available analyses for shear wall frame structures, reviewed in 
CHAPTER II, are limited to the elastic regime and tend to give rise to 
different load factors for different members in a structure. Consequently, 
the information from such an analysis is not sufficient to design by a limit 
design procedure and it does not constitute the true behaviour of the 


structure. 


The theory of "plastic hinges" is now widely accepted to estab- 
lish the behaviour of steel members at failure. However, this type of 
failure is still being investigated for multi-story structures. The 
trial and error method of choosing the plastic hinge configuration in 
the mechanism method of analysis has the serious disadvantage that there 
is no guarantee that the worst possible hinge configuration has been 
chosen. The recent development of electronic computers has made it possible 
for engineers to perform more satisfactory analyses of multi-story steel 


structures. 


Because of the complex nature of behaviour of reinforced concrete 
members, the development of the comparable theories for reinforced concrete 
structures has not been so marked. Since reinforced concrete columns lack 
ductility, an ultimate load theory for concrete structures, in which 
collapse is forced to occur due to beam hinging by creating a strong column 
(H2) 


and weak beam design could be used. This method has the disadvantage of 


the trial and error method of choosing the worst mechanism and leads to the 


Overdesigning of the columns. The most satisfactory and probably eventual 
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analysis method for concrete structure seems to be the technique of analys- 
ing a mathematical model on a computer through each stage in its loading 
history. This method has the advantage that the load history dependence 


of the stress-strain curve can be included in the program. 


In tall buildings, the axial loads in columns are significant 
and produce additional overturning moments commonly referred to as P-A 
moments; where P is the axial load and A is the displacement due to side- 
Sway. The effect of these "secondary" moments is to reduce the capacity 
of the structure to resist lateral loads. These moments must be considered 
in the analysis of tall buildings (D3) . In an unbraced structure the P-A 
effects are significant and increase as the inelastic action progresses 
Causing a rapid softening of the structure. If the members are sufficiently 
Slender,it is conceivable that the failure load could be initiated by elastic 
instability before any portion of the structure reached the inelastic 
behaviour range. On the other hand in a laterally braced structure the 
P-A effects are reduced due to lateral stiffness of the structure and 
hence less danger of elastic instability failure exists. The response of 
a structure is significantly influenced by the shear wall, which provides 


(C5) 


a sort of bracing system. Although the problem has been studied , there 


is no specific way to differentiate between braced and unbraced structures. 


(A2 ) 


It is mentioned elsewhere that "what constitutes adequate bracing in 


any given case must be left to the judgement of the engineer, depending 


on the arrangement of the structure in question". 


For shear wall-frame structures, it has been common practice to 


design the frame for the vertical load and shear wall to resist the lateral 


(B15) 


loads only. .. For a building of relatively small height this assumption 


is adequate but considerable error results when it is applied to high 
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buildings. Therefore this assumption seems to be an over-simplification. 
The behaviour of a building depends on the way it is built and not in the 
way it is arbitrarily assumed to behave. When the frame portion of the 
building is fairly rigid by itself, the interaction between the shear wall 
and frame can result in considerably more rigid and efficient design. The 
redistribution of the lateral forces resulting from the inelastic action 
would influence the interaction behaviour in a shear wall frame system. 
Therefore, a realistic approach to assess the load sharing between shear 


wall and frame is an obvious necessity for a rational and economical design. 


An attempt has been made in this thesis to present an approximate 
computer analysis to trace the second order elastic-inelastic behaviour of 
shear wall frame structures. A brief review of the previous work on such 
Structures is given in CHAPTER II. CHAPTER III presents the proposed 
method of analysis. The assumptions made in CHAPTER III are critically 
examined in CHAPTERS IV, V and VI. The cross-sectional response is dis- 
cussed in CHAPTER IV. CHAPTER V describes the response of the members 
and joints. CHAPTER VI is devoted to the discussion of the simplified 
Structure assumed in the analysis. CHAPTER VII describes the fabrication 
and testing of a four story, one bay reinforced concrete shear wall frame 
Structure. The test and the processing of test data are discussed in 
CHAPTER VIII and the actual and predicted behaviour are compared in CHAPTER 


IX. CHAPTER X summarises and concludes the dissertation. 
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CHAPTER ITI 


REVIEW OF THE PREVIOUS WORK 


7.1) Introduction 


For the investigation of the problem outlined in CHAPTER I, it 
is necessary to employ a method of analysis to trace the response of shear 


wall-frame structure as loading progresses to failure. 


A brief review of the extent of the present information on the 
elastic and inelastic action of multi-story structures is included in this 
CHAPTER. The review considers five different aspects which influence the 
Study of shear wall frame structures. The areas of this survey of litera- 


ture are as follows: 


1. Cross-sectional behaviour in the form of moment-curvature 


relationship and material properties. 


2. Analysis of frames. 


3. Analysis of shear walls. 


4. Analysis of shear walls combined with frames. 


5. Tests on concrete structures. 
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2.2 Moment Curvature Relationship 


An extensive study on plain and reinforced concrete columns was 
carried out in the period 1929 to 1933 under the sponsorship of the American 


(RT) | The most important conclusion of this study was 


Concrete Institute 
that the full concrete cylinder strength cannot be taken as the failing 
strength of a reinforced concrete column,and a factor, 0.85, was established 
to be multiplied to the cylinder strength of concrete. This value of 0.85 


(HT) For analysing his 


was confirmed by the studies made by Hognestad 
test results, Hognestad proposed a new stress-strain relationship for 
concrete in flexure. There is a discontinuity in the slope of the curve 
at its maximum value.and the concreté is assumed to fail at an ultimate 


Serain of 0.00388 


The stress-strain relationship for concrete is the subject of much 
debate, since the mechanism of failure of concrete is not completely under- 


(NI ,D1 ,K3,B125W3) 4 ogregor proposed 


stood. Based on various literature, 
a stress-strain curve for eaneraran which has been used in this thesis. 


This stress-strain curve is described in CHAPTER IV. 


For cross-sections which consist of linearly elastic material 
it is possible and convenient to write a single expression relating moment, 
thrust and curvature. Unfortunately, the response of reinforced concrete 
section is extremely complex and it is difficult to write an M-P-¢ 
expression aS a simple function. Different authors have tried different 
approaches. Most methods use a trial and error procedure, involving con- 
siderable computations, and define the relationships among moment-thrust- 
curvature by a number of discrete points. Such presentation have been 


facilitated considerably by the use of electronic computer. Two methods 
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of such analysis, suggested by Ernst, Hromadik and pivevand ec! un 1953, 


are as follows: 


1. The values of axial load, P, and moment, M, can be obtained 
for a given value of outer fibre compressive strainye 
and using various values of E The moment-curvature and 
axial load-curvature curves for values of E, can be plotted. 
From a series of such curves for different values of E a 
moment-curvature curve at constant axial load can be ob- 


tained. 


2. For a given axial load, P, and a given outer fibre compress- 
ive strain, ae it is possible to find the value of a by 
trial. The process can be repeated for various values of 
- to get a direct moment-curvature curve at constant axial 


load, P. 


(B2) used the Hognestad stress-strain curve and 


Broms and Viest 
derived four equations relating the axial load and the moments to the outer 


fibre strains, for four different possible strain configurations. 


Chang‘ °2) rederived these equations and used them to derive the 
moment-curvature relationship using the first approach. Breen’8!1) 


considerably improved this procedure. 


sree e established contours of equal curvature to define the 
combination of axial loads and moments for a given curvature. By inter- 
polating a sufficient number of contours the moment-thrust-curvature 


relationship can be established. 


All of the above analyses used the Hogenstad stress-strain curve 
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(F2) 


and neglected the tensile stress of concrete. Fowler has suggested 
an analysis using the second approach for getting direct moment-curvature 
relationship for a given axial load. Different stress-strain curves were 


used. The ultimate strain limit was not taken as the criteria for failure. 


(F2) 


An analysis similar to that presented by Fowler using the 
Stress-strain curve proposed in reference including the effect of 
tensile stress of concrete is presented in CHAPTER IV, to define the 


behaviour of reinforced concrete cross-sections. 


2.3 Analysis of Frames 


Several theories and methods have been developed for the analysis 
of single and multi-story frames. The structural action of a multi-story 


frame can be classified as follows: 


1. The frame can resist gravity load without lateral deflection 
of the floors until loads are increased to the load corres- 


ponding to an elastic sway buckling failure. 


2. The frame can resist gravity load without lateral deflection 
at floor level and does not fail in sway mode. The frame 


fails by plastic hinging. 


3. The frame can resist gravity load with lateral deflection of 


floor level. Failure occurs by inelastic buckling. 


4. The frame can resist both gravity and lateral load with 
lateral deflection of floor level and the frame fails by 


instability. This is characterized by a gradually increasing 
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lateral deflection, first under a growing load and later 


after reaching a maximum load, under a diminishing load. 


The various methods for evaluating the elastic buckling load 


(B1)_ 


have been summarized by Bleich The elastic and inelastic solution 


of buckling load of symmetrical three story structures subjected to 
uniformly distributed gravity loads have been described by pyfb2) 

He has also suggested that the solution described for a three story frame 
may be applied to taller frames by first dividing them into several 
three-story tiers. The case, where the failure results from a trans- 
lational instability of the structure is the most practical problem. In 


most cases, much of the frame is strained into the inelastic range before 


the buckling loads are reached. 


Merchant MV) has proposed a Rankine type formula, 


ee +. Ey for the evaluation of inelastic buckling capacity of 

he Ap he 

the structure, hes where hp is a function of rigid-plastic collapse 

load and he is the function of elastic critical load. home cue? has 

justified this formula and it was experimentally verified by Magid (M7). 
Wood! "2) illustrated the loss in structure stiffness in the 


presence of plastic hinges by inserting real hinges in the structure and 
computed "deteriorated elastic critical load". His illustrations indicated 
the tendency for structures to become unstable prior to the formation of 
enough hinges to render the structure into a mechanism. However, it is 

not possible to predict the deteriorated state of the structure without 
advance knowledge of the actual history of the order of hinge formations. 


Even then, only an upper and lower bound to the collapse load can be 
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11 
obtained. 


Recently, computer oriented elastic-plastic analyses of unbraced 
plane frames have been developed by assuming point plastic hinges and 


elastic-perfectly plastic section behaviour. 


(J1) 


Jennings and Majid developed an iterative procedure to per- 
form an elastic-plastic analysis of unbraced frames loaded by static, 
proportional concentrated loads. The solution is derived using a matrix 
technique having joint displacement as unknowns. This analysis yields 
the complete load-deflection characteristic and the order of plastic 


(D2) extended this analysis for cyclic loading 


hinge formation. Davies 
and considered the effects of hinge reversals by a locked hinge with a 


rotational discontinuity. 


Parikh P2) 


formulated a second order elastic-inelastic analysis 
employing slope deflection equations. An iteration procedure was used 

for solving the simultaneous equations. Korn(K5) has also developed a similar 
Solution except that he uses matrix techniquesto solve the simultaneous 


equations. 


Heyman ‘42) has proposed a design method for multi-story frames 
in which beams were designed by plastic methods and columns were designed 
to be elastic up to the design ultimate load. An approximate method to 


calculate the deflection has also been proposed. 


Daniels and Ly623) proposed a subassemblage method to analyze 
unbraced frames subjected to vertical and lateral loads. Isolated parts 


of a multi-story frame were analyzed as independent units with the point 


beosrdau to esevlens afieslq-oftesl> batns ero votNamed ey ltei3g 


brs gapntd ott2zatq anteg pntmuees yd boqofeveb Ased°ov 2 ‘ | 
Wwortversd nordose oftealq ylas | 


i 
p i 7 i. 
( ink = a 
J ie, in 3 \ 
i | ) 
1 — 
i 


wraq 03 swwhso07g ov tisvadt as beqolavab’ cL i 1 


a edt dese yd bebsol 25nsyt boos ydriv bo leteylens aid 26 Tg: } otSe6 Te 
xin tam 6 pneu bev Fiab er notutoz ‘nT ebsof bareriiea ) 

a eblary ef aulsns inl Jenworiny 26 Insmessigeth tntog 
atteslg +o "abv ait bib ottetysiobyedo notgaeitsbet 
($4) soy tvaG an 
6 ai epafd badlaol 6 yd elezraver sent to 2etostts ont 


ae) vet 


entbao| at faye Yoh efeylans vind bobnetxs 


2teyians attesloni- ai te5 fs WSb10 broose 6 bots lume? ud 
uy beau: esy swubso0vg nottsyott aA ano Feu notice feb snore 
i ’, satis 6 beqolsvab oe! 6 agit A prvoy ‘enot Supe Frisia Hele oitt 
i”  gwoonbifynre ots. av foe od zouptnnoed zee adel) 2n sits ‘in 


™~ 


\ 


eens sanep-teiam OF boditem Ap Tesb*5 bs20qexq: 26d (cmt 


a c... saw enmfoo bre abortan ofsept q yd Rege reab ovew amoad 
‘hada odiembxorags nA. fy) seni aty: Fs ae art nae itn 


| sloal — iat 


\ 


of contraflexure in the columns being assumed to be at mid-height. Sub- 
sequently, this model is reduced to a number of subassemblages, each 
consisting of a column restrained by the adjoining beams. Both inelastic 


action and the second order P-A effects are included in this analysis. 


(KT) 


Kloucek suggested a simplified solution for elastic analysis 
of the sidesway of a multi-story closed system. The true multi-story 
system is replaced by a substitute cantilever with the same loading. The 
Substitute member stiffness factors of the cantilever were obtained by 
addtion of all the column stiffness factors on the same floor of the 
original frame. The beams were represented by "knots" at each floor with 
the stiffness factors of the "knots" determined by adding an imaginary 
member having a rotational stiffness equal to sum of all the beam stiff- 
ness factors on one floor multiplied by a factor "A". The factor "A" 

is a fairly insensitive quantity. Many examples were solved to show the 


(L1) 


rapidity and validity of the method. Lightfoot has illustrated this 


method and has solved two examples. 


2.4 Analysis of Coupled Shear Walls 


This section is devoted to some of the literature available for 


the analysis of groups of shear walls coupled by beams. 


The analysis of coupled shear walls is due in large part to 


Rosman and Beck (R3,R4 ,B9) Rosman established the solution in the form 


of a trigonometric series using the principle of least work. In subse- 


(R3) | He derived 


quent work he used a direct mathematical solution 
Solutions for a wall with two symmetric bands of openings, with various 


conditions of support at the lower end (piers on rigid basement, on separate 
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foundations and on various forms of column supports). Two loading cases, 
a uniform wind load, and a point load at the top of the building were 


considered. Beck (89) 


» who used a similar approach as Rosman, presented 
an approximate method of analysis where a continuous system replaces the 
discontinuous system. He assumed that all connecting beams had 

the same distance from each other and, except the end beam, had 

the same stiffness. The end beam had one half the cross-section and 

one half the moment of inertia of a normal connecting beam. He treated 
the single case of two uniform coupled shear walls on a rigid foundation, 


Subjected to a uniformly distributed lateral load. Useful graphs to 


enable bending moment, axial force and deflection have been presented. 


LE) 


Frischman, Prabhu and Toppler have used two different 
methods, both based on the principle of analysis for a rigidly jointed 
framework, to deal with the coupled shear walls. Equivalent column 
method has been used, by which a multi-bay structure can be reduced to 

a single column by the lumping together of all column stiffnesses, with 

a single beam restraint at each story level, given by the addition of 

the beam stiffnesses, to analyze the coupled shear walls. The problem 
was reduced to a single second order differential equation, by using a 
continuous distribution of restraints and loads, to evaluate the bending 
moments in the continuous system. A conventional method using the flexi- 
bility approach was also used. A similar approach was used by Coull 

(C4) 


and Choudhury who have presented curves for the moments and maximum 


deflection. 


2.5 Analysis of Shear Walls Combined with Frames 


Several investigators have considered the problem of load 
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distribution in the combined frames and walls system. Because of the 
high in plane stiffness of floor slabs, most investigators have been 


concerned with the plane frame wall system. 


(R2) 


Rosenblueth and Holz present an approximate analysis of the 
interaction between single shear wall and a rigidly jointed frame 
structure. The shear carried by the frame in any one story was assumed 
to be proportional to the average slope in that story whereas the moments 
and shears in beams supported by the wall were proportional to the flex- 
ural slope of the wall. Newmark's method of successive approximation 


was suggested as a solution. Column shortening and foundation deformations 


were neglected. 


(C1) 


Cardan reduced the problem to the solution of second order 
linear differential equations, giving the slope of the shear wall. The 
properties of frames and walls were assumed to be constant ‘throughout 
the height. Axial deformations were neglected and the wind loads were 


assumed to be distributed throughout the height of the structure. 


Bandet (810) 


suggested the analysis of shear trusses, to replace 
an equivalent shear wall, combined with frames. Power series solution 
was suggested which requires the solution of a set of simultaneous 
equations. Column shortening was considered but axial deformations of 


the beams were neglected and a fixed end foundation condition was assumed. 


A computer program, using the stiffness method of analysis, 


taking flexure, shear and axial deformation into account, was developed by 


Clough, King and Witson'©3) , The axial deformation of beams Was neglected. 
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It was assumed that the building was laid out in a regular rectangular 
grid pattern, with each floor level constrained to translate but not 
rotate under the action of lateral forces. It was further assumed that 
the shear walls were of uniform width throughout the entire height, 


although variations in stiffness of shear walls were allowed. 


(K2) 


Khan and Sbarounis presented a method of successive approx- 
imation on a substitute structure derived for the purpose. of. analysis.The 
derivation of the substitute structure resembled that used by Kloucek‘K1) oy 


Lightfoot ‘t!). 


Forced convergence techniques have been used to arrive at 
an elastic solution. In this way, convergence is tied directly to the 
physical reality that the deflections of the entire structure must lie 
between zero and the free deflection of the shear wall. This method 

will. converge for any combination of structural stiffnesses although 

at different rates. A number of influence curves were presented for 


use in preliminary design. The analysis proposed in this thesis is on 


Somewhat similar lines. 


Goula'&!) 


replaced the frame by a system of rotational and 
translational springs. The finite difference technique was used 
to solve the fourth order governing beam equations giving the wall 


deflection at each floor level. 


A computer analysis, which traces second order elastic-plastic 
behaviour of planar reinforced concrete structures as loading progresses 
to failure, was presented by Clark, MacGregor and Adams (65 +66) | The 
analysis considers axial shortening of the columns and shear walls and 


includes the effect of finite width of shear wall. Time effect and shear 
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deformations on the members and joints are neglected. 


(W7) 


Winokur and Gluck presented a method of analysis for 
asymmetric multi-story structures. They have shown that the rotational 
constraints of the stiffening elements in their planes at the floor 
diaphragms are important factors which are to be taken into account in 


analysis. 


A number of other papers (M4 »B14 sG2 ,J2 ,K4 ,P3,S2 471,72 ,W5 ,W6 ) 


are also available which illustrate the problem and solutions have 


been suggested with various assumptions. 


2.6 Tests of Shear Walls and Shear Wall-Frame Structures 


Though much has been said about the analysis of multi-story 
structures, experimental studies should not be forgotten. Although 
many experimental studies have been carried out on multi-story frame 
Structures, relatively few experimental studies of shear wall-frame 


Structures have been carried out. 


The first major tests were carried out in the United States in 
full scale by atomic explosions at Eniwetok. Abbreviated reports were 


(WI) Approximate methods based 


presented by Whitney, Anderson and Cohen 
On simple strength of materials theories were proposed to predict the 
stiffness and ultimate capacity of simple shear wall structures subjected 


to lateral loads. 


A large number of tests were made by Benjamin and Williams 


(B3 ,B4,B5 ,B6) | on single story brick and reinforced concrete shear walls 
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with and without openings and on one or two story parallel reinforced 


concrete shear wall assemblies connected by diaphragms. 


A model test on a reinforced mortar eight story shear wall- 


frame was described in reference (A4). The stories were each one foot 


high. Photoelastique methods were also used in this investigation. 


A model scaled at 1:64 was tested by Barnard and Schwaighofer 
(B14). establish the width of "slab strip", which acts as the connecting 
media between shear wall coupled by slab floors, and to determine the 
validity of an analysis presented. The walls for the model were cut 
from 1/4 inch thick epoxy sheets and the floor slabs were slotted to fit 
around the walls and then glued to the walls using an epoxy glue. The 


model was tested in horizontal position with the base being rigidly 


fixed. 


Twenty three story models constructed from 1/4 inch thick per- 


(J2) 


Spex sheets were tested by Jenkins and Harrison to compare their 


analysis under bending and torsion. 


2./ Conclusions 


Extensive literature is available on different aspects of the 
problems encountered in dealing with multi-story structures. Some of 
these papers are reviewed briefly inthe previous sections. The review 
indicates that many questions regarding the behaviour of shear wall frame 
structures are yet unanswered. Most of the studies were limited to 
elastic cases only. In many cases the analytical procedures are cumber- 


Some requiring extensive computation work and not suitable for design 
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office use. 


Very few actual laboratory tests have been reported especially 
on reinforced concrete shear wall-frame structures. Many tests have 


been performed to describe the behaviour of the members. 


In view of the above limitations it was felt necessary to derive 
an approximate rapid, elastic-plastic analysis of shear wall frame 
structures, and to check this analysis by tests on large scale reinforced 


concrete models. This is undertaken in this thesis. 
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PART II 


APPROXIMATE ANALYSIS OF REINFORCED CONCRETE 
SHEAR WALL-FRAME STRUCTURES 
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CHAPTER: ITI 


APPROXIMATE ANALYSIS 


3.1 Introduction 


A computer analysis of shear wall-frame structures under 
incremental lateral loads and constant vertical load is described in 
this CHAPTER. The analysis traces the second order elastic-plastic 
behaviour of the structure as the loading approaches the ultimate 
value. To simplify the problem the analysis deals with the reduced 
Structure shown in FIGURE 3.1. The formulation of the analysis and 
the assumptions made are outlined in this CHAPTER. A design example 
is calculated using this analysis in APPENDIX A. The validity of the 
assumptions will be discussed more fully in CHAPTERS IV, V and VI. 

The analysis will be compared to the results of a test of a four story 


frame in CHAPTER IX. 


3.2 Assumptions in the Analysis 


The assumptions made in deriving this analysis fall into 
several major categories. The first group includes those dealing with 
the behaviour of the structural members, joints and plastic hinges in the 
structure... These assumptions are discussed in detail in CHAPTERS IV 


and \V. 
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It is assumed that the moment rotation relation- 
ship for the individual frame members can be 
idealised as elastic-plastic for the beams and 
columns and elastic-strain hardening for the 


shear walls, as shown in FIGURES 3.2(a) and (b). 


All members in the idealised structure are assumed 
to be prismatic. That is, the stiffness and ulti- 
mate moment capacity may vary from member to 
member but not along the length of a given member. 
It is possible to develop moments of opposite sign 
at the two ends of the member, The 

ultimate moment capacity, My? and stiffness. El, 


of the members are derived in CHAPTER IV. 


The stiffness factors C and S commonly used in the 
slope-deflection equations have been assumed to 
remain sufficiently close to 4 and 2 that these 
values can be used. In deriving the M, value for 
the cross-section, however, the axial load is taken 


into account. 


Shear deformations are neglected. 


The joints are rigid prior to hinging of the members 


and do not undergo shearing distortions. 


The plastic hinges in the columns and beams are 


assumed to be point hinges at the ends of the members. 
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Plastic hinges in the walls are assumed to form 


in a segment of preset small finite length. 


The hinges are assumed to have an infinite rotation 


capacity. 


The second and most important group of assumptions deals 


with the derivation of the simplified model used in the analysis. This 


group of assumptions is derived and discussed in detail in CHAPTER VI. 


ce 


It is assumed that the load-deflection behaviour of 


multi-story shear wall-frame structures can be represented by the 


Simplified or "lumped" structure shown in FIGURE 3.1. This assumption 


implies that: 


(a) 


Each floor is assumed to act as a diaphragm of 
infinite rigidity and it is assumed that the 
structure does not rotate about a vertical axis. 
Thus, all points on a given floor level are 
assumed to undergo the same lateral deflection. 
This assumption allows the entire building to be 


considered as a plane frame. 


The shear walls in any one story can be represented 
by a single wall with a moment curvature relation- 
ship defined by superimposing the moment curvature 
relationships of each of the individual walls in 

that story. The resulting moment-curvature relation- 


ship is approximated by a bi-linear curve. The 
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Slope of the second branch of this curve may be 
small but not zero and the ratio of the slope 

of the second branch to that of the elastic portion 
of the curve must be the same in all stories. The 
neutral axis of the wall is assumed to coincide 


with the centroid of the uncracked wall. 


The columns in any one story are assumed to be 
represented by a single column as shown in FIGURE 
3.1. The lumped column is assumed to have a stiff- 
ness and a plastic moment capacity equal to the sum 
of the corresponding values for all the columns 


being lumped together in a given story. 


The beams linking the shear wall to the rest of the 
structure (link-beams) in each story are assumed 

to be represented by a single link beam as shown in 
FIGURE 3.1. The lumped link beam has a length equal 
to the average length of all the link beams in the 
story except that the lumped link beams must have 
the same length in all stories. The stiffness, 
EI/L, of the lumped link beams shall be taken as the 
sum of the link beams in the story in question. The 
plastic moment capacity of the lumped beam shall be 
taken as az M, for all the corresponding beams. The 


term a, a reduction factor, is discussed in section 
Orie se: 


The remaining beams connecting adjacent Columns in 


each floor are assumed to be represented by a single 
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beam as shown in FIGURE 3.1. The two ends of this 
beam are assumed to undergo the same rotations. That 
is, the point of contraflexure in the beam is assumed 


to be at mid-span. 


(f) The lumped beams are assumed to have a length, Lp» 
equal to the average length of all the corresponding 
beams under consideration. The moment of inertia, 

I p> and the plastic moment capacities of the lumped 


beams , Moi p> are given by equations: 


lp oliget lpia) \ oe (351) 


M ra(M M 


PLB ~ part Mppo) tee 


These equations are derived in sections 6.2.1 and 


642.262. anrnerterm-a 1S discussedsinisection-6:242.2. 


(g) The foundation conditions are represented by elastic 
rotational springs at the bases of the column and 
the wall of the lumped model. The spring constants 
are obtained by adding the known or assumed foundation 
spring constants of all the columns or walls in the 


real structure. 


The third group of assumptions deals with the loads to be applied 


to the lumped structure. These assumptions are discussed in CHAPTER VI. 


3. (a) It is assumed that the total vertical load in a story 
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can be applied as a concentrated load. Thus, there 
are no dead and live load moments on the beams in 
the lumped structure. The vertical loads remain 


constant throughout the analysis. 


(b) The lateral loads are assumed to be concentrated at 
the floor levels. The lateral loads are incremented 


during the analysis. 


(c) The loads are assumed to be applied statically and 


act in the plane of the structure. 


The fourth group of assumptions concern special aspects of 


behaviour of entire structure. These are as follows: 


4. (a) It is assumed that there is no out-of-plane behaviour 
of the members such as lateral-torsional-buckling 
etc. and it is assumed that none of the members fail 
due to buckling prior to the attainment of the ulti- 


mate load. 


(b) Temperature changes and differential axial settlement 


or shortening are neglected. 


3.3 Method of Analysis 


To determine the deflected shape of a structure subjected to a 
Static load, the conditions of equilibrium and compatibility must be 


Satisfied. The usual methods for such an analysis are the displacement 
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(or stiffness) method or the force (or flexibility) method. 


In the displacement method the joint displacements are treated 
as unknowns and equilibrium conditions of external and internal forces 
are used to solve for these unknowns. In the force method, the redundant 
forces are treated as unknowns and consistent deformation conditions are 


then used to solve for these unknowns. 


The present analysis uses an iterative technique in which both 
the above procedures have been used to arrive at a solution. For this 
purpose the analytical model has been divided into the frame system and 
the wall system shown in FIGURE 3.1. The procedure used requires the 


following steps: 


A. Apply all the lateral load to the free wall system, The 


deformations are then computed(FIGURE. 3.3). 


B. Force the frame system through the deflections computed 
in step A at each floor level and compute the shear 
induced in the frame due to this force-fitting (FIGURE 


rf 


C. Compute the shears and moments in the wall system in each 
story due to the applied loads and the force fitting 


process. (FIGURE 3.5). 


D. Calculate the deflections, A_:, and rotations, Oi of 


WI 
the wall system due to the shears and moments computed 


in step C. A “forced-convergence correction" was applied 


to obtain the initial trial deformations of the (n#1) 
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cycle. The correction used in this analysis is described 


in reference (K2), 


E. Steps A through D are repeated in order until the desired 


accuracy has been achieved. 


This basic procedure was originally described by Khan and 


Sbarounis 62) 


for structures in the elastic range. It is extended in 
this report to trace the second-order elastic and inelastic behaviour 


of the structure. 


A computer program was written in Fortran IV language for 
IBM system 0S/360 to carry this out. The flow diagram, program 
nomenclature and the listing of the program are presented in APPENDIX 


B. 


The program takes into account the inelastic action of the 
structure due to the formation of plastic hinges and the secondary 
moments caused by the P-A effect. To take this effect into consider- 


ation the equilibrium equations are formulated on the deformed structure. 


The steps in the program are outlined below: 


1. The program starts by reading in all the data necessary 
to describe the problem. This includes the number of 
structures to be considered, the dimensions and properties 
of the structural members in the lumped frame shown in 
FIGURE 3.1, the loading and load increments, and the 


convergence limit. 
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The stiffness, EI/L, of all the beams and columns in the 
lumped frame are computed and the vertical load acting 

on the frame is also computed. This load remains constant 
while the lateral load on the structure is incremented to 
determine the load-deflection characteristics. Each story 
of the shear wall is divided into a desired number of 


segments and systemized for the calculations. 


The entire lateral load is applied to the shear wall and 

the moment at the centre of each segment is calculated 

by considering the wall as a free cantilever. The rotations 
and deflections are calculated at every segment using moment- 
area principles. The vertical displacement at the junction 
of the shear wall and the frame is computed as the product 
of the rotation of the wall and one-half the wall width. 
FIGURE 3.3 shows the deflected shape of the shear wall 
a0) 


under the entire lateral load, where Ow and 


gy, 6! are the rotation, lateral displacement and vertical 
displacement of the wall. The initial subscript W repres- 
ents the shear wall, while i represents the floor under 
consideration and the superscript (1) represents the 
initial iteration cycle. The story height is Hous In 

the latter stage of the analysis a reduced value of the 


wall stiffness is used for any wall segment where a hinge 


has been detected. 


The frame system is forced into the deflected shape of the 


wall as defined by the deformations computed above. This 
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step is shown in FIGURE 3.6. 


Next the frame joint rotations are computed. In FIGURE 3.6 
B. is the beam-to-column joint, W. is the joint at the 
connection of the wall and the right side beam and F. is 
the left end of the left side beam. Under lateral loads 
beam B.F. will have a point of contraflexure approximately 
at midspan. This is enforced in the structural model by 
assuming that the joints B. and F. have equal rotations 


and that the left end support, Fs is free to translate. 


At each joint the slope-deflection equations are written and 
by using the moment equilibrium equations, an expression is 
developed which relates the joint rotation to the deformations 
of the adjacent joints and the member properties. The 
detailed derivation of this equation is given in APPENDIX B. 
A set of equations is generated by applying the moment ~ 
equilibrium equation to all the beam-to-column joints in 
turn. The joint rotations are then computed using the 
Gauss-Seidel Iteration Method. The moments at potential 
hinge locations in the members of the substitute frame can 
now be computed. Next the vertical and horizontal shears 


at the junctions of the wall and frame system are computed. 


At each floor level, the moment to be applied to the wall 
system for the next cycle of iteration, Migs is the 
algebraic sum of the moment, Mpa ° at the end W. of the 
right beam and the moment produced by the beam shear 


multiplied by half the wall width. FIGURE 3.5 shows typical 
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forces at floor level i. In the initial iteration, the 
moment , Mia is assumed to be zero. This is corrected 


in subsequent cycles. 


The net out-of-balance forces are applied to the wall and 


(2) (2) 


the resulting deflections, Aw; » rotations 8.4 


(2) 


and 


vertical displacements 6 are computed as described 


Wi 
in Step (3). 


In order to speed convergence, the frame is not forced into 


(2) 


the deformation mode defined by Avs ,» etc. Instead the 


frame is forced into the position defined by: 


Shee : 
Qf (2) : Sea pny, ee) 
as =. Ui 
Pree) (1) (1) 
AW See. a Wis? (Ne ee (3.4) 
eg aim 0. 
ay ROMER 
ond Sine i a Be es a ee (3.5) 


where Di is the width of the shear wall at the i th floor. 
The analysis now returns to Step (4) and the steps (4) to 
(7) are repeated. For any cycle (after the initial cycle), 


the forcing Equations (3.3), (3.4) and (3.5) are replaced 


by: 
aneiy 
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where n denotes the cycle in’ progress. 


For any cycle (except the initial cycle) convergence tests 
for rotations and deflections are performed at the end 
of Step (5) to determine whether further iteration is 


required. 


The system has converged if the quantities 


( ) \ ( =| ! 1 -| 
aa | Dai ie i a ‘te ; are less 


than a specified limit. 


If the system has not converged as defined above, then 


force F. and moment M (FIGURE 3.5) at each floor 


Bi 
are applied on the shear wall and the Steps from (6) to 


(8) are repeated. 


On the other hand, if the system has converged, the P-A 
effect is considered next. The total axial load at each 
floor level of the structure has been read into the pro- 
gram. The increased story shears due to the P-A effect 

are simulated by increasing the applied lateral loads. 

The additional force Ho is calculated on the basis of the 


deformations obtained from the analysis under the lateral 
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loads H.. The derivation is given in section 3.4. Next 
the lateral loads (H, + H',.) are applied on the shear wall 
system and Steps (3) to (8) are repeated. The deformations 
of the structure obtained from this analysis are compared 
with the deformations obtained from the analysis under 


lateral loads H. only. 


If the deformations obtained from the two analyses do not 
agree (within a specified limit) new additional forces 


H',. are computed as above using the deformations of the 


27 
Structure taken from Step (9). Steps (3) and (10) are 
then repeated until the desirable agreement in deforma- 


tions at each floor level are achieved. 


If the deformations agree, the shear resisted by the 
frame is then deducted from the total lateral load to 


determine the shear to be resisted by the wall. 


The program then proceeds for detection of hinges by 
comparison of computed moments at critical sections to 
that of their plastic moment capacity. If no hinges 

have formed the lateral loads are increased by a preset 
factor times the original starting load and step (3) 
through (11) are repeated. One incrementing factor is 
used in the elastic range and a second smaller increment- 
ing factor is used after hinging has occurred. If, on 
the other hand, hinges are detected for the first time 

in the history of the analysis the locations and moment 


values are printed out and the loads are decreased to the 
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previous values increased by the load increment to be 
used for the inelastic range. When hinges are detected 
in the first loading the previous load is assumed to be 
zero. Following this, steps (3) through (12) are 


repeated. 


From now on the loads are incremented by the inelastic 
factor and steps (3) through (12) are repeated to get 


the load-deflection characteristics. 


The analysis ends if: 1. the desired number of points 
on the load-deflection curve has been obtained; 2. the 
lateral deflections have become larger than a preset 
value; 3. convergence within given limit is not achieved 


in a given cycle. 


The manner in which the inelastic action has been included, 
implies that the stiffness deterioration lags the actual 
stiffness deterioration by one load increment. To mini- 
mize this effect, it is desirable to use smaller increments 


in the inelastic range. 


3.4 Inclusion of the P-A Effect 


The P-A effect plays a dominant role in the behaviour of flexible 
Structures. The vertical load ,P, on the structure, acting through a side- 
Sway displacement, A, produces an additional overturning moment commonly 
known as the P-A moment. Consider a single story structure, which has 


been displaced by an amount A and is subjected to vertical loads, P, 
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as shown in FIGURE 3.7. 


of PA/H, is required in the columns of the frame. 


In this analysis, the P-A effect has been taken into account 


by analyzing the structure under an equivalent lateral force. 


FIGURE 3.8(@), H, 


is the applied lateral load at the i th fbhoor. 


To balance the P-A moment an additional shear 


In 


Due 


to this applied lateral force system, a lateral deflection A; wil] 


be produced. 


The story moment between the i th and (i-1)th floor due to 


the P-A effect is obtained by taking moments about the i th floor 


of all the axial loads above this level. 


at each level, the other story moments are obtained. 


due to the P-A effect, Vas is given by: 


pee eeabe Tone EO m, 


; 
Hej 


where He; is the story height and P. is the total vertical load on the 


Structure above the i th floor. 


The additional lateral load, Haye simulating the P-A effect 


can then be determined and is given by: 


ile LSE ' _ Ps/A,. “4 
Heese lqureigey um aiokl 
_ Paea(Aa (itt) - 
HS (441) 


Similarly, by taking moments 


The extra shear 
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The structure is now reanalyzed under the lateral force system, 


Hyot Hi, 


(without considering vertical loads). 


In FIGURE 3.8(b) 
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the lateral deflection corresponding to the new lateral force system. 
The deflections Ao, etc. are then compared with the deflections, A; 
etc. If these deflections do not agree (within a specified convergence 


limit) the additional horizontal force, H: due to the P-A effect 


a 
for the next cycle is given by : 
i MET ae ZC eT erga eee 
a Nei Ns (441) 


The structure is again analyzed under the lateral force 
System H. + Ho. The process is continued until the changes in 


lateral deflections are within the convergence limit. 


The vertical load at each floor level has been read into the 
program. In the present method of simulating the P-A effect a constant 
vertical load has been used throughout the analysis. In the actual 
structure the girder shears will introduce tension in one column and 
compression in the other. This will not change the total vertical 


load or the gross P-A effect. 


3.5 Functions of the Main Program and Subroutines 


To perform the computer analysis, several subroutines were 
employed in conjunction with the main program. A brief description 


of these subroutines are presented here. 


In the 'MAIN' program, the number of problems to be solved 
is read in. The Subroutine "SR" is called. This subroutine reads in 
the data required for the analysis of the first structure. It computes 


the co-ordinates of each floor from the base and the stiffness of all 
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the members. All the input quantities are printed out. Control is then 


returned to the 'MAIN' program. 


The maximum number of load increments to be used in the 
analysis is specified. The initial load system (i.e. first increment) 


is applied to the shear wall system. 


Subroutine 'BAKA' is called. 'BAKA' computes the distri- 
bution of lateral load between the frame and the shear wall systems. 
The moments and deformations at each segment of the wall are computed. 
The convergence formula is then applied (except in the first cycle) 
to the deflections and rotations. The deformations as computed by 
the convergence formula are enforced on the frame system. The joint 
rotations of the frame system are computed in Subroutine 'FRAME' by 
the Gauss-Seidel Iteration Method. In 'BAKA', the moments at all 
locations are compared with the plastic moment capacities (except 
during the first load increments). If a hinge is detected, the moment 
is set equal to the plastic moment capacity for all subsequent load 
increments. The deformations are compared with those of the previous 
cycle of iteration. If the deformations are within preset limit of one 
another, control returns to the 'MAIN' program. If the deformations 
are not within the preset limit, the process is repeated within ‘BAKA’. 
In this manner the forces and deformations (elastic) for a given loading 


condition are completely known. 


Equivalent horizontal forces simulating the P-A effect are 
computed. These are added to the initial horizontal load and the total 


forces are then applied to the structure. The structure is analyzed by 
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recalling Subroutine 'BAKA'. The deflections obtained from the two 
analyses (i.e. with and without axial load) are compared; if they are 
within preset limit, the results of the last analysis are printed out 


by Subroutine 'ROFA' and 'SR3', . 


On the other hand, if the deflections are not within preset 
limit, new equivalent horizontal loads, are computed on the basis of 
the latest set of deflections. The process is repeated until the deflec- 


tions have converaed. 


Inelastic action of the shear wall is considered in the 'MAIN' 
program. At each segment of the shear wall, the computed moment is 
compared with the plastic moment capacity. If the moment exceeds the 
plastic moment capacity of the story, the moment of inertia of the 
plastified wall segment is reduced according to the M-o diagram. This 


reduced stiffness is used in subsequent steps in the analysis. 


In Subroutine 'SR1', the formation of plastic hinges in the 
members of the frame system is detected. If a hinge is detected, the 
location of the plastic hinge is printed out with the magnitude of the 


moment. Control then returns to 'MAIN' program. 


If the structure is still elastic, the horizontal load is 
incremented and the structure is reanalyzed. On the other hand, if a 
hinge has formed at some location in the structure, the horizontal load 
is decremented so that the structure is still elastic. The structure 


is then reanalyzed with this reduced load. Subsequent increments in 
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load are reduced in magnitude to trace the inelastic range as closely 
as possible. The program stops for the reasons specified earlier. 


It then returns to beginning to analyze the next structure. 
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CHAPTER IV 


RESPONSE OF CROSS-SECTIONS 


4.1 Introduction 


The theoretical analysis presented in CHAPTER III assumes 
that the behaviour of cross-sections of the members can be represented 
by bi-linear moment-curvature relationships. A method of tracing the 
theoretical load-moment-curvature (P-M-o) response of reinforced 
concrete cross-sections is presented in this CHAPTER. This analysis 
predicts the ultimate capacity of the cross-section and the response 
of the cross-section at loads less than the ultimate load. It differs 
from previous analysis by considering the tensile stresses in the 
concrete. The theoretical load-moment-curvature diagrams were used 
in the computation of moments and thrust from the curvatures measured 
in the frame test. Approximate load-moment-curvature diagrams based 
on the theoretical diagrams are also presented in this CHAPTER. These 


are intended for use in design. 


Throughout this CHAPTER the abbreviation P-M-¢ will be 


used to refer to the term "load-moment-curvature". 


4.2 Analysis of Load-Moment-Curvature Response of Reinforced Concrete 


Cross-Sections 
For the purpose of this report, moment-curvature relationships 
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will be computed for the case of constant axial load and changing moment 
and curvature. For a known value of axial load and curvature, it is 
possible to find the value of moment by trial and error. This can be 
repeated for different values of curvature to get the entire moment- 
curvature relationship for the given axial load. The procedure can be 
repeated for as many different values of axial loads as necessary to 

get a set of moment-curvature characteristics for the section. If 


desired, interaction diagrams can be prepared from these curves. 


The following assumptions were made in deriving the P-M-o 


relationships for reinforced concrete section: 
1. The section is rectangular. 


2. The stress-strain curves for concrete and reinforcing steel 


are those presented in section 4.2.1 and 4.2.2, respectively. 


3. The maximum stress developed in concrete is 0.85 times the 
Standard 28 day compressive strength. This has been veri- 


fied by extensive tests (H1) 


4. At all loads, plane section perpendicular to the axis of 


the member remain plane. 


5. Creep and shrinkage deformations are neglected. If 
desired, creep deformations can be approximated by 
increasing the short-time strains. This approximation 
is discussed in references (M2) and (M6). 


4.2.1 Concrete Stress-Strain Curve 


For the purpose of defining the stress-strain 
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response of concrete, the stress-strain pieen ae defined by equations 


(4.1) and (4.2) are assumed to define the behaviour in compression and 


tension, respectively, in members loaded in flexure or combined axial 


load and bending. The relationships are illustrated in FIGURE 4.1 for 


concrete having a compressive strength of 4000 psi. 
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f is compressive strength of concrete loaded in flexure 
and is assumed to be 85 percent of the concrete strength in direct 
compression #2), 


The stress-strain curves were derived to satisfy the following 


requirements: 


The stress-strain curve should be continuous in tension 
and bentraaet ones In applying the stress-strain 
curve in this analysis the ACI code value of the modulus 
of elasticity was used although any other convenient 


value could be used. 
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2. Concrete loaded in combined bending and compression is 
assumed to fail at a limiting strain of 0.0035 'K35C7) 
Under these loadings, concrete will generally fail at 
strains greater than this value. If it does this 
assumption will always result in computed ultimate 
moments which are equal to or less than the greatest 
moment, the section can carry. This can be seen from 


an examination of P-M- diagrams such as those shown 


in FIGURE 4.3. 


3. There should be reasonable agreement between measured 
and computed compression stress-strain curves for strains 


less than 0.0035. 


4. There should be reasonable agreement between the measured 


and computed tension stress-strain curves (W3) 


5. The maximum tensile strength in a split tension test is 
(N1) 


assumed to be Ls = 7/7 i ; 

6. The modulus of rupture, computed as lie = 6M/bt* is 
about 40 percent greater than the tensile strength of 
concrete. That is, the modulus of rupture should be 


about 10 v ie (N1) : 


4.2.2 Reinforcing Steel Stress-Strain Curve 


The assumed steel stress-strain for reinforcing 
steel is shown in FIGURE 4.2. It is assumed that the yield point 


and modulus of elasticity are identical in tension and compression. 
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The effect of strain hardening was neglected. In view of the large 
deformations attendant to the development of the stresses in the 


Strain hardening range, this assumption is fairly reasonable. 


4.2.3 Computer Analysis for Load-Moment-Curvature Relationship 


As the axial loads and moments are varied the cross- 


section can have three different strain configurations: 


1. The full section is under compression. 


2. One face of the section is under compression and the 


other is under tension but uncracked. 


3. One face of the section is under compression and the 


other is under tension and cracked. 


The relationships between axial load, fibre strains, 
moment and curvature for each of these cases are derived in APPENDIX 


e. 


A computer program was written to generate the moment- 
curvature relationship for a given axial load. The program was 
written for IBM, 0S/360 system in Fortran IV language. The flow 
diagram, the nomenclature used in the program and a listing of the 
program are aiven in APPENDIX C. The principal steps in the pro- 


gram are described below: 


1. The number of cross-sections to be analyzed and the 
material properties, cross-sectional dimensions and 


axial load, P, on each section are read. If the ten- 
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Sile stresses in the concrete are to be neglected, the 
ultimate tensile strain is specified as zero. The area 
of steel in each layer and the distance from the tension 
face to the respective centroids is also read. The 


pertinent data are then written out. 


The ultimate capacity of the section under pure axial 
load is then calculated in the form P/btf .". This is 
the upper limit of axial load on the section. If this 
axial load is less than the axial load under consider- 


ation, computations are terminated. 


A value of curvature times depth, ot, iS assumed. 


A trial value of extreme compressive face strain, Eqs 
is assumed and the extreme tensile face strain corres- 
ponding to the curvature assumed in step (3) is then 


computed. 


The subroutine "PM" is then used to find the corres- 
ponding strains in each layer of steel and the strain 
domain for the calculation of axial load and moments. 
This subroutine then calculates the load, P/btf and 


the moment, Mbt", for the assumed strain distribution. 


The calculated value of the load, P/btf is then 
compared with the value of P/btf ." for which the 
P-M-¢ curve is being computed. If the difference is one 


percent or less, it is assumed that a solution has been 
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found which satisfies P, M, ot, and Eq: The com- 
puted moment and curvature represent one point on 

the load-moment-curvature diagram for the cross- 
section. For values of P/btt. less than 0.05, an 
acceptable solution is assumed to exist when P/btf 
differs from the given value by 0.001 or less. If 
the convergence is not satisfactory the extreme 
compressive fibre strain is increased by 0.000001 and 
steps (4) through (6) are repeated in order until 


the loads converge. 


7. Following convergence in Step (6) the value of «curvature 
is then increased to find another point on the P-M-¢ 
diagrams and steps (4) through (7) are repeated in 
order. This is continued until the maximum compression 
strains reach 0.0040 at which time the program moves 
on to consider another section. The limiting strain 
of 0.0040 was taken so that a point near the strain 


of 0.0035 could be found in the data. 


This program was also modified to generate P-M-o diagrams 
for a given cross-section subjected to a series of values of P/btf 
ranging in increments of 0.1 from 0.0 to the ultimate axial load 
capacity of the cross-section. The maximum moment from each of these 
curves is stored and used to construct an interaction diagram for 


the column cross-section. 
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4.3 Effect of Various Parameters on Load-Moment-Curvature Relationships 


Prior to discussing the approximate P-M-o relationships assumed 
in the analysis the main parameters affecting this relationship will be 
Studied. For this study a symmetrically reinforced, rectangular cross- 
section, 2 percent longitudinal reinforcement in each of the two faces, 
4,000 psi concrete strength and 50,000 psi steel strength with a cover 


ratio of 0.1 was chosen. The basic section is illustrated in FIGURE 4.3. 


4.3.1 Effect of Axial Load 


P-M-@ diagrams for the basic section with axial loads 
Ratios. P/btf » equal to O, 0.2, 0.4, 120, IZ and 134 are plotted 
in FIGURE 4.3 to show the effect of axial load on the moment 
curvature relationship. The latter two curves correspond approx- 
imately to eccentricity ratios, e/t, of 0.1and 0.05 respectively 
while P/btf equal to 0.4 is approximately the balanced load. 
The ultimate axial load capacity, P /btt. , of the basic section 
was 1.55. The solid lines in FIGURE 4.3 correspond to compression 
failures, the broken lines to a balanced failure and the dashed 
lines to tension failures. The interaction diagram for this 
column is plotted in FIGURE 4.4 with a solid line. The horizontal 
lines in this figure correspond to the moment-curvature diagrams 
in FIGURE 4.3. It can be seen that the stiffness and ultimate 
moment capacities depend on the axial load on the cross-section. 
For columns failing in tension the P-M-o diagrams are essentially 
elastic-plastic. This is not true for columns failing in compression. 
It can also be seen that the rotation capacity of the cross-section 


decreases as the load is increased. For high values of axial load 
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the rotation capacity is very small. 


4.3.2 Effect of Cross-Sectional Properties 


The basic cross-section was reanalyzed assuming the yield 
Strength of the reinforcement was 40,000 psi. The moment-curvature 
diagrams for P/btf. equal to 0.3 and 0.8 has been plotted in 
FIGURE 4.5 for comparison. As can be seen the yield strength of 
reinforcement does not have any effect on the moment-curvature 
relationship prior to the yielding of reinforcement. After the 
reinforcement yields the change in moment capacity is approxi- 
mately constant and varies roughly in proportion to the yield 
Strength of reinforcing steel. It should be noted that for com- 
pression failures the initial non-linearity in the moment-curva- 
ture diagram frequently results from the yielding of the compression 


reinforcement. 


To study the effect of concrete strength the chosen 
section was reanalyzed for 3000 psi concrete strength. Moment- 
curvature diagrams for the two concrete strength have been 
compared in FIGURE 4.6 for values of P/bt equal to 1020 and 
2040 psi. It can be seen that the reduction in moment capacity 
due to decrease in concrete strength is more for high value of 
curvatures and axial loads. By reducing the concrete strength 
the wa hue. of qs Prt y/F, is increased so that the axial stress of 
1020 psi is much closer to the balanced load in the case of 3000 
psi concrete than it is for 4000 psi concrete. As a result a 


reduction in rotation capacity was exhibited as can be seen in 
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FIGURE 4.6. The intial slope of the curves is reduced a little 
as the concrete strength decreases. This corresponds approxi- 
mately to a reduction in the portion of the EI term which depends 
on the modulus of elasticity of concrete. The rotation capacity, 
stiffness and moment capacity are reduced with a reduction in 


material strength. 


To study the effect of the reinforcement ratio the 
chosen section was analyzed for 4 percent and 1 percent total 
longitudinal reinforcement. The comparison of the two analyses 
is shown in FIGURE 4.7. It is noted that the ultimate capacity 
of the section is greatly effected by reinforcement ratio. The 
change in the slope of initial part of the curve shown in 
FIGURE 4.7 is partly due to the reduction in percentage of 
reinforcement. The effect of cracking of concrete will be 
discussed later. It can also be seen that rotation capacity 
for the beam section (P/btf? = 0.0) is larger for smaller percen- 


tage of reinforcement. 


To study the effect of the cover ratio the chosen 
section was analyzed for ratios of d'/t equal to 0.15 and 0.10 
for fifteen different axial loads. The values of Mu correspond- 
ing to each cover ratio can be compared in FIGURE 4.4 in the 
form of interaction diagram. It can be seen that the cover 
ratio has a significant effect on the moment capacity of the 
cross-section for low axial loads. Under pure moment the ultimate 
moments vary essentially in the ratio of the distances between 


the layers of reinforcement. The effect of cover on moment 
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Capacity diminishes as the axial load increases. The cover ratio 
does not have any effect on the ultimate capacity under pure 


axial load © 


P-M=> diagrams for the values of P/btf equal to 0.3 and 
0.8 and cover ratio of 0.1 and 0.15 are compared in FIGURE 4.8. 
It is observed that the increase in cover ratio reduces the stiff- 
ness of the section especially at low axial loads. This is due in 
part to the reduction in the moment of inertia of the steel about 


the centroid of the column. 
4.3.3. Effect of Tensile Stress in the Concrete 


Moment-curvature diagrams computed with and without 
tensile stresses in the concrete are compared in FIGURE 4.7(a) 
and (b) for two sections having one and four percent longitudinal 
reinforcement for axial loads, P/btf equal to zero and 0.5 . 
The tensile stress in the concrete affects the initial part of 
the curve prior to cracking but has no significant effect on the 
ultimate capacity of the sections considered. The kink in the 
M-d curve of FIGURE 4.7(b) for P/btf¢ equal to 0.5 is observed 


due to the transition of axial load from inside to outside of the 


kern of the section. 


The effect of the tensile stress diminishes as the 
axial load and/or the steel percentage increases. The behaviour 
portrayed in FIGURE 4.7(a) for the pure bending case (P/btf? = 0) 


is quite typical of the behaviour expected from a typical under- 
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reinforced beam section. In the case of extremely under-reinforced 
cross-sections, the cracking moment may exceed the ultimate moment 


of the cracked section resulting in an unstable P-M-¢ relationship. 


4,4 Load-Moment-Curvature Relationship for Use in Approximate Analysis 


of Reinforced Concrete Shear Wall-Frame Structures. 


In the analysis presented in CHAPTER III it is assumed that 
the load-deflection curve of the. strweture can be computed by using an ideal- 
ized elastic-plastic moment-curvature (M-¢) relationship for the beams 
and columns and an elastic-strain hardening moment-curvature diagram 


for the shear walls as shown in FIGURE 3.2(a) and (b) respectively. 


4.4.1 Derivation of Approximate Load-Moment-Curvature Relation- 
ships for Use in the Analysis 


The assumed P-M-¢ diagram for a cross-section may be 


defined by establishing four parameters: 


1. The initial slope or EI value. 
2. The moment, Mp, at which the curve bends. 
3. The slope of the second branch of the curve. 


4. The limiting rotation capacity of the section. 


If the theoretical P-M-¢ curves are available, each of 
these terms can be established by fitting the approximate curve 
to the theoretical diagram. This procedure was used in deriving 
the P-M-o curves used in the analysis of the test frame. In the 
more general design case, however, each of these values must be 


estimated at least approximately. 
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1. Value of EI - The stiffness or effective EI 
of a structural cross-section at any given curvature will be 
given by the ratio M/d. Using P-M-¢ curves similar to those 
developed earlier in this CHAPTER, it is possible to define 
secant or tangent EI values. As shown in FIGURE 4.7, the actual 
EI varies significantly from this at any stage prior to 
yielding. This variation may be due to the reinforcement ratio, 
Pye the stage of cracking, the axial load level, or a number 
of other variables discussed in SECTION 4.3.2. Traditionally, 
structural analysis of concrete structures is based on the EI 
of the uncracked concrete section. It can be concluded from 
FIGURES 4.3, 4.5 to 4.8 that in the normal frame, where beams 
are cracked and have one percent tension reinforcement while 
the generally uncracked columns have a total of about two to 
four percent longitudinal reinforcement, the EI of the un- 
cracked concrete section is not sufficiently accurate. This 
problem has also been discussed by Okamura et ar, (01), 

For the design of slender columns the ACI column 


committee (9) has proposed the following value of EI: 


where E. and E. are the modulus of elasticity of concrete and 
Steel respectively. Je and I. are the moments of inertia of 
the concrete and the steel, respectively, about centroidal axis. 
As shown in reference (M9), however, the equation (4.3) tends to 


underestimate the effect of the concrete term and overestimate 
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the steel term. A somewhat better approximation would be: 


Piha ot eS ees Wie ent ee aed) 
: SuaS 
This value is plotted with fine radial lines on FIGURES 4.3, 4.5 
to 4.8. The examination of these lines reveals that the equation 
(4.4) gives the same trends as the main variables discussed in 


section 4.3. 


For the preliminary calculations, the EI values of the 
members of the structure could be approximated, using the equation 
(4.4), assuming rectangular beams with one percent reinforcement 
in each face and assuming the columns have two percent reinforce- 
ment in each face. Shear walls may be assumed to have one percent 


reinforcement. 


2. Value of M_ - The ultimate moment capacity of the 
cross-section varies significantly as shown in FIGURES_4.3 to 4.8. 
This variation may be due to axial load level, reinforcement 
ratio, Py» or a number of other variables discussed in section 
4.3. The ultimate moment capacities, My? for a given reinforced 
concrete section can be obtained by drawing interaction diagram 
of the form shown in FIGURE 4.4. However, it is not always 
possible to draw such diagrams easily. For simplicity in pre- 
liminary calculations, the M values of the members of structure 
could be taken as equal to the ultimate moment of the section 
for the loading considered as computed by ACI code. This is 


Somewhat a conservative estimate. 
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3. Slope of Second Branch of the Curve - For beams and 


columns the slope of the second branch of the P-M-¢ curve is 

assumed equal to zero. This approaches the truth for loads less 
than the balanced load as shown by the dashed lines in FIGURE 4.3. 
For columns failing in compression, plotted in solid lines in 

FIGURE 4.3, it is necessary to arbitrarily approximate the theoreti- 
cal curve with two straight lines. For the shear wall this slope 

is assumed to be positive, non-zero and small. This will generally 
be a reasonably good assumption for practical shear walls near 


the base of a building. 


4. Rotation Capacity - For simplicity an infinite 
rotation capacity will be assumed for the section. This is 
reasonably true for sections failing in tension in view of the 
large rotational capacities of such sections as can be seen in 
FIGURES 4.3, 4.5 to 4.8. It is not good assumption for columns 
failing in compression. For this reason hinges should not be 
allowed to occur in columns unless special binding is provided 
in the region of the hinge. This is discussed more fully in 


reference (C6). 


4.4.2 Discussion of Assumed Load-Moment-Curvature Relationship 


The load-moment-curvature relationship was assumed to 
be elastic perfectly plastic for the columns and beams and an 
elastic-strain hardening relationship was assumed for the wall. 
The slope of the strain-hardening range will generally be assumed 


to be small. 
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Several moment-curvature diagrams showing the effect 
of various parameters have been presented in SECTION 4.3. In 
general it is concluded that a normal beam section closely 
approximates the assumed elastic-perfectly plastic moment- 
curvature diagram with a long range of rotation capacity. 
Therefore the idealization made here is quite reasonable since 
the error involved is very small. Furthermore the rotational 
Capacity of a normal beam is very large and will not limit the 
strength of the ee nucture, (Ce), 

For axially loaded sections failing in tension, that 
is, for section in which the reinforcement yields before the 
concrete crushes, the moment-curvature diagram was found to 
approach the idealised diagram. Therefore this assumption seems 
to be quite reasonable for the axial loads below the balanced 


load of the section. 


For column sections it was found that the moment- 
curvature diagrams tend to move away from the elastic-perfectly 
plastic relationship as the axial load is increased. The 
sections failing in compression exhibited curved moment-curvature 
diagrams rather than bi-linear diagrams. In addition, these 
sections had a very limited rotational capacity. For the columns, 
the sections could be so chosen that a weak beam and strong 


column design is achieved for the building. 
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CHAPTER V 


RESPONSE OF MEMBERS AND JOINTS 


The shear wall frame analysis presented in CHAPTER III makes 
a number of assumptions which deal with the behaviour of the individual 
members in the structure. These assumptions are listed and discussed 


in this CHAPTER. 


5.1 Effect of Axial Load on Member Stiffness 


The stiffness factors C and S commonly used in slope deflection 
equations have been assumed to remain sufficiently close to 4 and 2 that 
these values may be used. In deriving the Mo value for the cross-section 


however, the axial load should be taken into account. 


The slope deflection equations for a generalized member, 


Subjected to end forces, P, only, can be written as: 


a z ecu Ohl 
where My is the moment at the end A of the member AB. Bq and 8, are 


the slopes at the end A and B, respectively, pis the sway angle of 
the member AB and L is the length of the member AB. The C and S are 


given by equations (5.2) and (5.3). 


SinKL___- _KLcosKL WOT ey 
C= KL (3 %coskL=—KLSinkL’ 22) 
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For zero axial load the terms C and S have the values of 4 and 


2 and equation (5.1) reduces to: 


El 


Mequs (46, + 26, - 6?) See (rd) 


A ple 


When the members are subjected to axial loads the moments 
induced by the deflections affect the distribution of moments. For 
an elastic member the values of the terms C and S vary as a function 


of the parameter L / P/EI. 


For the beams the assumption that C and S remain constant 
is essentially true since the axial load in the beam generally remains 
very close to zero. In the case of shear walls the ratio L vy P/EI is 
Small since the stiffness is high and the change in the values of C 
and S are not of significant importance. However, this assumption 


may be of concern in the case of columns. 


For a rectangular reinforced concrete column the value of 


Ba is given by: 


Pa ee t bt 


where K is a constant which varies between 1.0 to 1.8 for pure axial 


load. For eccentrically loaded columns K will be lower. For this 
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discussion assume that K equal to 1.0 and EI is approximately equal 


to 1000 fr. btr® where r is the radius of gyration. Thus: 


SOU Ry: ' supe —_ 
T To0der" bt r a 10S, 
c Yr 
Therefore L / P/EI = 0.0316 = bee Gy 


Assuming that r equal to approximately 0.3 t for the rectangular section 


of depth, t, the equation (5.6) reduces to: 


Ais Saale : onan (oie) 


The value of L /Y P/EI given by equation (5.7) increases as the 
value L/t increases. In a multistory building L/t seldom exceeds 20, (M9) 
The value of L Y P/EI will be equal to 2.10 and the corresponding values 
of C and S are 3.3745 and 2.1699. Therefore for the case of equal 


rotations at the two ends of the member and zero sway angle the assumption 


will be unconservative by the ratio 6/(2.1699 + 3.3745) = 1.08. 


More than 90 percent of the columns in buildings will have 
L/t less than 8 (M9) The value of L / P/EI given by equation (5.7) will 
be equal to 0.84 and corresponding values of C and S are 3.9050 and 
2.0240. Thus for most building columns this assumption will be uncon- 
Servative by less than one percent. Therefore this assumption will be 


quite reasonable for all practical purposes. 


5.2 Effect of Shear Deformation 


The shear deformation may be of concern in the case of members 
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developing large shear such as shear walls. The relative magnitude of 
the flexural and shearing deformations can be estimated in the following 


analysis. 


upon 
For a shear wall of height, h, acted by a shear force, V, the 


deflections due to bending and shear will be given by equations (5.8) 


and (5.9): 
3 
_ Vh 
dm = 3ET eS 
— 1.2Vh 


where A is the cross-sectional area and I is the moment of 
inertia of the section. The term 3 in equation (5.8) assumes that 
the wall acts as a cantilever and the term 1.2 in equation (5.9) 
is based on a rectangular section. E and G are the modulus of elasticity 


and modulus of rigidity respectively. 


Substituting G = E/2(1+u) and I = Aré where r is the radius 
of gyration and can be taken as 0.3t for rectangular section and where 
u is Poisson's Ratio which can be assumed to be equal to zero for 
reinforced concrete, the ratio of the deflections given by equation 


(5.8) and (5.9) is: 
— 2 eee 8 5.10 


For a variation of h/t from 1 to 2 for a one story wall the 
ratio of flexural to shear deflections given by equation (5.10) will 
vary from 1.54 to 6.17. In this case the shear deflections are Sig- 


nificant. In a multi-story building, however, since the wall generally 
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acts as a cantilever for several stories the ratio h/t will generally 
exceed 5. In this case Ay/ A, will exceed 38. In this case the shear 
deflection is very small in comparison to flexural deflection and can 


be ignored. 


In the case of columns the factor 3 in equation (5.8) should 
be taken as 12. The coefficient 1.54 in equation (5.10) reduces to 
0.385. For h/t equal to 8 and 20, respectively, the ratios of flexural 


to shear deflection are 24.7 and 154. 


Thus the shear deformations will generally be small in 
comparison to bending deformations in columns and will tend to be smal] 
for walls. It is reasonable therefore to neglect the shear deflections 
in approximate analysis. The comparison of test results also shows 
that the shear deflections of columns and walls are not of considerable 


importance. 


In the examples of 10 and 20 story buildings presented in 
reference (G2), the shear deformation was found to be greatest in 
the bottom part of the shear wall. These examples are also analyzed by 


approximate analysis and have been compared in CHAPTER VI. 


9.3 Behaviour of Plastic Hinges 
In CHAPTER III the following three assumptions were made to 
define the properties of plastic hinges in the members: 


(a) Plastic point hinges form at the centre of the beam- 


column joint and at the face of the wall in the case of 
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beam to wall joints. 
(b) In case of shear wall the hinges form in pre-specified 
small finite width. 


(c) Plastic hinges have infinite rotation capacity. 


In the frame test described in PART III of this report the 
plastic hinges in the beams occurred at about 0.25d from the face 
of the columns and the face of the wall as shown in FIGURE 5.1 which 
Shows the third floor beam in the test specimen. As a result the 
moments at the centreline of the column were about 8 percent higher 
than the computed plastic moment. The assumed hinge location leads 
to a conservative estimate of the hinging load. The degree of con- 
Sservatism is not serious for normal frame structures, but can be in 
the case of the exterior wall of a "tube-in tube" type of building 
where the clear span of the beam is generally less than 80 percent 
of the beam span measured centre-to-centre of the columns. In this 
case it would seem reasonable to increase the plastic moment capacity 
of the beams to be used in the analysis by the ratio of the centre- 


to-centre span divided by the clear span. 


The assumption that the hinges can undergo any needed 
rotation is reasonable for the beams which can generally develop very 
large curvature, at hinges, as shown in FIGURE 4.3 and 4.7. The 
assumption may be far from the truth in the case of columns failing 
in compression, however, since the rotation capacity of Such members 
is small as shown in FIGURE 4.3. Walls fall in between beams and 
columns but generally will fail in tension and thus generally will be 


able to develop large curvatures and rotations. Thus the analysis will 
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lead to a reasonable prediction of frame behaviour provided column 
hinges do not form. This will occur if a strong column and weak 


beam design is used. 
5.4 Behaviour of Joints 


In CHAPTER III it was assumed that the joints were rigid prior 
to hinging of the members at the joints and did not undergo any shear- 
ing deformations. The shears induced in the joint by the compression 
and tension forces at the end of the beams frequently lead to diagonal 
cracks within the joint as shown in FIGURE 5.1(a). The forces causing 
such a crack are shown in FIGURE 5.2. For such cracks to occur the 
joint must undergo a shearing distortion which leads to the reduction 
in the joint stiffness. As a result of assuming stiff joint, the 
analysis will lead to an unconservative estimate of the ultimate loads 
and deflections. If ties are provided in the joint to transmit the 
Shears after inclined cracking in the joint, the shearing distortion 
of such joints will not be significant. Ties were provided in the beam 


(H4 ) have dis- 


to column joints in the test frame. Hanson and Conner 
cussed the strength of such joints. To take account of the flexibility 
of the beam to wall joint, an equivalent length method, in which beam 
length is extended in the shear wall beyond the face by 1/2 of the 

beam depth, has been suggested in reference (M10). As will be seen 
later, in the analysis of test frame in CHAPTER IX, extending the beam 
length of the test frame by 1/2 the beam depth resulted in less than two 
percent reduction in the beam moments and little difference in the 
overall behaviour is found. Therefore it can be concluded that the 


assumption of rigid joint is reasonable for an approximate analysis. 
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FIGURE 5.1 


(a) 


(b) 
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FIGURE 5.2 FORCES ON BEAM - COLUMN JOINT 


VY 
=i 
a 
i 


—— 


—- 


~ 


NE ee A eR 
eS ee ee me 


a on 


— — | oe 2 


x 


“oars it so 4 hi As te 
rm 
mn a i nen 


CHAPTER VI 


DEVELOPMENT OF THE LUMPED EQUIVALENT 
STRUCTURE 


6.1 Introduction 


In this CHAPTER, an attempt has been made to determine the 
arrangement of a lumped or equivalent structure, so that the analysis 
developed in CHAPTER III can be used to predict the ultimate load 


of the real structure with reasonable accuracy. 


CHAPTER III contains a preliminary discussion of the simpli- 
fied or equivalent structure. The lumped model used for the analysis 


(K1) (G2) 


was Similar to the elastic model used by Kloucek , Goldberg and 


Lightfoot ‘t1) 


for the analysis of plane frame. An extensive discussion 
of the lumping procedure has been presented by Kloucek. Also a large 
number of examples are presented to show the validity of the procedure 


(K2) 


in the elastic range. Khan and Sbarounis used a similar approach 


for elastic shear wall-frame structures. 


In developing the lumped structure, the chosen examples will 
be analyzed by the approximate method developed in CHAPTER III and the 
response so obtained compared with that obtained from a more rigorous 
analysis of the complete structure, using the method developed in 
reference (C6). The load deflection curves for the building or frame 
under consideration will be compared to show the accuracy obtained by 
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using the lumped model. 


6.2 Development of the Lumping Procedure 


In the discussion of the lumping procedure the frame members 
will be considered first. It has been assumed that the shear wall 
is forced to undergo the same sway deformation as the structural 
frame. FIGURE 6.1 shows the floor plan of a typical building with 
the vertical supporting structure consisting of walls and frames. The 
lateral force is assumed to be applied parallel to the plane of the 


frames. 


A symmetrical rectangular portal frame subjected to a hori- 
zontal load is shown in FIGURE 6.2(a). Each column has a moment of 
inertia, Ic, and the beam moment of inertia is represented by Ip, 

Mc and Mpp are the ultimate moment capacities of the column and beam 
respectively. The lumped model for this frame is shown in FIGURE 


pec (b).. I B and I,. are the moments of inertia of the lumped beam and 


Y Ke 
column in the simplified structure. Moic and Mor are the moment 
capacities of the lumped column and beam respectively. Both the original 
and the lumped structures were analyzed by slope-deflection procedure 
for a horizontal load, H, applied at the top of the frame. For the two 


frames to have identical deformations the following relationships must 


be satisfied: 


& BAe 
lic 21, (6.1) 

e A se 
lip 21, (Cr) 
M = 2M oR ee) 
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In the model shown in FIGURE 6.2(b) the length of the beam 
is Lp/2. Alternately, it is possible to have the length of the beam 
in the lumped model equal to that in the real frame provided both 
ends of the beam are rotated through the same angle. This later 
procedure has been followed in the analysis in CHAPTER III to avoid 


confusion in the input data. 


6.2.1 Lumping of Multibay Frames under the Action of Lateral Load 
Only 


For a multibay frame having equal bay lengths and 
Subjected only to a lateral load, the following relationships 
Similar to equations (6.1) to (6.4) must be satisfied for the 


lumped model shown in FIGURE 6.2(b) 


le = ZI, (6.5) 
lip Mig tests (6.6) 
MLC Rbincr et Pe (6.7) 
Mote = 2 Mp, (6.8) 


where I, and I, represent the moments of inertias of 
columns and beams respectively, whereas Moc and Mop are the ulti- 
mate moment capacities of the columns and beams respectively. In 
the case of unequal bay lengths, the beam in the lumped structure 


can be assigned a length equal to the average length of all the beams 
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in the original structure. The moment of inertia of the lumped 


beam is given by: 
Ip = 2 Lip 2 Z(1,/Lp) pa. - (6.9) 
where Lip is the average length of the beam. 


In the event of unequal plastic moment capacities at 
the two ends of the beam in a story, the plastic moment capacity 
of the lumped beam is given by: 


M = =(M + M 


PLB PBI cea deed 138 


pp2) 


where M 1 and Mops are the plastic moment capacities 


PB 
of the left and right ends, respectively, of each beam in the story. 


6.2.1.1 Check of Basic Lumping Procedure 


To justify the above procedure, several one story 
Structures, one bay to seven bays in width were analyzed. 
The results of the approximate analysis developed in CHAPTER 
III were compared with those of the more rigorous analysis 
developed in reference (C6). In each case the story height 
was 11 feet and the bay width was 18 feet. The beams were 
of reinforced concrete, 10 inches x 20 inches in cross- 
section with 2 percent steel; the columns were 10 inches x 
10 inches with 4 percent reinforcement. Each frame was 
analyzed under the action of a horizontal load applied at 


the top of the frame. A typical comparison of the two 
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analyses is given in FIGURE 6.3 for a four bay, one story 
frame. The dashed line shows the results given by the 
rigorous analysis developed in reference (C6) and the solid 
line shows the results of the approximate analysis. The 
initial slopes of the load-deflection curves are almost 
identical. The difference in the two curves becomes notice- 
able as inelastic action of the structure becomes severe. 
The sequence of hinge formation as predicted by the two analyses 
is shown in FIGURE 6.3 and explains this difference. The 
rigorous analysis predicts that the first hinges will form 
at the bases of the three interior columns, whereas the 
second hinges, at the bases of the exterior column form 
later in the loading history. By lumping these columns, 

the approximate analysis predicts that all hinges form 
together, that is, at the base of the lumped column; causing 
the change in the hinging load. However, the difference is 
not significant. The ultimate load for this structure was 
75.2 kips by the approximate analysis whereas the rigorous 
analysis predicts an ultimate load of 74.4 kips. Thus for 
this frame, the approximate analysis overestimated the 
ultimate load by about one percent. All other frames 
analyzed in this section, exhibited similar trends in load- 
deflection curves and the sequence of formation of hinges. 
TABLE 6.1 compares the ultimate loads predicted by the two 
analyses, for all the frames. In all cases the approximate 
analysis tended to overestimate the ultimate load from 


1.0 to’ 2.2 percent: 
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6.2.1.2 Effect of Variations in Relative Stiffness of Beams 


eens eee teeny 


and Columns 


To investigate the effect of varying the beam to 
column stiffness, the four bay, one story structure, des- 
cribed in section 6.2.1.1 was selected as a model. Seven 
different sets of columns, varying in cross-section from 
10 inches x 10 inches to 34 inches x 34 inches with four 
percent reinforcement, were considered. All frames were 
analyzed under the action of a horizontal load applied at 
the top of the frame, by both methods; that is, the 
approximate analysis developed in CHAPTER III and the rig- 
orous analysis developed in reference (C6). No vertical 
loads were considered. A typical example is shown in 
FIGURE 6.4. In this case the columns were 20 inches square. 
The dashed line gives the results of the analysis developed 
in reference (C6), whereas the solid line is that of the 
approximate analysis. Both analyses predicted similar load- 
deflection curves and hinging patterns. The initial slope 
of the load-deflection curves given by the two analyses 
are almost identical as can be seen in FIGURE 6.4. The 
ultimate load predicted by the rigorous analysis was 557.5 
kips whereas approximate analysis predicted an ultimate 
load of 568.0 kips. Thus the approximate analysis over- 
estimated the ultimate load of this frame by 2 percent. 

All other frames analyzed in this paragraph yielded 
Similar comparisons. FIGURE 6.3 also shows a frame 
in this series. The hinging sequence differed 


for the two frames. TABLE 6.2 compares the ultimate 
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loads for all the frames. The ratio of column to beam 
stiffness varied from 0.24 to 42.06. By varying this 
quantity it was possible to have different hinging sequence 
for different frames in the rigorous analysis. In all cases 
the approximate analysis tended to overestimate the ultimate 
load by 0.5 to 4.5 percent. The errors tended to increase 


as the column stiffnesses increased. 


6.2.1.3 Effect of Unequal Beam Spans 


In the previous sections only those frame having 
equal bay legnths have been considered. However, many 
Structures have unequal bay lengths and subsequently beams 
of different length in the same story. In the approximate 
analysis, however, all the beams in a floor have to be 
lumped together causing the hinges to form simultaneously 
in all the beams of a given floor. To compensate for this 
effect the use of average bay length has been suggested 
above. The moment of inertia of the lumped beam will be 
given by equation (6.9) whereas the ultimate moment 
capacities will be given by equation (6.8), or (6.10) if 


no transverse loads are present on the beams. 


To study the effect of span variation, two of 
the four bay, one story, frames discussed above were con- 
sidered. The first frame had 10 inch square columns and the 
second 20 inch square columns. The beam length was 11 feet 


in one, two or three of the four bays whereas all other bay 
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lengths were 18 feet. Various possible combinations were 
considered, resulting in 11 different frames for each 

column size. All 22 frames were analyzed under the action 

of a horizontal load, applied at the top of the frame, by 

both the approximate and the rigorous methods. No vertical 
loads were applied to the structure. In all cases, both 
analyses predicted similar load-deflection curves and 

hinging patterns for the structures. The initial slope of 

the load deflection curves for the frames were found to be 
almost idential by both methods. Two typical load-deflection 
curves are presented in FIGURES 6.5 and 6.6. In both FIGURES 
6.5 and 6.6 the dashed lines represent the results of the 
rigorous analysis whereas the solid lines represent the 
results obtained by the approximate analysis. FIGURE 6.5 

is representative of the first set of frames, which had 10 
inches square columns. These frames failed in a sway mechan- 
ism. The ultimate load predicted by the rigorous analysis for 
the frame shown in FIGURE 6.5 was found to be 74.5 kips 
whereas the approximate analysis predicted an ultimate load 

of 75.2 kips. Thus the approximate analysis overestimated 

the ultimate load by approximately one percent. FIGURE 6.6 
shows the results of a similar frame, chosen from the second 
set, which had 20 inches square columns. These frames failed 
in a combined mechanism. The ultimate load for the frame 
shown in FIGURE 6.6 was found to be 522 kips by the approximate 
analysis whereas the rigorous analysis predicted an ultimate 
load of 525 kips. In this case the approximate analysis under- 


estimated the load by 0.6 percent. The ultimate loads for frame 
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considered in this section are given in TABLE 6.3. In all 


cases the error was found to be within + 1.4 percent. 


6.2.2 Lumping of Multibay Frames Subjected to Combined Loads 


In section 6.2.1 the lumping procedure for frames sub- 
jected to lateral load, has been derived. However, the structure 
will carry vertical load in addition to lateral load. The vertical 
load in a structure could be applied as concentrated loads on the 
tops of columns or as the uniformly distributed loads on the beams. 
It was assumed in CHAPTER III, that the effect of the total vertical 
load in a story could be simulated by placing concentrated loads 


on the column tops. 


6.2.2.1 Effect of Column Top Loads 


To study the effect of column top loads, the four 
bay, one story frame having 20 inch square columns, shown 
in FIGURE 6.4, was again analyzed. The ultimate axial 
load capacity, AP: for this column is 2110 kips. Column 
top loads, P, varying from 25 kips to 2000 kips, producing 
8 different loading conditions, were considered. The frames 
were analyzed by both the rigorous method and by the approxi- 
mate method. The frames were lumped by the procedure developed 
in section 6.2.1 for the purpose of the approximate analyses. 
All frames were analyzed under the action of incremental 
horizontal load applied at the top of the frame, while the 
vertical load remained constant. Typical load-deflection curves 


from the two analyses are shown in FIGURE 6.7. In this case 
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the load on the top of each column was 500 kips, giving 

Ee = 0.237. The dashed line in FIGURE 6.7 represents 

the results obtained from the rigorous analysis whereas 

the solid line represents the results of the approximate 
analysis. The initial slopes of the load deflection 

curves obtained by the two analyses are almost identical. 
Similar hinging patterns were also predicted. The ultimate 
lateral load predicted by the rigorous analysis was 633 kips, 
whereas the approximate analysis predicted an ultimate 

load of 642 kips. Thus the approximate analysis overestimated 
the ultimate load by about 1.5 percent. All other frames 
exhibited similar load-deflection curves and similar hinging 
patterns were obtained by the two analyses. TABLE 6.4 
compares the ultimate loads for all frames. The error was 
about 2 percent for ages less than 0.25, increasing 
exponentially to 11 percent for Pike = 0.95. These errors 
are probably due to the fact that the effect of axial 

load on the coefficient C and S commonly used in slope- 
deflections equations are ignored by the approximate 
analysis. This problem has also been discussed in section 


Sea 
6.2.2.2 Effect of Uniformly Distributed Loads on Beams 


To study the effect of applying a uniformly dis- 
tributed load to the beams, the four bay, one story structure 
was analyzed under the influence of incremental lateral load 


and a constant uniformly distributed vertical load on the 
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beams. The lumped models, derived using the principles dev- 
eloped in section 6.2.1, were analyzed by approximate method 
Whereas the corresponding actual frames were analyzed by the 
rigorous analysis. The results show that the approximate 
analysis tended to grossly overestimate the ultimate load 
depending upon the magnitude of uniformly distributed load. 
Accordingly, it is necessary to reduce the ultimate moment 


capacity of the beams in the approximate analysis. The 


following analysis was performed to find the reduction factor 


which can be applied to the ultimate moment capacities 
of the beam. In such cases the equation (6.8) and (6.10) 
can be redefined as: 


Moi p = 220. Mp, Pe eieeg (se 


or M = Xe! (Mop + Mop) ys. « (612) 


where ais the reduction factor. 


A sway mechanism will require two beam hinges plus 
hinges at the bases of the columns. If no load is applied 
between the ends of a restrained beam in a laterally loaded 
frame, and if the point of contraflexure is at mid-span of 
the beam, the moment diagram for the beam will be as shown 
in FIGURE 6.8(a) and the beam plastic hinges will occur 
simultaneously as shown in FIGURE 6.8(b). When gravity 
loads are applied between the ends of the beam, the 


negative moments due to the vertical load add to the wind 
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moments at the leeward end of the beam, as shown in FIGURE 
6.8(c) reducing the lateral load required to cause hinging 
at that point; and subtract from the wind moments at the 
other end tending to develop a hinge configuration as 


shown in FIGURE 6.8(d). 


Consider a restrained beam of span, L, and ulti- 
mate moment capacity My» in both positive and negative 
moment regions, subjected to a constant uniformly distri- 
buted load, w per unit length, and a wind moment distribution 
characterized by the end moments, My As the wind moment, 
Mv is increased, the first hinge forms at the leeward 
end of the beam; on further increase in the wind moment, a 
second hinge will form either at the windward end or in 
between the ends of the beam. The results istdtited from 
such an analysis are plotted in FIGURE 6.9 for various ratios 
y of midspan moment capacity to end moment capacity. The 
term y will be explained in the next paragraph. The ratio 
B= MyM is plotted on the vertical axis whereas the values 
of wl? 7M, are plotted on the horizontal axis. For a known 
value of wl? 7M, the ultimate moment capacity of the beam should 
be multiplied by the corresponding value of 8 obtained 
from FIGURE 6.9 to get the correct mechanism load. If the 
value of 8B is taken corresponding to the first hinge 
configuration, the result obtained by the approximate 
analysis of the frame will be too conservative since the 
Second hinge is forced to form at the first hinging 


load. If the higer value of 8, given by the 
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second curve is used, the result will be unconservative since 
the first hinging load is increased. Therefore, it is 
Suggested that the quantity a, in equation (6.10) should lie 
close to the mean of the two values of 8 obtained from 


FIGURE 6.9. The mean value of 8 is plotted in FIGURE 6.10. 


In reinforced concrete structure, however, the 
positive moment capacity is normally less than the negative 
moment capacity. The extreme case would correspond to the 
elastic gravity load moments which for a fixed ended beam 
lead to a mid-span plastic moment capacity equal to half the 
plastic moment capacity at the ends. Let the plastic moment 
Capacity in the positive moment region is given by yM. where 


p 
M. is the plastic moment capacity in the negative moment 


p 

region. Normally the value of Y will vary from 0.5 to 1.0. 
The results of similar analysis described in previous para- 
graphs for various values of Y varying from 0.5 to 1.0 are 
shown in FIGURE 6.9. The mean value of 8 obtained from 
FIGURE 6.9 are plotted in FIGURE 6.10. From the FIGURE 

6.10 it can be seen that the value of a can be approximated 
in straight line relationship with WLM, and a sufficiently 


accurate approximation to the values of a will be given by 


the equation: 


Hee eOeeO 01S CI-at) (WL 7M, ) on Gas 


For the case of beams with equal positive and nega- 
tive My S: that is, Y = 1.0, the equation (6.13) will 


result in: 
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and for the case of beams having half the value 
of negative M53 in the positive moment region the equation 


(6.13) will result in: 
eee TAO = 0.0825 (WL 7M, ) eee 6 ae) 


The other case of the reinforced concrete struc- 
ture will be the case of a beam which will have different 
moment capacities at the two ends, !he leeward end of the 
beam will have greater value of moment capacity than the 
windward end. If n is the ratio of the windward end 
moment capacity, MP» to the leeward end moment capacity 
the value of a given by equation (6.13) will be sufficiently 
accurate as long as the hinge does not form at the wind- 
ward end of the beam. For most of the structure this is 


true and this condition is satisfied if: 


aul Soucty.< in Bok ilGralf) 


If the value of a obtained by equation (6.13) 
does not satisfy the equation (6.16), a conservative esti- 
mate can be obtained by taking a equal to 0 or n depending 
upon whether the value of a obtained from equation (6.13) 


is less than (1-n)/2 or greater than n. 


In the analysis of a large building structure it 
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would become tedious if a has to be computed for every 

beam in the structure. Accordingly, it is desirable to 
have one general value of a which could be used at least 

in the preliminary analyses of multi-story structures. 

Such a factor can be derived from FIGURE 6.10. In the 
design of reinforced concrete buildings by the 1970 ACI 
code proposed, it will be necessary to design the structure 


for ultimate loads of: 


Ce O33) Oe ae ed Ree G 17) 
or = O75 C1 4D eel, TE eee) P62 18) 
where U = required ultimate load capacity of 
section. 
DL = dead load 
LL = live load 
and WL = wind load 


A well designed shear wall-frame structure will 
come very close to satisfying both equations (6.17) and 
(6.18) simultaneously. Normally the beams will be designed 
for the moment capacity, Mo? equal to wh 2/10 to satisfy 
the equation (6.17). At the ultimate state under lateral 
loads the beams will carry a gravity load of 0.75(1.4DL + 
1.7LL) or 75 percent of the design gravity load. For this 
load the value of WL? 7M, is 7.5 and the corresponding value 


of a, given by FIGURE 6.10, should vary from 0.62 for the 


case of equal positive and negative plastic moment capacities 
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to 0.44 for the case of beams having half the value of 


end moment capacity in the positive moment region. 


The correction factor a will lead to a better 
prediction of the load at which the beam mechanisms 
form in the structure and hence a better prediction of the 
load at which significant softening of the structure begins. 
On the other hand, however, the work done in deforming the 
beam will be underestimated considerably with a resulting 
decrease in the stiffness predicted for a structure con- 
taining beam hinges. This is illustrated in FIGURE 6.11. 
The work done in deforming the columns in mechanism (a), 


Shown in FIGURE 6.11, through an angle 6 is: 


a a _ wb.b.a a 
Mo (e+ Fo) + M, (6+ 7s ce bilan 7H Oo + Wa. 5 @ + HyhO 
ae M(4Y De " war + Hyh ae Sen.) 


Combining equation (6.19) and (6.13) and solving for H, 


will result in: 
3 (1-a)a/L 
Hyh = Mo LCT ay 5 - 0.03(7-37 } ve UG 220) 


The work done in deforming the columns in mechanism 


(b), shown in FIGURE 6.11, through an angle 6 is: 


2am ,° Hohe 


= M sido LO 
or Hoh 2a : ( ) 
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The ratio, Ho/H, » will be given from equations (6.20) 
and (6.21) as 


2 "Tae 
ai = vara F -O a/L Manne Once) 
Hy NT 0st) 


For the value of y equal to one and suggested 
value of a equal to 0.62 the ratio Ho/H, given by 
equation (6.22) will be 0.785. In other words, approximate 
analysis will underestimate the ultimate horizontal load 
for the structure shown in FIGURE 6.11 by 21.5 percent. 
However, this error will be greatly reduced for the multi- 
Story structure, since many of the beams will not contain 
hinges at the ultimate load. Thus it will be seen in 
example presented later in this CHAPTER, for large structures 
the load at the first major decrease in stiffness is pre- 
dicted fairly well by the approximate analysis but the 


stiffness, after hinging,is underestimated by it. 


To study the effect of a uniformly distributed load 
on the beam, the four bay, one story model analyzed earlier, 
was considered. The ultimate load capacity BAe for each 
column was 2110 kips whereas each beam required a uniformly 
distributed load, Ws of 11.44 kips per foot to form hinges 
at the ends under the vertical load only. Four different 
uniformly distributed loads varying from 0 to 76 percent of 
Wh were chosen. In each case, however, the total vertical 
load was 2500 kips corresponding to a ratio of B/ Pg of 


approximately 0.24 for each column. All frames were analyzed 
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by both the rigorous and approximate methods. The values 
of a used in the approximate analysis were taken from 
FIGURE 6.10. Vertical loads remained constant while 

the horizontal load applied at the top of the frame was 
incremented in both the analyses. A typical example, in 
which the uniformly distributed load was taken as 76 
percent of We is shown in FIGURE 6.12. The initial slope 
of the load-deflection curves obtained by the two analyses 
are almost identical. The first hinge in the rigorous 
analysis was detected earlier than in the approximate 
analysis causing a discrepancy in the load-deflection curve 
at this point. However, the approximate analysis predicts 
a more rapid softening of the structure following the 
formation of the first hinge producing approximately the 
Same ultimate load as that predicted by the rigorous analysis. 
In this case the ultimate load predicted by the rigorous 
analysis was 544 kips, whereas the approximate analysis 
predicted an ultimate load of 522 kips. Thus the approxi- 
mate analysis underestimated the load by 4 percent. TABLE 
6.5 compares the ultimate loads obtained by both analyses 
for all the frames. In all cases the approximate analysis 
under-estimated the ultimate load by less than 4 percent 
except for the first case where no uniformly distributed 
load was present. In this case the approximate analysis 


overestimates the ultimate load by about 1.5 percent. 


The example frame considered above was reanalyzed 
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with the beam length of bays 1 and 3 set at 11 feet. The 
results are shown in FIGURE 6.13. A comparison of FIGURES 
6.12 and 6.13 shows that for the frame having varying 
span lengths, the load-deflection relationship predicted 
by the approximate method agreed more closely with the 
results predicted by the rigorous analysis. The ulti- 
mate load predicted by the rigorous analysis was 584 kips 
whereas the approximate analysis predicted an ultimate 
load of 564 kips, underestimating the ultimate load by 


3.4 percent. 


6.2.3 Multi-story Structures 


In an attempt to check the validity of the lumping pro- 
cedure, derived in section 6.2.1 and 6.2.2, for multi-story 
frames linked to shear wall, several frames were analyzed by the 
program described in CHAPTER III and the results were compared 
with those obtained from the program developed in reference (C6). 
The member sizes of the basic four bay frame, used as an example 
in section 6.2.1.2, were used in frames of one, two, three and 
four stories. These frames were linked with a shear wall, 10 
inches x 120 inches in cross-section, having one percent rein- 
forcement. Hinged link beams were used to restrict the study 
to effects of the assumptions made in lumping of the frame itself. 
The frames were analyzed under the action of horizontal loads 
applied at each floor level. The roof level had one half the 
lateral load applied at each floor. No vertical loads acted on 


the structure. A typical three story structure is shown in FIGURE 
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6.14(a) and the lumped equivalent structure is shown in 

FIGURE 6.14(b). The four frames were analyzed by the rigorous 
method and by the approximate method and the results were 
compared. In all cases the initial slopes of the load- 
deflection curves were essentially equal. A typical load- 
deflection curve for the three story building is shown in 
FIGURE 6.15. The dashed line represent the results obtained 
from the rigorous analysis whereas the solid line is that of 
the approximate analysis. Both analyses predicted similar 
load-deflection responses and hinging patterns. Due to the 
Slight differences in the predicted hinging pattern the two 
curves do not match exactly over the inelastic portions of the 
curves. The ultimate load for this frame was found to be 263 
kips by the rigorous analysis and 260 kips by the approximate 
analysis. Thus the approximate analysis underestimated the 
ultimate load by 1 percent. The comparisons of the ultimate 


loads for all four frames are presented in TABLE 6.6. 


A 20 story, two bay structure, shown in FIGURE 6.16 
(a) was chosen as a second example and will be termed as 
frame "A". The structure consisted of a two bay rigid frame 
linked with a 10 inches x 60 inches wall by hinged link beams. 
The sizes of the frame members are taken from reference (C6). 
The equivalent lumped structure is shown in FIGURE 6.16(b). 
No gravity loads were considered except for the self weight 
and the value of a was taken as one. This structure was 
analyzed under the action of incremental loads applied at the 


floor levels. The load-deflection curves obtained from the 
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rigorous and approximate analyses are shown in FIGURE 6.17. 

The dashed line shows the results obtained from the rigorous 
analysis whereas the solid line is that of the approximate 
analysis. The hinge configuration of this structure is shown 
in FIGURE 6.18. The initial slopes of the load-deflection 
curves are identical and the hinging patterns are also similar, 
as can be seen in FIGURES 6.17 and 6.18. The load factor 
predicted by the rigorous analysis is 2.7 whereas the approxi- 
mate analysis predicted a load factor of 2.65. Thus the 


approximate analysis underestimated the load by about 2 percent. 


The above 20 story structure was also analyzed under 
the action of incremental horizontal loads at each floor level 
and constant vertical loads. This frame will be termed as 
frame "B". Uniformly distributed loads of 4 kips per foot 
were considered on all beams except the roof girders, which 
were subjected to loads of 2.2 kips per foot. The value of a 
was computed from FIGURE 6.10 for each beam. The load-deflection 
curves obtained from both analyses are presented in FIGURE 6.19. 
The dashed line was obtained by the rigorous analysis whereas 
the solid line represents the approximate analysis. Due to the 
differences in the hinging patterns predicted by the two 
analyses, the curves differ slightly in the inelastic range. 
hinging patterns obtained are shown in FIGURE 6.20. The approxi- 
mate analysis predicts hinges at the ends of the beam whereas the 
rigorous analysis predicts one hinge at the leeward end of the beam 
and-one between the two ends. The approximate analysis predicts 


earlier failure of the structure. 
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This frame was also analyzed with various wall sizes. 
The load-deflection curves obtained by both analyses for a 10 
inches x 80 inches wall, frame "C", and for a 10 inches x 96 
inches wall, frame "D", are shown in FIGURES 6.21 and 6.22. 
In these cases similar comparisons between the two predicted 
load-deflection curves and hinging patterns were obtained. In 
the case of frame "D" the ratio of the wall stiffness to the 
sum of column stiffness varied from 4.87 in the bottom story 
to 300 in the top story. The choice of this ratio was limited 
in the sense that the analysis presented in reference (C6) 
does not converge for large differences of wall to column 


stiffness. 


It can be concluded from the results of the frames 
A,B, C and D analyzed above, that the approximate analysis 
yields reasonable description of load-deflection curve for 
the multi-story structure when there is no vertical load present 
on the structure. The approximate analysis yield a slightly 
inaccurate but conservative description of load-deflection curve 
for the structure when the gravity loads are present on the 
structure. This is due to the way the approximate analysis 
takes into accountthe uniformly distributed load on the structure. 
As mentioned earlier the value of a is so chosen that the 
approximate analysis will always result in underestimation of 
the ultimate load. The designer may wish to use a different 
value of a to more closely predict some other part of the load- 
deflection curve such as the start of hinging. However, the 


discrepancy in the results obtained by approximate and rigorous 
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analysis are within tolerable limits as can be seen by the 


results of frames A, B, C and D. 


6.2.4 Lumping of Link Beams 


In deriving the analysis it was assumed that the 
beams linking the shear wall to the rest of the structure (link- 
beams) in each story are assumed to be represented by a Single 
link beam as shown in FIGURE 3.1. The lumped link beam has a 
length equal to the average length of all the link beams in 
the story except that the lumped link beams must have the same 
length in all stories. The stiffness, EI/L of the lumped link 
beams shall be taken as sum of the stiffnesses of all the link 
beam in the story in question. The plastic moment capacity of 
the lumped beam shall be taken as one, for all the corresponding 


beams. 


The assumptions concerning the use of the average 
beam length and the derivation of the EI and Mo terms used in 
the analysis have been discussed in sections 6.2.1, 6.2.1.3 
and 6.2.2.2. These explanations will also apply to the link 
beams. The general value of a derived in section 6.2.2.2 will 
often be too low since the major portion of the moments in 


the link beams often result from wind load rather than gravity 


load. 


6.2.5 Lumping of Shear Walls 


In deriving the analysis it was assumed that the shear 
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walls in one story can be represented by a Single wall with a 
moment-curvature relationship defined by superimposing the 
moment curvature relationships of each of the individual walls 
in that story. The resulting moment-curvature relationship 

is approximated by a bi-linear curve. The slope of the second 
branch of this curve may be small but not zero and the ratio 
of the slope of the second branch to that of the elastic 
portion of the curve must be the same in all stories. The 
neutral axis of the wall is assumed to coincide with the 


centroid of the uncracked wall. 


It is customary to assume that the wind shear resisted 
by the shear walls in a structure can be distributed between 
the various shear walls in each story in proportion to their 


(K4) has shown that this assumption 


individual stiffnesses. Khan 
will result in the correct moments and shears in the shear wall 


only if one of the following conditions is satisfied: 


1. Each shear wall has constant section properties 


throughout the height of the building. 


2. Where wall sections change, the relative stiff- 
ness of each wall remains unchanged throughout 


the height of the building. 


(K4) (T1) 


In the examples presented by Khan and Tezcan 
the errors in the wall shears and moments estimated by distri- 
buting the total shear in the walls in any story in the ratio 


of the stiffnesses of the walls ranged from about -300 percent 
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to about +200 percent. The errors appeared to be smaller when 
the wall moments were distributed in this way rather than the 


wall shears. 


Tezcan has shown, however, that the lateral deflection 
of a number of shear walls connected by hinged link beams is 
equal to the deflection of a single shear wall having a stiffness 
in each story equal to the sum of the individual wall stiffnesses 
in that story. Thus, the lateral deflections computed for a 
structure by lumping the wall stiffness as is done in this 
analysis will be essentially correct. To compute the shears 
and moments in the individual walls it is then necessary to force 


each individual wall through the computed deflections. 


The effect of summing the EI/L values of the shear 
walls was studied using a ten story building containing three 
Shear walls, inter-connected by linked beams at each floor 
level as shown in FIGURE 6.23(a). This structure was also 
discussed in reference (K4). A first order analysis of the 
actual structure was performed using a matrix method. In 
addition, the three walls were lumped together into the single 
cantilever as shown in FIGURE 6.23(b) which was also analyzed. 
In both cases analysis were limited to elastic range. The 
deflections at each floor level as predicted by the two 
analyses are compared in TABLE 6.7. The actual shears obtained 
by distributing the total shear in the ratio of the respective 
wall stiffnesses are presented in TABLE 6.8. An incorrect 


distribution of shears resulted from this procedure. Since the 
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lumped model yields almost correct deflections, the shears and 


moments can be computed from the deflected shape of the structure. 


The final assumption considered in lumping the shear 
wall was that the neutral axis of the wall was at mid-depth 
of the wall. Although the exact location of the neutral axis 
depends on the shape of the wall, the distribution of the wall 
reinforcement and the axial load and moment in the wall, it will 
generally lie between the mid-depth and the compression face 
of the wall. This assumption affects the ultimate load in two 
ways: 1. The moments induced in the link beams differ on 
the two sides of the core as shown in FIGURE 6.24. This effect 
tends to cause hinging to start earlier than predicted on the 
windward side and later on the leeward side of the wall. Since 
the link beams on both sides of the wall are lumped together 
into one link beam, this assumption has relatively little 
effect on the overall behaviour predicted by the analysis. 
2. If the two link beams are similar, more work will be done 
in raising the wall end of the windward beam than is dissipated 
by lowering the wall end of the leeward beam. This is not true 
if the neutral axis of the wall is at mid-depth. As a result 


the analysis will tend to underestimate the stability of the 


structure. 


6.3 Analysis of Planar Structure 


It is assumed in the analysis that each floor acts as a rigid 


diaphragm and that the structure does not twist about its longitudinal 
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axis. 


Benjamin’ oe 


has shown that a frame analysis based on the 
assumption that the floor diaphragm are rigid, is satisfactory if 


the ratio of the diaphragm to frame rigidity is greater than 1.0. 


The problem is also discussed by Goldberg ‘&2) 


who presents 
an exact first order analysis including shear and flexural deforma- 
tions of the floor diaphragm and the walls. It is interesting to 
see the comparison of elastic deflections and wall shears at each 


floor level of ten and twenty story building taken from reference 


(G2). The comparison has been presented in section 6.4. 


In an example presented in reference (W6), moving the shear 
wall from the middle of a 60 foot wide building to 7.5 and 15 feet 
from the middle led to reductions in deflection of the top of the 
wall by 4 and 15 percent, respectively. Generally, however, the 
non-symmetrical nature of the building tends to increase the maximum 


deflections of the most highly deflected frame. 


6.4 Comparison with other Published Examples 


The lumping procedure developed in this CHAPTER was checked 
in part by comparisons of the computed deflections with those of 
several published example structures. The example structures selected 
were reduced to their equivalent models and analyzed by the program 


developed in CHAPTER III. 


Example one is a ten story, three bay steel frame first 


i iy Cai ae 

ket 

ae Gas) 

a EDD eHb AY Noe Gime 

oh ae coy ee 

ie om fo baesd ef2avlfans omer? 6 jerd ni ie ‘dl ee ra 

h | ai Vos ostertse net sotpin a6 mosydgstb yoo fr’ oa sedt m 101 

bie Out nans Ed 84D at xbPbtpty smsy?. ot ‘ments a Ao 

Ae Maer 2 

a etnsesrg onw (58) nadine yd iaacuaett oot ar mt adt Ms 

* -emoteb {svxoft brs vesde ontbutadt ateyfsnstebIm, eat 21 

ties os onttesradnt et +1 .eblsw ads bag, mpsiiqetb watt “2 

vs dos +6 evesie iw bis 2nottoat'teb: otsesle to net 
390913T9Y MONT PUY ontbitud wrote yioswd bas, ‘oaad wt 

i Bs 9 Nortose ut botraesyg ged ead - noe ’ 

mi | ysane ont pnrvem - (obi) goneyetay nt pono spans i’ ; 

i Pot) #887 Gl bas 2.0 ot orth t tud aarw s90% Od 6 aie abe) Ad 1 m0 

i We ony to god sid to nofsaelteb nt enotdauben bad bat of | 

a ant .tevewor .vIlsysisd faut sagen re 


mul Bh say sessyont oF ebnad onibt fad ons %0. ie 
omen? bed et ior te 


hs i. ; ‘ Py 4 ee - a % i 
7 ay p satieanclh att uns dus TNO Paw 
; —- — = es 


nih 
A Pp vi ny : 


\ - 
a 3 Pe dau r my 
. Be } r ay 


ae on sara atnd: af + Sogo au swoon 
‘a nw | to ezon>. att rw 2notsoai tab b re ibs 


dy " ah asnudoinse iat Ting 


be [teal 
ae 
| 
~ 


105 


presented in reference (P2). This structure had unequal bay lengths 
and unequal stiffness in the four column lines. The actual struc- 
ture and the lumped model are shown in FIGURE 6.25. The loading 
and member sizes used for this frame are given in reference (P2). 
All the beams had uniformly distributed load and the corresponding 
values of a were computed from FIGURE 6.10. FIGURE 6.26 compares 
the load deflection curve obtained by the approximate analysis 

to that presented in reference (P2). The dashed line 

shows the curve presented in Reference (P2) and the 

Solid line is that obtained by the approximate analysis. Good 
agreement was found between the initial part of the two curves. The 
discrepancies in the inelastic range of the curve are due to the 
different hinging patterns predicted by the two analyses. Later 
softening but earlier failure was predicted by the approximate 


analysis. 


Example two is a ten story building taken from reference 
(G2). The structure is shown in FIGURE 6.27(a) whereas the lumped 
model is shown in FIGURE 6.27(b). No axial load was considered. 
Comparison of elastic deflections at each floor level and the wall 
Shears are presented in TABLE 6.9 and 6.10. When only bending de- 
formations are taken into account, the error in the wall deflection 
is within one percent in all floors whereas the error in the shear 
force is two percent or less in all floors except for the first and 
top, where the errors are about 7 and 35 percent respectively. The 
comparison of the deflections obtained by approximate analysis with 


the results of combined bending and shear deformation obtained in 
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reference (G2) are in greater error. The error in the bottom part of 
the building is more than the top part of the building. This is 
because of the fact that wall shear is greater in magnitude at the 
bottom part of the building. However, it should be remembered that 
small numbers have been compared. In terms of the shearing forces 

in the wall, however, the approximate analysis gave 


reasonable results. 


Example three is a twenty-story building taken from refer- 
ence (G2). The structure is shown in FIGURE 6.28(a) and the lumped 
model is shown in FIGURE 6.28(b). No axial load was present on the 
Structure. The comparisons of elastic deflections and wall shears 
at each story level are presented in TABLE 6.11 and 6.12. The 
error in the wall deflection in all stories was within one percent 
in comparison to the bending deformation computed in reference 
(G2) and the error in wall shear was less than 5 percent except for 
the first floor and floors near the top of the structure. In com- 
parison to the combined bending and shear deflections obtained in 
reference (G2), the approximate analysis resulted in greater error. 
This error was greater in the bottom part of the building since the 
bottom part of the building attracted greater wall shear in comparison 
to the top part of the building. Here also small numbers have been 
compared. The reasonable comparisons with Goldberg's bending 
deflection analysis suggests that for most normal buildings it jis 


sufficiently accurate to consider the floor as a rigid diaphragm. 


In example two the ratio of the wall to sum of the column 


stiffnesses varied from 514 at the bottom story to 5310 in the top 
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story. In example three the ratio of wall to sum of column stiff- 
nesses varied from 219 at the bottom story to 788 in the top story. 


The comparison was limited to elastic range only. 


6.5 Discussion of Additional Assumptions 


The final group of assumptions presented in CHAPTER III 
dealt with specified aspects of member and frame behaviour, such 
as there is no out of plane behaviour of the structure and there 


is no axial deformation. 


The out of plane behaviour is not a problem except for 
thin walls. Many designers provide columns at the end of walls 


to prevent this tendency. 


Studies of concrete building designs have shown that 
columns in multi-story shear wall frame structure will seldom 
exceed L/r = 35. According to the long column sections proposed 
for the 1970 ACI Code this length would have no effect on the 
majority of the building columns. Therefore there will be no 


buckling prior to mechanism failure. 


In the braced structure H-50, studied in APPENDIX E of 
reference (C6), the failure loads and deformations were influenced 
very little by axial deformations. However, when differential 
shortening of the columns and walls was severe, extensive hinging 
occurred at loads as low as 90 percent of the service load. The 
hinging caused an increase of 14 percent in the roof sway deflec- 


tion at working loads for this particular example. It is important 
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therefore, for serviceability considerations to investigate 
the possibility of relative axial deformations due to elastic 
shortening corresponding to temperature changes, settlement and 


other causes. 


Axial deformations may also cause "cantilever" 
deflections of the entire structure. Khan and Sbarounis ‘K2) 
have proposed a modification to their method of analysis, by 
which this can be accounted for, at least approximately. A 


Similar procedure could be used in the analysis proposed in 


CHAPTER III. 


6.6 Summary of Chapters IV, V and VI 


The assumptions made in the analysis presented in CHAPTER 


III have been discussed in CHAPTERS IV, V and VI. 


The moment-curvature diagrams presented in CHAPTER IV 
Suggest that for concrete structures a strong column-weak beam 
design is desirable in view of the limited rotational capacity 


of column sections. 


The study of the proposed lumping procedure and the 
other assumptions indicates that the approximate analysis yields 
a reasonable prediction of the load-deflection characteristics of 
the building. The ultimate load predicted by this analysis was 


close to the ultimate load predicted by various other analyses. 


Good agreement was found between the deflections and hinging patterns 
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predicted by the rigorous and approximate analysis. The test 
results presented in CHAPTER IX of this report show that the 
analysis closely predicted the load deflection response and 
the development of hinges, although it overestimated the 


measured ultimate load by about six percent. 
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COMPARISON OF ULTIMATE LOADS IN ONE STORY STRUCTURES, 


ONE TO SEVEN BAYS IN WIDTH 


fuliber Ultimate Load in Kips 


Analysis By Present 
Reference(C6) Analysis 


of Bays 


TABLE Dez 


COMPARISON OF ULTIMATE LOADS IN FOUR BAY, 
ONE STORY STRUCTURES 


ee RAEN Ultimate Load in Kips 
Column Beam 
Size Stiffness 
Inches I6/Le¢ by 


Reference(C6)| Analysis 


Analysis Present 


Percent 
Error 


Percent 
Error 
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TABLE 6.4 


COMPARISON OF ULTIMATE LOADS. IN 
FOUR BAY, ONE STORY STRUCTURES 
CARRYING LOADS AT THE TOP OF COLUMNS 


The Ultimate Load in Kips | Percent 
Ne soe Analysis By Present 
o |Reference(C6)} Analysis Error 
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TABLE: 6.5 


COMPARISON OF ULTIMATE LOADS IN 
FOUR BAY, ONE STORY STRUCTURES 
CARRYING UNIFORMLY DISTRIBUTED LOAD 


Vater Ultimate Load in Kips SehcaRe 


OF Analysis By Present 
1 Reference(C6)} Analysis Error 


U.D. Load 
Load to Cause 
end hinges 


TABLE 6.6 


COMPARISON OF ULTIMATE LOADS IN 
FOUR BAY AND ONE, TWO, THREE AND 
FOUR STORY BUILDINGS 


No. of Ultimate Load in Kips Percent 
Analysis By Present 
Stories | Reference(C6 Analysis Error 
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TABLE 6.9 


COMPARISON OF DEFLECTIONS AND WALL SHEAR 
IN TEN-STORY BUILDING (ONLY BENDING 
DEFORMATIONS INCLUDED IN GOLDBERG'S ANALYSIS) 


Wall Deflections Wall Shear 
Story aa Inches@ AN. VULUEERY (2 AIVILT sie) KIDS. eee 
‘Goldberg's | Present | Percent Goldberg's| Present | Percent 
Analysis | Analysis Error — Analysis Analysis Siperron. 


0. 0. +0.73 a5yai5e | fevdmaia | frvss 
2 0. 0. +0.74 231.869 | 236.955 | +2.19 
3 0. 0. +0.54 202.257 | 201.170 | -0.54 
4 0. 0. +0.39 {7¢.082 | I71-735 | -orz0 
5 0. 0. 40.24 | 142.234 | 141.950 | -0.20 
6 0. 0. +0.06 ir3445° 4 Tisez5s) 0.17 
7 0. 0. 40.01 | 85.338 85 1 
8 0. 0. -0.06 58.033 | 57 48 
9 0. 0. eo 1G 32.644 | 32 

10 0. 0. -0.19 9.654 
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TABLE 6.10 


COMPARISON OF DEFLECTIONS AND WALL SHEAR 
IN TEN-STORY BUILDING (BENDING AND SHEAR DEFORMATIONS 
INCLUDED IN GOLDBERG'S ANALYSIS) 


Wall Deflections, Inches 


Wall Shear, Kips) 


Present Percent 


Analysis Analysis EYYor Analysis Analysis Error 


0.0142 0.00554 -60.98 241.254 270.45 +11,96 
0.0369 0.02035 -44 .85 224.436 236959 = OE OY 
0.0591 0.03720 -37.05 197.465 201.170 12 os 
0.0838 0405/82 -31.60 168.545 17 wss +1389 
0.1100 0.07979 -27.46 139.743 141.950 +. ly 58 
0.1369 0.10386 -24.13 11 le 2e 113 a205 + 1.46 
0.1641 0.12891 -21.44 84.186 85.245 + 26 
0.1910 0.15450 -19.11 57.508 57 #5 - 0.58 
0.2175 0.18032 -17.09 32.658 32.85 - 1.40 
Q. 0.20620 -15.28 10.129 6.280 -38.00 
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TABLE 6.11 


COMPARISON OF DEFLECTIONS AND WALL SHEAR 
IN TWENTY-STORY BUILDING (ONLY BENDING 
DEFORMATIONS INCLUDED IN GOLDBERG'S ANALYSIS) 


Wall Deflections, Inches 


Goldberg's Present 

Analysis Analysis 
0.0156 OF01572 
0.0439 0.04394 
0.0778 0.07782 
0.1188 0.11874 
0.1656 0.16544 
Onc EPl 0.21684 
0.2643 0.26390 
O..3137. 0.31318 
0.3648 0.36418 
0.4171 0.41638 
0.4702 0.46935 
0.5237 0.52271 
0)..5473 O)..57 6.6 
0.6309 0.62946 
0.6842 0.68244 
0.7371 0.73499 
0.7896 0.78705 
0.8417 0.83863 
0.8935 0.88981 
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TABLE 6.12 


COMPARISON OF DEFLECTIONS AND WALL SHEAR IN 
TWENTY-STORY BUILDING (BENDING AND SHEAR DEFORMATIONS 
INCLUDED IN GOLDBERG'S ANALYSIS) 


Wall Shear, Kips | 


Goldberg's Present | Percent 
ERYOr 


Wall Deflections, Inches 


Percent 
Error 


Present 
Analysis 


Goldberg's 
Analysis 


Analysis 


Analysis 


1 0. 0.01572 | -43.85 248.448 284.740 
2 0. 0.04394 80 215.392 225.520 
3 0. 0.07782 15 184.025 190.735 
4 0. 0.11874 60 157.310 159.405 
5 0. 0.16544 44 137.365 141.730 
6 0. 0.21684 121.269 126.010 
7 0. 0.26390 108.240 104.895 
8 0. 0.31318 102.408 103.470 | 
9 0. 0.36418 98.585 | 110.700 
10 0. 0.41368 ; 86.810 87.305 
ih 0.5142 0.46935 | 8. 74.960 79.590 | 
12 0.5686 0.52271 8. 60.905 61.905 
13 0.6229 0.57616 oe 47.433 49.840 
14 0.6767 0.62946 6. 33.933 32.480 | 
15 0.7300 0.68244 6. 23.099 26.685 | 
16 0.7896 0.73499 7 9.889 9.520 | 
17 0.8345 0.78705 5. -2.504 | - 0.945 — 
0.8857 0.83863 | - 5. -15.798 | -16.110 
0.9362 0.88981 4. -27.554 07,605 ts 
0.9865 0.94074 4. -32.641 -57.740 | 
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FIGURE 6.1 


FLOOR PLAN OF A BUILDING 
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FIGURE 6.2 LUMPING OF A PORTAL FRAME 
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FIGURE 6.8 MOMENTS AND HINGE CONFIGURATION IN A RESTRAINED BEAM 
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FIGURE 6.14 THREE STORY STRUCTURE 
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(a) DEFLECTED SHAPE OF THE STRUCTURE 
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(b) MOMENT DIAGRAM 


FIGURE 6.24 EFFECT OF MOVING NEUTRAL AXIS FROM CENTROID 
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(b) LUMPED FRAME 
FIGURE 6.25 TEN STORY STEEL FRAME 
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(b) LUMPED MODEL 
FIGURE 6.28 TWENTY STORY BUILDING WITH SHEAR WALL 
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PART III 


TEST OF A REINFORCED CONCRETE 


SHEAR WALL-FRAME STRUCTURE 
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CHAPTER VII 


TEST SPECIMEN 


7.1 Introduction 


A 20 ft. high, one bay, four story reinforced concrete shear 
wall-frame was tested to failure under a constant vertical load on the 
columns and walls and lateral loads applied in increments, at the 
floor levels. The test was carried out to trace the actual behaviour 
of such a frame and to check on the accuracy of the analysis presented 


POOCHAPTER III. 


The test specimen was constructed and tested inside a steel 
testing frame designed for this purpose. The frame provided the 
reaction for the horizontal loading jacks and supports for observation 
platforms and lateral bracing. Sectional elevations of the testing 
frame, specimen and the loading apparatus are given in FIGURES 7.1 
and 7.2. FIGURE 7.3 shows a photograph of the testing frame, specimen 


and instrumentation during the test. 


7.2 Layout of Specimen and Reinforcement 


The test specimen was a four story, one bay frame, having one 
column and one shear wall interconnected with four beams as shown in 
FIGURE 7.4. Each story was 5 foot high center to center of beams. The 


columns, beams and walls had nominal cross-sections of 6 x 4% in., 
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6 xX 6 in.,and 6x 24 in. respectively. The span of the beams was 84 


inches from the center of the column to the inner face of wall. In 
addition all the beams had a projection of 7 inches from the outer 
face of the columns to ensure adequate anchorage for the beam rein- 
forcement. The various cross-section used for columns and beams and 
for walls are given in FIGURES 7.5(a) and (b) respectively. The 


measured dimensions of the members are listed in TABLE 7.1. 


The reinforcement in the wall and column was welded to the 
respective base plates and extended vertically far enough to be 
anchored past the point where it was no longer required. The fourth 
story column had an additional No. 4 bar in the two faces. The beams 
were reinforced with four bars, two at the top and two at the bottom 
as shown in FIGURE 7.4. Details of the reinforcement in the column 
a wall joints are given in FIGURE 7.6. FIGURE 7.7 shows the clevis 


used for transferring the horizontal jack load to the wall joint. 


The ties and stirrups in the columns and beams were made 
of No. 14 wires, wound in a spiral at 4% inches on centers. Stirrups 
for the first two stories of the wall were made out of U-shaped No. 3 
bars placed at 6 inches on centers. The legs of these stirrups were 
lapped on the width of the wall. In the third and fourth story walls, 
stirrups made of No. 9 wires were provided at 6 inches on centers. 
Two No. 3 bars were used as stirrups in the wall beam joints. The 


column beam joints had stirrups within the joint made of No. 9 wire at 


2 inches on centers. 
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7.3 Construction of Specimen 


The specimen was constructed one story at a time in the vertical 
position at the location where it was to be tested so that it was not 
cracked or stressed due to handling. The reinforcement cages for 
columns and walls were made in the horizontal position and lifted 
into place after they were welded to the respective base plates. The 
concrete was placed in stages. In a typical story the beam concrete 
was placed one day followed by the next higher column and wall the 
next day. The concrete was vibrated using a 5 foot long internal 
vibrator. TABLE 7.2 gives the date of casting and the age of the 
concrete in the various members on the day of testing. FIGURE 7.8 


shows the specimen under construction. 


7.3.1 Formwork 


The formwork consisted of boxes made of 3/8 inch thick 
plastic coated ply-wood. These were braced at several points to 
prevent deflections due to the pressure of the concrete during 
casting. The column and wall forms could be opened at the two 
diagonally opposite corners to break them away from the sides of 
the member. The forms were designed so that the column and wall 
for one story supported the forms for the next higher beam which 
in turn supported the forms for the next higher story. The 
threaded inserts used to hold the rotation meters were bolted 
to the inside of the forms in the correct locations before 
casting. The forms were oiled before placing in position. Transits 
were used to align the formwork vertically. The forms were braced 


against the testing frame to keep them in vertical position during 
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casting. 


Pax. anCuring 


The forms were loosened one day after casting to destroy 
the bond between forms and concrete. Once this had been done 
the formwork was put back in place and left until one week after 
the entire frame was completed to support the forms for higher 
stories and protect the concrete from drying. The tops of the 
beams were covered with wet burlap and polythene. After the 
forms were removed the frame was cured in the laboratory air 
at room temperature for 5] days before the test was carried out. 
The control cylinders made during the various casting stages 


were also cured in the same way. 


Material Properties 


7.4.1 Concrete 


The concrete mix was designed to have an average strength 
of 3000 psi in 28 days. The slump used was 3 inches. The mix 
proportions were: cement 170 lbs., sand 530 Ilbs., coarse aggre- 
gates 370 lbs., water 100 lbs. The coarse aggregate was 3/8 
inch maximum size. This mix was enough for one column and one 
wall plus six control cylinders. For the beams, half of this 
mix was enough to make one beam and six control cylinders. The 
cement used was high early strength portland cement and was 


fresh at the time of use. The concrete was mixed in a 9 cubic 


foot Eirich mixer. 
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Six standard control cylinders were made from each batch. 
Three of these were tested in compression two days before the 
test and one was tested in compression, three days after the 
test. The strength reported in TABLE 7.2 is the average strength 
of the four cylinders from each batch. The remaining two cylinders 
were used for split tension tests whose average strengths are also 


reported in TABLE 7.2. 


7.4.2 Reinforcing Steel 


Intermediate grade deformed reinforcing bar meeting ASTM 
specification A-15 was purchased from a local supplier. All the 
reinforcement of any one size came from the same heat. The heat 
number, composition and standard mill test results of these bars 
were supplied by the producer. Two specimens of each of the No. 3, 
No. 4 and No. 5 bars were tested in tension. The loads and the 
extension on a 2 inch gauge length were recorded by electrical drum 
type recorder mounted on the testing machine. The properties of 
the reinforcement are tabulated in TABLE 7.3. Typical stress- 
strain curves for the bars tested are shown in FIGURE 7.9. The 
yield strength of No. 6 bar was taken from the mill test. The 
cross sectional areas of all the bars were taken from the mill 
test. The value of modulus of elasticity, E, of these bars were 


derived from the stress-strain curves shown in FIGURE 7.9. 


7.5 Base Connection 


The base connections used for the bases of the column and 


wall are shown in FIGURES 7.10 and 7.11 respectively. The column 
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reinforcement extended through holes in the base plate to the 
bottom of the plate. The holes were enlarged at the bottom of 
the plate and the bars were welded in these holes. The base 
plate for the wall was a built up channel and the wall reinforce- 
ment was welded on the inside face of the flanges of the channel. 
In each case the welding was adequate to develop the strength 


of the bars. 


Each base plate was connected to the floor of the laboratory 
by means of two 14; inch diameter high strength steel threaded 
rods, 24 inches apart. The rods had a guaranteed yield strength 
of 128,000 psi and ultimate strength of 142,000 psi. The bolts 
extended through the test floor of the laboratory and were an- 


chored by nuts bearing on the underside of the floor. 


The construction of the wall and placing of curvature 
meters on the column limited the extension of these bolts above 
the floor. As a result the bolts used were too long and a sleeve, 
54; inches long, of 2 inch standard pipe was used as a spacer under 
the nuts. Although the bolts used were high strength steel and 
were adequate to take the loads anticipated at the base, the pipe 
Sleeve in the wall base connection yielded at about half the 
anticipated failure load causing a premature failure of the frame. 
An identical piece of pipe, cut from the same pipe, was tested in 
compression in the laboratory to find the stress-strain characteristics 
of the pipe. From this information a moment curvature relationship 
was established for the base of the wall as discussed in section 


8.10. This predicted moment curvature relationship was used in the 
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theoretical analysis and good agreement was obtained between 


the measured and computed behaviour. 


7.6 Loading Apparatus 


The loading devices were shown in FIGURES 7.1, 7.2 and 
7.3. Separate hydraulic systems were used for the vertical and 
horizontal loads. For each system oi] was pumped from a central 
supply by a pump operated by air pressure. The pumps automatically 
maintained a given oi] pressure. Bleeder valves were used to maintain 


the desired load. 


Horizontal loads were applied by tension jacks connected 
to the frame by the clevises shown in FIGURE 7.7. The jacks reacted 
against the testing frame. The three lower jacks were connected on 
the same hydraulic line so that equal loads could be applied by all 
three jacks. The fourth floor jack was connected to another hydraulic 
line. Hydraulic pressure gauges with a range of 0 to 10,000 psi were 
used to read the pressure between the pump manifold and jack and to 


facilitate the controlling of loads. 


The apparatus used to apply vertical loads on the top of 
the columns and walls is shown in FIGURES 7.1 and 7.2. A loading beam, 
placed parallel to the beams in the frame, applied loads through rollers 
and plates to tops of the columns and walls. The rollers were positioned 
over the center of the wall and the column. A cross-beam, mounted on 
top of the center of the loading beam, was connected in turn to two tie 
rods of 3 x 1 inch cross-section. The lower ends of the tie rods were 


connected to two 100 ton hydraulic jacks acting in tension. These jacks 
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reacted, in turn, against gravity load simulator of the type developed 
at Lehigh University! !), The gravity load simulator is a mechanism 
which can transmit vertical load but cannot resist horizontal load. 
Thus the line of action of the applied load remains vertical even if 
the structure sways. This was checked during the test using scales 


mounted on the simulator and loading beams read with a transit. 


7.7 Lateral Braces 


Five lateral braces were used to prevent the out of plane 
movement of the test specimen. Three were mounted on rods embedded 
in the projecting ends of the top three beams. Two more lateral 
braces were used to prevent out of plane movement of the loading 
beam at the top of the structure. The lateral braces were similar 


(YT). The lateral braces 


to those developed at Lehigh University 
consist of a mechanism which permits the specimen to deflect in its 
plane but does not allow lateral movement of the braced points. 

FIGURE 7.12 shows the details of one half of a lateral brace. The 
measurement of lateral movement onthe fourth story wall indicated that 


no significant out of plane movements of the test specimen occurred 


during the test. 


7.8 Instrumentation and Measurements 


a EERE en 


7.8.1 Introduction 


The objective of this test was to observe the actual 
behaviour of the structure and try to relate it with the analysis 
described in CHAPTER III. To aid in this correlation, the 


following measurements were taken: 
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1. The applied vertical and horizontal loads were 


measured. 


2. The concrete strains were measured at various 
positions along the length of the wall so that 


curvatures and moments could be computed. 


3. Curvature were measured at a number of points 
along each of the columns and beams to enable 


the computation of moments. 


4. The rotations were measured near the ends of the 


member. 


5. Horizontal deflections were measured at each floor 
level as were the vertical deflections of the top 
of the column and wall. Lateral deflections were 


measured near the top of the wall. 


The techniques adopted for the above measurements are 


described briefly in the following sections. 


7.8.2 Measurements of Loads 


The horizontal and vertical loads were measured both by 


oil pressure and by strains measured in previously calibrated bars. 


The graduations on the oil pressure gauges, for horizontal jacks, 
corresponded to approximately 350 lbs. of load per division. It 
was possible to estimate the dial reading to the nearest third of 
a division which corresponded to about 115 lbs. of load. Wherever 


possible the loads were chosen to fall on a division mark. -Between 
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the end of the jack ram and the device used to transmit the 
loads to specimen was a 2 inch long section, one inch in 
diameter on which strain gauges were mounted to measure the 
jack load. Two longitudinal and two transverse strain gauges 
were mounted on the circular piece to form the four arms of 
wheat-stone bridge circuit. In this measurement the effect 
of any bending or temperature change were compensated for (B7) 
To measure the vertical load, one strain gauge was 
mounted longitudinally on two faces of each tie rod. These 
strain gauges were connected to form two opposite arms of 
wheat-stone bridge circuit. The other two arms of the bridge 
were provided by mounting two dummy strain gauges on similar 


Steel plates to provide temperature compensation. 


The two jacks used to apply the vertical loads were 
connected to a control panel with pressure gauges to measure 
the pressure in each jack. The dial readings could be estimated 
to the nearest third of a division which corresponded to about 
650 lbs. in each jack. Loads were chosen to fall ona 


division mark. 


The strain gauges were all connected to a Dymec data 
acquisition system which read the voltage difference between the 
output terminals and printed this out on an automatic printer. 
The jacks and pull rods were calibrated before the test to find 
the correlation between the oil pressure, voltage difference 


reading and the applied loads. The calibrations were found to be 
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linear. 


7.8.3 Concrete Strain Measurements in the Walls 


Strains were measured on the face of the wall using 
Whittemore and Demec gauges. The first two stories of the wall 
had Whittemore points at 10 inch intervals and the upper two 
stories had Demec points at 8 inch intervals. The gauge lines 
were used, located at half inch and 4 inches from the outer face 
of the wall and at 1% inches from the inner face of the wall, 
respectively. A diagram showing these points is given in 
FIGURE 7.13. The points were placed on opposite faces of the wall 
in alternate stories to facilitate taking the readings from the 
working platforms. Two groups of three measurement points on 
two adjacent cross-sections made up a curvature measuring station. 
In FIGURE 7.13 the curvature measuring stations are labelled by 
numbers in circles placed in the center of the curvature measuring 


station. 


The Whittemore and Demec gauge readings could be repro- 
duced to the nearest 10 and 12.5 micro-inches/inch of strains, 
respectively. This will result in a maximum error of 21.4 and 
26.6 micro-inches/inch in the estimation of curvatures times depth, 

dt, for the wall sections, from Whittemore and Demec gauge 
readings, respectively. The concrete strain readings, in the first 
story, were found to be very erratic presumably due to the cramped 


working space between the gravity load simulator and the frame in 


this story. 
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7.8.4 Curvature Measurements in Columns and Beams 


Curvatures were measured at two locations at each end 
of each column and one location at each end of each beam using 


(C2) The details of a 


curvature meters developed by Chang 
typical curvature meter are shown in FIGURE 7.14. The meter 
consists of two steel reference frames which are attached to 
the members by four cap-screws screwed into the concrete surface 


and have 0.0001 inch dial gauges mounted at mid-depth of the 


Specimen to read the relative deflections of the ends of the arms. 


The gauge length was 8-3/8 + 3/8 inches and these gauges gave 
an average curvature over this length. In order to mount the 
meters at the joints, the curvature meter was modified for the 
attachment to the member. Similarly, the curvature meter used 
for top beam was modified to ensure against interference of the 
loading beam. The maximum error in the estimation of curvature 
times depth,dt, of the column and beam sections would be in 
order of 18.7 and 25 micro-inches/inch, respectively, if the 
dial gauge readings on the curvature meters are reproducable 


to nearest 0.0005 inch. 


The locations of the measuring stations are shown in 
FIGURE 7.13. Two adjacent reference frames form a curvature 
measuring station. The numbers in circles, shown in FIGURE 7.13, 


indicate the measuring station numbers. 


7.8.5 Rotation Measurements 


The rotations of points on the columns and beams near 


each joint and at the center of the wall joints were measured 
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electrically. The location of the rotation meters are shown in 
FIGURE 7.15. The measuring stations were labelled and are given 


by the numbers in circle, shown in FIGURE 7.15. 


The basic tool used to measure rotations is shown in 
FIGURE 7.16. This rotation meter is essentially the same as 


SU paeteayiy 


the rotation meter developed at Lehigh University 
the rotation meters consisted of a 0.03 inch thick spring steel 
strip with a 9.56 lbs .weight clamped to the lower end. The 

upper end of the strip was connected to two square bars which were 
firmly screwed to a base plate. The base plate was fastened to 
the specimen with 1/4 inch bolts screwed into inserts cast into 
the members during construction of the specimen. The weight 

tends to remain vertical as the structure rotates and bending 
strains are induced in the spring steel strip which acts as a beam 
under tension. By measuring the bending strain induced at the 


end of the strip the rotation can be determined using the formulas 


for a beam under tension. 


Four active strain gauges were used to form the four 
arms of a wheat-stone bridge to measure the strains as shown in 
FIGURE 7.17. Gauges 1 and 4 were on one side of the strip and 2 and 
3 on the other. This circuit compensated any temperature change and 
the strain measured is four times the average strain induced in 
the strip at the gauge point??? All the rotation meters were 
calibrated using a mechanical rotation meter, consisting of a level 


bubble and micrometer. The calibration was found to be linear and 


was approximately 525 micro inches/inch per degree of rotation. 
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7.8.6 Deflection Measurements 


The horizontal deflections were measured using scales, 
divided in 50ths of an inch, mounted on the beams. The scales 
were read with a transit. The scales mounted on the top three 
beams were 12 inches long whereas the scale on first beam was 
6 inches long. A scale was also mounted on the laboratory floor 


to check for movement of the transit. 


To measure the vertical deflections, magnetic bars 
were attached to the plates below the rollers on which the loading 
beam sat. Vertical wires attached to these bars extended to dial 
gauges mounted on magnetic stands on the floor of the laboratory. 
The deflection measured in this way included a horizontal compon- 
ent due to the lateral deflection of the specimen. This was 
corrected for in computing the deflections although the maximum 


error did not exceed 0.05 inch. 


A scale was mounted on one end of the fourth story wall 


and the lateral deflection was read with a second transit. 
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FABLE 733 


MEASURED DIMENSIONS OF THE SHEAR WALL-FRAME SPECIMEN 


; Steel 
Bier Cross-section Layer No. 
g Width | Depth From 
Member ‘ > : 
Inches ension 
Inch Inch Face 
60.5 6.0 4.5 ] 
2 


Beam 1 84.0 
(Col. end) 


(Wall end) 


Beam 2 84.0 
(Col. end) 


(Wall end) 


Beam 3 84.0 
(Col. end) 


(Wall end) ] 


Dist. of 
Steel 
Areas of | Centroid 
Steel From 
2 Tension 
in. Face 
Inch 
Oo222 0.750 
O.227ee | Sia oD 


4.5 ] 0.22270; | 0.750 
2 : Beek! 


0.875 
33625 
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TABLE 7.1 (contd. ) 


Beam 4 84.0 6.0 650 ] 
(Col. end) : 
(Wall end) ] 0. 
Bs 
Wall 1 60.5 24.0 ] OR 
2 OF 
3 0. 
4 0. 
5 0. 
6 ok 
7 0. 
8 0. 
Wall 2 6.0 24.0 ] 0, 
2 0. 
3 0. 
4 04 
Wall 3 60.0 6.0 24.0 0. 
and 4 0. 


be 
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er 
$32. 8f° 
2X8. oT 

avers 
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& 
- 
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TABLE 7 2 


AGE AND CONCRETE STRENGTH OF THE MEMBERS 


Concrete 
Strength 
in Comp. 


Age on 
The Day 
of Test 
Days 


Date of 
Casting 


Member 


Psi 


TABLE 7.3 


PROPERTIES OF REINFORCEMENT 


Ultimate 
Strength* 
Psi 


Strength 
_Psi 


Bar Size 


55,700 78,900 
53,000 80,800 
50,700 76,500 


54 ,900* 


* Taken from Mill Test 


Modulus 


82,800 


Split 
Tensile 


Strength 


of 
Elasticity 
Psi 


27;7x10 
27.4x10 
26.3x10 

27.9X10 
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SECTION IN FIG. 72 SCALE IN FEET 


FIGURE 7.1 SECTIONAL ELEVATION OF TEST SPECIMEN INCLUDING 


TESTING FRAME AND LOADING DEVICE 
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SECTION IN FIG. 7.) 
LOADING BEAM 
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LAB. FLOOR 


SECTION IN FIG. 7.1 1 aE ab DA BERS 
SCALE IN FEET 


FIGURE 7.2 SECTIONAL SIDE ELEVATION SHOWING SPECIMEN, TESTING FRAME, 


OBSERVATION PLATFORM AND VERTICAL LOADING DEVICE 
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FIGURE 7.3 SPECIMEN DURING TEST 
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FIGURE 7.4 TEST SPECIMEN SHOWING MAIN REINFORCEMENT 
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FIGURE 7.5 CROSS-SECTIONS OF FRAME MEMBERS 
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FIGURE 7.7 CLEVICE TO TRANSFER HORIZONTAL JACK LOAD 
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FIGURE 7.8 SPECIMEN UNDER CONSTRUCTION 
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FIGURE 7.13 LOCATIONS OF CURVATURE MEASURING STATIONS 
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CHAPTER VIII 


TEST AND ANALYSIS OF TEST DATA 


8.1 Introduction 


The testing apparatus, procedure and instrumentation used 
in the test of the shear wall-frame specimen were described in 
CHAPTER VII. The dimensions, structural details and material properties 
of the speicmen are also described in that CHAPTER. The computations 
carried out in reducing the test data are presented here. The test 


results are compared with the theoretical analysis in CHAPTER IX. 


Sa2 Description of Load Application and Measurements 


The shear wall-frame specimen was loaded on two subsequent 
days. After three load increments the Dymec electronic data recorder 
readings were found to be unreliable and the structure was unloaded 
and the test was postponed for a day while the Dymec was repaired. 


All the beams were cracked near the ends during this first test. 


The test was started again on the second day and was continued 
to failure. TABLE 8.1 presents a brief history of the second test 
including time and application of loads. All the deflections, curvatures 
and other readings have been referenced to the beginning of the second 


test. 


The testing and gauge reading was carried out by five groups 
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of two persons under the direction of one additional person. One 
group was fully responsible for applying and recording the loads, 

and printing the readings from the Dymec data logging apparatus. 

Each of the other four groups was responsible for reading and 
recording all dial gauges, Whittemore or Demec gauges and marking 
cracks on one story of the frame. The eleventh person was responsible 
for reading and recording the frame deflections and for the overall 


supervision of the test. 


Following the test the frame was photographed and the actual 


cross-sectional dimensions recorded in TABLE 7.1 were measured. 


8.3 Loads 


The horizontal and vertical loads were computed using 


hydraulic pressure-load calibration charts,established before the test. 


The electrical resistance strain gauge readings were found to be 
completely undependable due to inconsistencies in the Dymec data 
logging system. In computing the vertical loads, the weight of the 
loading beam, cross-beam and tie rods, totalling 1360 lbs., was 


added to the vertical jack loads. 


The measured loads from all six rams are presented in TABLE 
8.2. At load number 7 the load could not be held constant. For this 
reason the loads were measured before taking the other measurements 


(load 7) and again, after all the other measurements had been finished 


(load 7A). 
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8.4 Cross-Section Response of the Members in the Test Frame 


Theoretical load-moment-curvature relationships were computed 
for all the members of the test frame using the program described in 
CHAPTER IV, for the axial loads applied during the test. The load- 
moment-curvature characteristics for all sections were idealised to 
elastic-inelastic diagrams to comply with the assumptions made in 
CHAPTER III and the stiffnesses and ultimate moment capacities used in 
the theoretical analysis of the test frame were computed from these 
diagrams. On the other hand, the theoretical load-moment-curvature 
diagrams were used in computing moments from the curvatures measured 
during the test. The theoretical load-moment-curvature diagrams for 
all the columns, beams and walls are presented in FIGURES 8.1, 8.2 
and 8.3 respectively along with the idealised P-M-¢ diagrams used in 


the frame analysis. 


8.5 Curvature 
-8.5.1 Curvatures in Columns and Beams 


Curvatures were measured in the columns and beams using 
the device described in section 7.8.4. The average curvature and 
fibre strains can be computed from the formulae presented in 
APPENDIX D, using the dial gauge readings and parameters related 
to the measured dimensions of the curvature meters and the cross- 
sections. The initial dial gauge readings at load 1 were adjusted 
assuming that a linear variation of dial readings occurred under 
loads 1, 2 and 3 when only vertical loads were on the specimen. 
This was done to correct for slack, friction or initial movement 


of the dial gauges. TABLE D.1 in APPENDIX D lists the measured 
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curvature times the depth, ot, at all measurement stations, 
for all loads. Positive curvature indicates compression in 


inner face of the columns or the top face of the beams. 
8.5.2 Curvatures in Walls 


As mentioned in section 7.8.3. strains were measured, 
in the wall, using Whittemore and Demec gauges. Three continuous 
lines of these gauges were placed along the entire height of the 
wall, giving the average strain over 10 or 8 inch gauge length 
at three positions across the depth of section. Following the 
assumption that the strain distribution is linear over the depth 
of the wall, a least square line was fitted to the three readings 
to describe the strain configuration. The derivation of this 
line is described in APPENDIX D. At high load some of the Demec 
points were lost on the tension face of the wall, so that readings 
could be taken only at two positions. In these cases the strain 
distributions were based on these two readings. The curvature 
times depth, ot, from the fibre strains and geometry are listed 
in TABLE D.1 in APPENDIX D for all the sections in the wall and 


all the loads. Positive curvature indicates compression on the 


outer face of the wall. 


8.6 Moments 


Moment-curvature relationships were presented in section 8.4 
for all the frame members for all the applied axial loads. These 


curves were used to find the average moment in the gauge length, at 


each curvature measuring station. A linear interpolation was used to 
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estimate the moments falling between the computed points in these load- 
moment-curvature diagrams. The measured moments at the various 

stations are tabulated in TABLE D.2 in APPENDIX D, whereas the predicted 
moments are listed in TABLE 9.2 in CHAPTER IX. The error in the esti- 
mation of curvature times depth of section, ot, due to error in reading 
of dial gauges are given in CHAPTER VII. Using an error of 25 micro- 
inch/inch in the value of curvature times depth of section, ot, and 
FIGURES 8.1 and 8.2 the error in the estimation of moments are given 

by about 0.6 and 0.8 K-in. in the case of beams and columns, respectively 
in the intial part of the loading. Such error in the case of the wall, 
using FIGURE 8.3, is about 20 k-in. in bottom story to 12 K-in. in top 
story. These errors are reduced (see FIGURES 8.1, 8.2, 8.3) as the 
curvature of the section increases and the errors become very smal] 


when the hinges form. 


The moment diagram for the entire structure is presented in 
CHAPTER IX for each load. The moment diagram for each member was 
fitted to the individual moment values by eye. The measured and pre- 


dicted moments are compared in CHAPTER IX. 


8.7 Computer Analysis of Data 


The reduction of the test data, described in section 8.5 and 


8.6, were performed by computer at the University of Alberta. A 


computer program was written, for this purpose in Fortran IV language 
for IBM 0S/360 system. The flow diagram, nomenclature and the listing 


of the program are given in APPENDIX D. The various steps involved 


in the analysis are as follows: 
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1. The 'MAIN' program reads and writes all the data, such 
as cross-sectional dimensions and properties, member 
dimensions and properties, location and parameters of 
measuring stations and devices and all the dial gauges 


Or strain readings taken during test. 


2. The initial dial and strain readings were corrected as 


described in section 8.5 


3. The 'MAIN' program then calls the subroutine 'FSAC' 
which computes the fibre strains and curvatures at all 
the stations. These strains and curvatures are then 


printed out. 


4. The 'MAIN' program then calls the subroutine 'BMKP', 
which in turn interpolates the moments at all the 
stations for the curvatures computed in step 3 using 
the appropriate load-moment-curvature curves presented 


in section 8.4. These moments are then printed out. 


Another subroutine was also added to this program to compute 
the axial loads and moments on the section from the known fibre strains 
and curvatures computed in subroutine 'FSAC'. The axial loads and 
moments given by this subroutine were found to be unsatisfactory, since 
the dial gauge readings were accurate to 0.0001 inches while the axial 
loads and moments were found to be sensitive to gauge readings of one- 
hundredth of this amount. The moment computations based on curvatures 


are considerably more accurate as discussed in section 8.6, since the 
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curvature is the difference of fibre strains which tended to reduce 


the error. 


8.8 Shears 


The shear in the various members is equal to the algebraic 
sum of the moments at the ends of the member divided by th®-span of 
the member concerned. The shear for the various loads were computed 


from the measured moments. The P-A shears were calculated from 


measured deflections at all the floors and the measured vertical loads. 


The shear diagrams are plotted and compared to those from the analysis 


in CHAPTER IX. 


8.9 Deflections 


The deflection measurements at each floor level were 
reduced to get the actual deflections from the position of the frame 
at the start of the second test. The measured horizontal deflections 
at each floor level are tabulated in TABLE 8.3. The load deflection 
diagram for each floor is presented in CHAPTER IX and compared with 


the predicted deflections. 


8.10 Response of the Base of Wall 


As described in section 7.5 the base of the wall was anchored 


by two high strength bolts which reacted against pipe sleeves. One of 
these sleeves yielded during the test. A piece of the pipe of same 
length used in the wall connection was tested in compression. The 
measured load-shortening curve for the pipe is plotted in FIGURE 8.4. 
This curve was used to compute the wall base response to be used in 


theoretical analysis of the test frame. The connection is shown in 
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FIGURE 8.5(a). The total length of bolt which transferred the force 
from the wall base to the pipe, including the length of pipe and plate 
thickness was 32.5 inches before loading. The axial load applied to 
the wall was measured to be 50.2 kips. The wall width was 24 inches. 
If it is assumed that the center of the compression force, C, below 
the base plate coincides with one edge of the wall as shown in 

FIGURE 8.5(b), the relationship between the force in the pipe, T, and 
moment, M, in the base will be given by equations (8.1) and (8.2). 


Thele = M- 50.2x4 me. (oe) 


24T + 602.4 eae) 


or M 


The upward deflection, Ans of corner B of the wall is equal 
to the sum of the lengthening of the anchor bolt and the shortening 
of the pipe sleeve. The bolt diameter was 1.5 inches and the modulus 


6 psi. For each Value of T a 


of elasticity, E, was taken as 29.6 x 10 
value of M_ from equation (8.2) and a value of A, was computed. This 
was converted, in turn, into an equivalent curvature for use in the 
analysis by dividing by one tenth of the story height, since the 


analysis considered the wall divided into ten segments. 


The resulting moment-curvature-relationship for the base of 
the wall is plotted in FIGURE 8.6. The vertical initial tangent in 
the diagram corresponds to the period when the pipe sleeve was not 


subjected to any compressive force. 


For use in the shear wall frame analysis the moment-curvature 


relationship for the wall base was idealised as two straight lines as 
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Shown in FIGURE 8.6. The first line was chosen to correspond to 
first line of the moment curvature diagram presented in FIGURE 8.3 for 
the first story wall cross-section. The second line is an approxi- 
mation to the slope of the diagram plotted in FIGURE 8.6. This way 
the stiffness of the base was assumed to remain essentially the 
Same as that exhibited by the wall in the elastic range. The 
moment-curvature diagram for the pipe began to drop for curvatures 
in excess of 0.080 and this has arbitrarily been taken as the limit 


of the wall rotation in the analysis. 


8.11 Summary 


The procedure used to reduce the data from the test of a 
four story shear wall-frame structure have been outlined in this 
CHAPTER. Computer programs used in the data reduction are presented 
in APPENDIX D. The reduced data are summarized in APPENDIX D and 


will be discussed and compared to the analysis in CHAPTER IX of this 


report. 
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TABLE 8.1 
HISTORY OF TEST 


Load in Kips ae 
Load | 7... Total Total Defln. Remarks 
No. Vertical Horizontal | Inches 


Start of Test 


Fairly extensive beam 
cracking 


Base plate of wall lifted 
by 0.16 inch on west side. 
Crack at wall base ex- 
tended to 2nd Whittemore 
line 


Start of Readings 


End of Readings. Wall base 
plate lifted by 0.50 inch 
on west side. Crushing 
noted at wall end of 2nd, 
3rd and top beams. Load 
difficult to maintain. 


Wall base plate lifted by 
1-1/8 inches at west end. 
Top beam column joint 
failed. 


Horizontal loads released 
All loads released. 
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TABLE 8.2 


APPLIED LOADS DURING TEST 


Horizontal Lods in Kips Vertical Loads in Kips 
Load Ps 
0 0 0 0 0 0 
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CHAPTER IX 


COMPARISON OF PREDICTED 
VS. 


ACTUAL BEHAVIOUR OF FOUR STORY TEST FRAME 


9.1 Introduction 


A four story one bay reinforced concrete frame was tested 
under combined vertical and lateral loads. The test frame was des- 
cribed in CHAPTER VII and the results of the test were presented in 
CHAPTER VIII and APPENDIX D. In this CHAPTER a comparison is made 
between the observed behaviour and that predicted by the shear wall- 
frame analysis presented in CHAPTER III in order to illustrate the 
validity of the method of analysis. The comparisons have been made 


with respect to moments, shears and the sequence of hinge formation. 


The test results presented in this CHAPTER are those for the 
second application of load and have been referenced to the beginning 
of this test. The first loading was stopped because of failure of the 


electrical equipment. It should be noted that the beams were cracked 


in the first loading sequence. 


9.2 Behaviour of Test Frame 


The lateral load deflection diagram for the fourth floor of the 
test frame is presented in FIGURE 9.1. The data for both load applications 


are included in this diagram. The deflections have been referenced to the 
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beginning of the second load application. A positive deflection is 
one towards the applied loads (to the east). The loads referred 
to with letters and numbers in this figure are those for the first 


and second loadings, respectively. 


In each test, the first three load increments consisted of 
vertical loads of 0, 5.8 and 50.2 kips, respectively, applied to the 
center of the column and wall, at the top of the structure. Under 
these loadings the column shortened more than the wall and the struc- 
ture deflected to the west. Lateral loads were then applied to pull 
the frame to the east. Two increments of lateral load were applied 
in the first test. During the second lateral load increment, small 
cracks developed at the ends of the beams in the floors. Flexural 


cracks also developed in the wall in the bottom story. 


The residual lateral deflection at the fourth floor was 
0.18 inches to the east after all the loads had been removed. Of this, 
at least 0.10 inches was creep deflection which occurred during the 
time measurements were taken at loads D and E. The balance of the 
residual deflection was probably caused by the change in frame stiff- 


ness due to cracking of the beams. 


In the second test of the frame, the behaviour essentially 
paralleled that observed in the first test. More extensive cracking 
was observed in the beams at load 5 and some flexural cracks were 
observed in the bottom story of the wall. A diagonal crack was 
observed in the third and fourth beam to column joints. The cracks in 


the third floor joints are shown in FIGURE 5.1. Between loads 5 and 6 
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the pipe sleeve (see FIGURES 7.11 and 8.5) below the floor at the 
west end of the wall yielded and the west end of the wall gradually 
lifted off the floor, rotating about a point near the east end of 

the wall. For the rest of the test, the base of the wall probably 
behaved according to a moment-rotation relationship similar to the 
One shown in FIGURE 8.6. At load 6, diagonal cracks were observed in 


the second beam to column joint. 


At Toad 7, a diagonal crack was observed in the first beam 
to column joint and some bond distress appeared to be developing 
at the top beam to column joint. The flexural cracks in the bottom 
two stories of the wall extended to form inclined cracks which 
tended to cross pre-existing flexural cracks. The loads were very 
difficult to maintain and dropped off about 10 percent at this load 
increment. The measured rotations at the base of the wall indicated 
that the moment-rotation relationship for the base of the wall had 
reached its peak and had started to descend. Crushing was noted at the 


wall end of the second, third and top beams. 


As the deflections were increased after load 7, the loads 
dropped continuously. Horizontal loading was stopped when the maximum 
horizontal jack stroke had been reached. Between loads 7A and 8, 
crushing occurred at the bottom of the column and at all other hinges. 


The top beam-column joint failed in bond. 


FIGURE 9.2 shows the cracks in the frame after the test. 
FIGURE 9.3 traces the progress of hinging observed in the test. Plastic 
hinges were assumed to exist if very wide cracks developed and if the 


measured moments reached the computed plastic moments at any point ina 
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9.3 Analysis of Test Frame 


In the approximate analysis program the test frame can be 
represented by defining the stiffnesses of the left hand beams in 
the analytical model (FIGURE 3.1) as zero. The wall base response 
computed in section 8.10 was assumed to represent the behaviour 
of the first segment of the wall above the base of the wall. The 
idealized stiffnesses and ultimate moment capacities of various 
members used in the analysis were computed according to sections 8.4 
and 8.10 and are presented in TABLE 9.1. The ultimate moment capacities 
of the beams presented in TABLE 9.1 differ from those in TABLE 9.2 
Since the moment capacities of the beams at the wall end were used 


in the analysis. 


The shear wall-frame analysis described in CHAPTER III assumes 

that the rotation of the wall section occurs about its center line 

when computing the vertical deflection of points on the face of the 
wall. Since the base of the wall in the test frame rotated about a 
point close to the east edge of the wall section, the deflections of 
points on the west face of the wall were approximately equal to those 
assumed to occur in a wall of twice the actual width. For this reason 
the depth of wall section was assumed to be 48 inches for the purpose 


of this analysis. 


The computer program, described in CHAPTER III, assumes constant 
stiffness and ultimate moment capacities throughout the length of a 


member. Accordingly, it was necessary to modify the program to allow a 
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different moment-curvature relationship for the bottom-most segment 


of the wall. 


The resulting moments from this analysis and observed moment 
during the test are presented in TABLE 9.2 and are compared in section 


9.4. 


9.4 Comparison of Measured and Computed Moments 


The moment diagram for all the members are shown in FIGURES 9.4, 

9.5, 9.6 and 9.7 for loads number 4, 5, 6 and 7 respectively. The 
points and solid lines refer to the moments computed from the curva- 
tures measured in the test and the dashed line shows the predicted 
moments. In general the analysis underestimated the column moments 

in the top story and overestimated them in the bottom three stories 

in the initial stages of loading. For load number 7 the column moments 
are in fairly good agreement. As for the wall, the agreement in the 
measured and computed moments are fair except for the bottom story. 

The measured moments are shown for this wall but the points have not 
been joined because of the scatter in the points. This difference 

is probably due to errors in the readings resulting from the cramped 


working space adjacent to the bottom story. 


The beam moments at the wall end are overestimated, by the 
analysis, during the prehinging stages of loading and are in good 
agreement after hinges formed in the beams. The explanation of this 
may lie in the fact that the effect of local deformations at the beam 


to wall joints was ignored in the analysis. A method of correcting for 
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this in which the beam is assumed to extend into the shear wall, 

3 Of the beam depth beyond the face of the wall, has been suggested 
in reference (M10). The effective change in length of the beams 
may be increased by cracking of the members and the joint. To 

take these factors into account the test frame was re-analyzed 

by the analysis developed in CHAPTER III with the beam length 
extended by half its cross-sectional depth. However, the reduction 
in beam moments with respect to previous analysis was found to be 
less than two percent.at the ends of the two beams and 5 percent 

at the measuring stations which were located about 73" from the 


face of the wall or 104" from the 'end' of the extended beam. 


Another reason for the discrepancy in the moment diagrams 
lies on the fact that stiffnesses and moment capacities of the 
beams used in the analysis were derived for the sections near the 
wall end. In fact the stiffnesses and moment capacity of the beam 
varied along the length because the reinforcement was displaced 
during casting. The observed moments in the beams at the column 
to beam joints were found to be higher than the predicted moments 
Since the stiffnesses and moment capacities of the beams were 


higher at these joints than the wall to beam joints. 


It can also be seen from FIGURES 9.4, 9.5, 9.6 and 9.7 that 
the moment balance is better at the top two joints than at the other 
joints. This difference may be caused by difference in axial loads 
on the columns and walls. The beam shear produces a tensile force in 
the columns and a compressive force-on the wall, as a result the 
columns and walls carries different axial’ loads than the one measured 


by the hydraulic jacks. These discrepancies in the axial 
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loads are greater in the bottom story than at relatively higher 


stories. 


Part of the difference in the moments may be caused by 
ignoring the relative axial shortening of column and wall in the 
analysis. The measured axial shortening of the column and wall 
reveal that it can cause a moment ranging from 1.50 inch-kips at the 
ends of bottom beam to 1.80 inch-kips at the ends of the top 
beam. If this were included in the analysis however, the errors 


in the predicted moments would increase slightly. 


The other factor affecting the moments might be due to the 
way in which hinges are assumed to form in the analysis. The 
analysis assumes that point hinges form at the intersection of 
the members. at the joints but the test indicated that the hinges 
in the beams formed away from the joints by about half the beam 
depth and are spread over a finite length. This can be seen in 


PLGURES 5.1. and 9.2% 


The idealization of moment-curvature diagram can also 
produce some differences in the observed and predicted moment. As 
mentioned in CHAPTER IV, the moment-curvature diagram for reinforced 
concrete cross-section is not exactly bi-linear. The analysis uses 
a bi-linear moment-curvature diagram which is only an approximation 
of the exact moment-curvature relationship. In addition the exact 
moment-curvature diagram for all sections in a beam may be difficult 
to derive, since the various factors which affect this relationship, 


discussed in CHAPTER IV, may not be consistent at all cross-sections 
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along the length of a member. For example the strength of concrete, 
which depends on the compaction, void ratio, proportion of mix etc. 
cannot be maintained exactly constant along the full length of the 
member. These factors may give rise to different properties at 
different sections along the same member. Therefore, the moments 
derived from the measured curvature, using moment-curvature relation- 
ship for different members, presented in CHAPTER VIII, will give only 


an estimate of moments, not the exact value of the developed moments. 


9.5 Comparison of Shears 

The shear diagrams for all the columns and walls are compared 
in FIGURES 9.8, 9.9, 9.10 and 9.11 for load numbers 4, 5, 6 and 7, 
respectively. The dotted lines show the predicted shear whereas solid 
lines give the observed shear. The applied shear is calculated from 
the loads measured by horizontal jacks. The P-A shear has been computed 
from the measured horizontal deflections of the floor levels and the 
applied vertical loads measured by the vertical jacks. Total shears 
are the sum of the applied shear and the P-A shear. The frame and 
wall shears have been computed from the observed moments and are 
given by the sum of the top and bottom moments observed in a story 
divided by the story height. Since the observed moments in the bottom 
story wall are not known, as mentioned in section 9.4, it was not 
possible to compute the observed wall shear in this story. The P-A 


Shears and frame shears are plotted on an enlarged scale. 


It can be seen in FIGURES 9.8, 9.9, 9.10 and 9.11 that the 


sum of the observed frame and wall shears do not agree with the total 


shears. This is because of the fact that the observed moments in section 
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9.4 are not in good agreement which resulted in the discrepancy of 
the shears. If the P-A shears are computed as the sum of frame and 
wall shears minus the applied shear the result would be different 
from the P-A shears computed from the measured vertical loads and 
measured deflections. Since the measured deflections are in good 
agreement with the predicted deflections by the analysis, at least 
for loads 4, 5 and 6, the P-A shears plotted in FIGURES 9.8, 9.9 and 
9.10 are in good agreement and hence the total shear in the structure. 
For load number 7 the P-A shears are not in agreement, as can be 
Seen in FIGURE 9.11, because the difference in the measured and 
Predicted deflections are large. It should be remembered, however, 
that load 7 is close to the peak point of load-deflection curve 

(see FIGURE 9.1) and little difference in the load can cause a 
greater change in the deflections. Since the shears have been 
computed from the moments, the discrepancies in the moment diagrams 


are liable to be carried and echoed in shear diagrams. 


It can also be seen from FIGURES 9.8, 9.9, 9.10 and 9.11 
that the shear wall attracts a greater percentage of the shear in 
the bottom story whereas the frame seems to attract the greater 
portion of the shear in the top stories. This is partly because the 
column is stiffer in the fourth story than in the other three stories. 
The major reason, however, is the interaction between the shear wall 
and the frame. In tall structures, the wall will frequently have 
negative shears at the top of the building and as a result the frame 


may have very high shears in top stories (K2, M8). 
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9.6 Comparison of Deflections 


The predicted and observed deflection diagrams for the various 
floor are given in FIGURES 9.12, 9.13, 9.14 and 9.15. The theoretical 
analysis seems to overestimate the deflection in the beginning and 
under-estimates it in the later stage of loading. The explanation of 
this lies on the idealization of the "elastic" portion of the moment- 
curvature relationship for the section. As shown in FIGURE 8.2 for 
example, the idealization of this curve leads to underestimate the 
stiffness at low loads so that the predicted deflections initially 
are in excess of the measured values and tends to overestimate the 
stiffness at later stages of loading resulting in an underestimation 
of the subsequent deflections. The observed deflections include the 
shear deformations whereas the analysis consider only the bending 
deformation. This may cause little discrepancies in the measured and 
predicted deformations. However, the agreement in general is satis- 
factory. The falling branch of the load-deflection curve is not possible 
to find by the analysis since it does not consider the falling branch 


of moment-curvature relationships for the cross-sections. 


9.7 Formation of Hinges 


FIGURE 9.3 compares the plastic hinges observed in the test 
to those predicted by the analysis. No hinges were found in the 
theoretical analysis, or in the test up to load 5. The observed and 
predicted hinges for load 6 and 7 are compared in FIGURE 9.3. No hinges 
were observed at load 6 in the test although the analysis detected 


hinges at the wall end of alt the beams. The observed moments were very 
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close to the hinging capacities at this load, however. At load 7 
all the hinges are comparable. In general the formation of plastic 
hinges detected in the analysis were close to that observed in 

the test. It is also observed that the hinges form away from the 
joints by about half the member's depth and spread over a finite 


length (see FIGURE 9.2). The hinges form with severe cracking of 


concrete. 
9.8 Summary 


The predicted and observed behaviour of the test frame 
was compared in this CHAPTER. The comparisons of moments, shears, 
deflections and formation of plastic hinges were ria in sections 
9.4, 9.5, 9.6 and'9.7., respectively. In general the analysis seems 
to give a good description of the behaviour. It was found that the 
various assumptions such as the idealization of moment-curvature 
relationship, idealization of material properties, consideration of 
point hinges, neglecting shear deformations etc. have very little 
effect on the analysis. However, the method of testing does not 
focus on the assumptions of lumping and loading made in CHAPTER III. 
These assumptions have been discussed in CHAPTER VI and a lumping 


procedure has been derived. 
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STIFFNESS AND MOMENT CAPACITY OF VARIOUS 
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FIGURE 9.2 SPECIMEN AFTER TEST SHOWING CRACKS 
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FIGURE 9.3 HINGE PATTERN IN TEST FRAME 
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CHAPTER X 


SUMMARY AND CONCLUSIONS 


10.1 Summary 


A method of analysis which traces the second order elastic- 
inelastic behaviour of shear wall frame structures has been presented. 
The analysis considers a simplified lumped model derived in CHAPTER VI. 
The iteration procedure has been used for the analysis. Equilibrium 
was formulated on the deformed structure to consider the effect of the 
P-A moment. Constant concentrated vertical load on each story were 
considered. The horizontal loads were assumed to be concentrated at 
the floor levels and were applied in increments. The analysis allows 
for the finite width of wall but neglects the axial and shear deform- 
ations and flexibility of joints. A computer proaram in FORTRAN IV 
language for IBM system 0S/360 was developed to perform the computa- 
tions and is presented in CHAPTER III. The assumptions in the analysis 
are discussed in CHAPTERS IV, V and VI. A design example has been 


presented in APPENDIX A. 


Elastic-perfectly plastic and elastic-strain hardening 


moment-curvature relationships were assumed for columns and beams, 


and for walls respectively. A method was developed to find the load- 


moment-curvature relationships for reinforced rectangular cross-sections. 


A computer program was written for the purpose and was presented in 


CHAPTER IV. 
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The lumping procedure for the structure has been derived 
in CHAPTER VI. The total uniformly distributed load on each story was 
replaced by a concentrated load. In order to take this into account, 
the plastic moment capacities of the beams were reduced by a reduction 
factor a. The analytical checks were also presented to see the 


validity of the approximate analysis. 


A full scale test on 20 feet high, 4 story, one bay 
Structure was described in CHAPTER VII. The test was carried under 
fixed vertical load and incremental lateral load. The machine pro- 
cessing of the test data was described in CHAPTER VIII. Reasonable 
correlation was obtained between the results of the test and analysis. 


The comparison was presented in CHAPTER IX. 


10.2 Conclusions 


The analytical checks in CHAPTER VI and the frame test 
Suggest that the analysis gives a good indication of the behaviour 
of the shear wall frame structure. The model used is relatively 
Simple to develop. The test and analysis have clearly shown the 
necessity of considering the interaction between the frame and the 
Shear wall. The shear developed in the top of the frame was found to 
be much larger than the shear in the wall. The test does not show the 
development of negative shear in top of the wall for the model tested, 
but it did show that tendency. The analysis has shown the development 
of negative shear on the top of the wall in high rise structures, in 
CHAPTER VI, which suggests that in many structures the frame has a 


tendency to limit the deflection of the shear wall. Therefore it is 


Boe 5 ' ny ; 5 
‘ mn ; 
Py 5 \ : a i ° ent in 
i a | q te 
, fs i Tie he i id mf 7 
7 ” ~ 
4 f wr A oe ' 
? ® ; 4 x Aes: : 


hind 
baviveb asad bs a ust DUNFe std ot arubspore onto " 


an ; 
26W HOS 2" IBS 0 bsot botudtyte hh: ‘ylmiot tau [stod oAT VI | AATSP 


40. 


one ee | 
| einyoose ofnt ett Shed ot “oby0. al: baot badsrdn99n02 6 va be 
j oe | R ty , com ave, 
| nottoubsy 6 vd baoubay ayew emsad Say tag gota 58969 ‘Snemam 2 tt26 
ri if te 7 
r, ie ont 443 ot DSTNsS2arvg ‘oats ovew. 24 sora Tot tyfens aaT x 
: Be 
me -ateyens. ‘9emtxorgas ont fo wat 


a 


: yy 
ysd sto .vyot2 f aes 08 Os nO fe st 5 feat: vt Ae 


nu botyyeo 26W hess ST My 9aTIAHD at badiroaeb: 26H sy 
fis 
1 ¥ i 
-org onftdosm sat .,.beol tev asf i lin Tero if bas oath ts a 
é bi. fe 
g{dsroessh \iTIN AF apy) nh bedtnozsb 25 tsb: dead, ow 
A 
afeylsas brs teat) odd to at [yeay ons nsswted pontaddo 26 ag ots 
; Ripa 
iki AaTIONY- nt begneze1q 284 moat ‘ 
: hws a7) » 
a 4 4 j 7, cy " 
Na y : y 
j mI Bei ae 4 a Poe - " oe 
¥ | lad f i, 1 ee OT ad a a 
| Feat omeyt ond ‘bie IV s3TAAHD nt sien tot stand” ‘sit 
yuotvsded sit to nofseaibint vest abate atavtond ed, wl is 
Pa ar oe > sa i ies ibe ek 
Pal ¥ viavitslet ef beeu fabom sat. smi sme ‘(sw "16s ie a 
ks) snd wore ulipato. vad steels es ‘teed ot avab of ‘sf 
urs nee Mas ; oo SDM Py 
rin 1) is bore stort sits negivtad ‘notio6 Bont ot ents abr ea 240 stat 


want, a Ry > re Pe) 
ei Ae “ot bows? 260 eine ‘wt to et. sty nF ano reveb 68 att wink 
| Pa foe 
i om wore ton 2a0b Seay eit Sen 16 am yep1s! doum 9 
Bi qe b 


| oases fobom orld 10} i hsw oi oO rf 
) a 
ffeasaren aid nworte ann etagtans 


+ 
at 1 


ve exw nde A i ; 5 ii ead ‘ nt no. * soa evt por 
ve u 2 92h pte at Ife mw sit TO got 9 ff. x6 ove $6 8 


a 
ai oe aa 


ie 2 En S@iBnt “an Ht mont $6 k ' So 2e0pU2 ih hele. , 7 ee 
he 4. if : iw - % 2 } _ " ‘ ’ ; 


ths ameter Te 


2pe 


necessary to evaluate the interaction between the frame and the shear 
wall if true behaviour of the structure has to be established. Neg- 
lecting this factor could lead to excessive deformations and possible 


failure of the affected frame members. 


The test indicated some differences in the observed and 
predicted moment and shear diagram. The beam moments were over- 
estimated by the analysis during the prehinging stage of the loading 
and are in better agreement after the hinges have formed. The hinges 
in the test frame were observed to form away from the joint and were 
spread over a finite length. The analysis tended to overestimate the 
deflection of the test frame in the beginning and underestimated in the 
later stage of loading. However, the test has clearly indicated that 
the assumption of point hinges and rigid joints and ignoring the shear 
deformations has little effect on the analysis. The analysis slightly 


overestimated the overall stiffness of the structure. 


The analysis does not consider the unloading branch of 
load-moment-curvature relationship of the cross-section. The load- 
moment-curvature diagram shows the unloading of the cross-section after 
reaching its ultimate capacity. The development of load-moment-curvature 
diagram indicated that maximum compressive strain corresponding to 
ultimate moment capacities varies and depends on the axial load, per- 
centage and placing of reinforcement. The effect of the tensile stress 
in the concrete was found to be considerable especially in the case of 
under-reinforced concrete cross-sections. However, the effect was 
considerable when the tensile strain of the extreme fibre is close to 


the cracking strain. The stiffnesses and ultimate moment capacities 
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were found to be dependent on axial load on the cross-section. At 
high axial load, the section did not have sufficient rotation 
Capacity to permit plastic hinging. However, the load-moment- 
curvature relationship developed here was found to be satisfactory 


in estimating the stiffnesses and ultimate moment capacities. 


The load-deflection characteristics of the structure was 
found to resemble well with the lumped model. The lumped model 
gives accurate description of behaviour of the buildings when sub- 
jected to lateral load only. For the buildings subjected to combined 
load the lumped model gives the latter softening but earlier failure 
of the structure. In such cases the predicted ultimate loads were 
found to be underestimated but the degree of conservatism was not 
in considerable error. The predicted ultimate load by the approxi- 
mate analysis was in good agreement with the other analyses (C6,P2, 
G2). Thus, the procedure developed for lumping the building is found 


to be satisfactory. 


In the development of the lumping procedure for walls it 
has been found that the moments and shears in the individual shear 
wall cannot be determined simply by the ratios of their respective 
stiffnesses. The moment diagram in individual members of the building 


should be found from the deflected shape of the structure. 


10.3 Recommendations for Future Research 


The method of analysis can be extended to consider the 


axial and shear deformations. The analysis can also be extended 


to consider the unloading branch of load-moment-curvature relationships 
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and insufficient rotation capacity of the section. 


Much more extensive large scale laboratory testing of 
Shear wall frame structures is required to Study the behaviour more 
completely. Studies of multi-story and multi-bay structures will 


illustrate the effect of lumping procedure used in this analysis. 


Test on cross-sections are required to accurately 
determine the behaviour of the sections. These should include tests 
On cross-sections with small percentage of reinforcement to determine 
accurately the influence of the uncracked concrete between tension 
cracks on the stiffness of the member. Test on cross-sections 
composed of high strength concrete and reinforcing steel are also 
required, since the multi-story construction may require high 


strength of the materials. 
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APPENDIX A 
EXAMPLE OF A REINFORCED CONCRETE 
BUILDING DESIGN” 
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A.1 Introduction 


A preliminary design of 20 story building 160' x 62' in plan, 
shown in FIGURE 6.28, has been presented here. The objective was to 
design the building for vertical load and investigate the behaviour 
of the structure, under combined lateral and vertical load, using the 
approximate analysis program developed in CHAPTER III. The axial 
loads were computed from the normal assumptions of tributary areas. 

The checker board loadings were used for the consideration of live 
load in the design. The members were proportioned by ultimate strength 


design procedure (A!) , 


A.2 Design for Vertical Loads 


The building consisted of seven frames and two shear walls 
in the direction of wind as shown in FIGURE 6.28. Each frame had 
four column lines interconnected with beams at each floor level. There 
was no cross girders. The thickness of the slab was assumed to be 
4 inches. The assumed vertical loadings were 100 psf live load on 
the floors and 50 psf live load on the roof. For simplicity it was 
assumed that slab and beams are not monolithic in construction so that 


the beams could be designed as a simple rectangular beam rather than 


as a T-beam. 


The chosen strength of concrete and steel were 4 and 60 ksi 


respectively. 


The beam moments were computed using coefficients given in 


section 904 of ACI Building Code (1963). The gravity load on the beams 
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were taken as (1.5D1 + 1.8LL) (Al). Clear span length of 19 feet were 
used for the computation of moments in negative and positive moment 
region. However, to maintain the uniformity and simplicity all the beams 
were proportioned, for the maximum moment of w1'*/106, by ultimate 
Strength design procedure. @ is the capacity reduction factor defined 

in section 1504 of ACI Building Code (1963) and was taken as 0.9. 

While proportioning the beams,equal amount of reinforcement were 

assumed in the tension and compression face and the width of beam 

was taken approximately half the depth. The percentage of reinforce- 


ment in the beams were selected to eatisfye | 


p,f 
Giases ere = sD ol 8 a (Aa) 
c 
where P, = percent of reinforcement based on gross area 
ur = yield strength of reinforcement 
and te = concrete strength 


The section of the beams designed are presented in TABLE A.1. The 
moment of inertias and plastic moment capacities are also presented 


in TABLE A.1. The moment of inertias were computed using equation 


(4.4) suggested in CHAPTER IV. 


The analysis of column loads were based on the assumed tri- 
butary areas. The moments were computed for the full dead load and 
the live load placed as a checker board. The columns were designed 
for the axial loads and moments as per requirement of section 914 of 


ACI Building Code (1963). The columns were chosen as square tied 
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column and had approximately four percent reinforcement. The 
Capacity reduction factor,, defined in section 1504 of ACI 

Building Code (1963), was taken as 0.7. Also the strength reduction 
factor for the length of compression members, defined in section 

916 of ACI Building Code (1963) was used. All the columns were 
proportioned by the ultimate strength design procedure H)). 

The sections, thus obtained, are presented in TABLE A.1. The 
moment of inertias of the column sections were computed by equation 


(4.4) suggested in CHAPTER IV. TABLE A.2 gives the axial loads on 


each column. 


The walls at the two ends of the building were taken as 
6" thick and 62' wide. These were found to be quite adequate for 
the axial loads and moments on the wall giving much extra capacity 


than required. A nominal reinforcement of one percent were used. 


A.3 Design for Combined Loads 

The building designed for vertical loads in section A.2 
consisted of seven similar frames and two walls in the windward 
direction. This building was reanalyzed to test its performance 
under the combined vertical and lateral loads. Lateral loads were 
assumed to be the wind load of 25 psf and were assumed to be con- 
centrated at the floor levels. For the purpose of combined loads 


(Al) 
analysis the vertical loads were taken as 1.25 (D1,.+ -LL) f 


The building was lumped by the procedure developed in 
CHAPTER VI. The two walls were lumped together to give the wall 


system in lumped model, shown in FIGURE 3.1. Since the wall wil] 
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have a lot of openings for windows etc., an equivalent depth of only 
45' was considered for the purpose of computing the stiffnesses and 
moment capacities. The effect of reinforcement were neglected. 
However, for computing vertical loads, solid wall of 6" x 62' were 


considered. 


The frame system in the lumped model shown in FIGURE 3.1 was 
obtained by lumping all the seven frames. All the columns and beams 
were lumped by the procedure developed in CHAPTER VI. The values of 
the factor a, defined in CHAPTER VI for the case of uniformly dis- 
tributed load on beams, were taken as 0.66 for roof level beams and 
0.61 for all other beams. These factors were computed from FIGURE 


62103 


The frame system and the wall system were connected by the 
hinged link beams at all floor levels. Thus the building could be 
represented by the lumped model for the purpose of approximate 
analysis developed in CHAPTER III. The moment of inertias of the 
various members of the lumped models are given in TABLE A.3 and the 
corresponding moment capacities are given in TABLE A.4. ‘The vertical 
loads acting on each story and the horizontal loads acting at each 


floor level of the lumped model are also given in TABLE A.3. 


The lumped frame thus obtained was analysed by the program 


developed in CHAPTER III under fixed vertical load of 1.25 (D1 + LL) 


and incremental horizontal loads. The load-deflection curve, obtained 


for the top of the building, is presented in FIGURE A.1. A load factor 


of 1.35 was found for this building. Under the action of combined 
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loads the ACI Building Code (1963) requires a load factor of 1.25. 
The story deflection at working load factor 1.0 and at a load 
factor of 1.25 are plotted in FIGURE A.2. The hinges formed at the 
ultimate load are shown in FIGURE A.3. 


A.4 Discussion and Conclusions 


The vertical load design generally indicated that the 
section can be chosen in various ways to satisfy the stiffness 
and moment capacity requirement. The chosen relationship between the 
various parameters of the section seems to serve the purpose when 
the building was subjected to combined vertical and lateral load. 
Under (1.25 (D1 + LL) + 1.0WL the maximum story to story deflection was 
0.25 inches or 1/576 of the story height. This is near the upper 
limit of acceptable values. Under the service load condition of 
1.0 (D1 + LL + W1) the deflections would be smaller and it may be 
desirable to check these in an actual design. However, the building 
as designed has sufficient stiffness under the action of combined 
loads and any alteration in sizes are not necessary from the stiffness 


point of view. The economical aspects have not been investigated. 


The behaviour portrayed under the action of combined loads 
are typical of a shear wall structure with the upper floors of the 
frame restraining the shear wall from additional deflection and the 


lower floors of the frame being restrained from additional deflection 


by the shear wall. 


All hinges formed only in the beams thus giving a strong 


column and weak beam design. The first hinge in the structure was 
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indicated at a load factor of 1.05. The first wall hinge at the 
base was indicated at a load factor of 1.30. Due to the importance 
of the wall to the overall stability and strength of the building 
the load factor at first hinging of the wall should exceed the 


desired load factor for the structure. 


As a result of the above design it can be said that the 
structural system employed for this building does not require extra 
stiffness in the frame or shear walls, to aid in the resistance of 


lateral loads, over and above that required by the gravity loads. 
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TABLE A.1 


SUMMARY OF THE BEAMS AND COLUMNS SECTION 
FOR A TYPICAL FRAME 
Beam Sections 


Area of 
Size Reinf. Moment of Plastic Moment 
Inch Sq .& din. Inertia, in Capacity, K-in. 
Roof level 12x24 3.60 7970 2070 


All other floors 14x27 4.40 12790 27/0 


7 Column Section 
Interior Column Exterior Column 


Plastic Area Plastic 

Moment Moment] Size of Moment Moment 

of Inertia | Capacity In. Reinf. Jof Inertia | Capacity 
bait K-in. in’ in’ K-in. 


106800 27x27 29.24 40900 26400 

106800 27x27 29.24 40900 26000 

80000 24x24 22.50 25100 18800 

80000 24x24 22490 25100 18650 

62200 225x224 20.00 19550 15400 

62200 225x225 20.00 19550 15380 

285X285 50800 21x21] 17.90 15060 12600 
283X285 50800 21x21 1/7 £90 15060 12500 
26x26 35100 1925x1935 15324 112g 10350 

10 26x26 35100 19%sx 1 9s 15.24 11120 10250 
1] 24x24 25200 18x18 12.70 7810 7860 
12 24x24 25200 18x18 12270 7810 7820 
13 22x22 18020 16x16 10.24 5030 5570 
14 22x22 18020 16x16 10.24 5030 5500 
75 19x19 9990 14x14 8.00 3020 3700 
16 19x19 9990 14x14 8.00 3020 3620 
17 16x16 5030 12x12 5.76 1595 2280 
18 16x16 5030 A ow Ae) 1595 2180 
19 12x12 1595 12x12 ave 1895 2020 
20 i2x12 1595 12x12 270 1595; 123 1800 
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TABLE Ad 
VERTICAL LOAD ON COLUMNS OF A TYPICAL FRAME 


Interior Column Exterior Column 


Dead Load Live Load Dead Load Live Load 
in kips in kips in kips in kips 
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TABLE A.3 


LOADS AND MOMENT OF INERTIA OF VARIOUS MEMBERS 
OF LUMPED MODEL 


Moment of Moment of Moment of | Moment of 
Horizontal | Vertical |Inertia of | Inertia of | Inertia of | Inertia of 
Load in Load in Lumped Lumped Hinged Lumped 
Kips Kips i Column Link Beam Wall 
4 4 


in in in? in 
537180 2067800 
537180 2067800 
537180 1471400 
537180 1471400 
537180 1144500 
537180 1144500 
537180 922040 
537180 922040 
537180 647080 
537180 647080 
537180 462140 
537180 462140 
537180 322700 
537180 322700 
537180 182140 
537180 182140 
537180 92750 
537180 92750 
537180 44660 
334740 44660 
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TABLE A.4 


MOMENT CAPACITIES OF VARIOUS MEMBER 
OF LUMPED MODEL 


Moment Capacity Moment Capacity Moment Capacity 
of Lumped Beam of Lumped Column of Lumped Wall 
K-in. K-in. K-in. 


[e>) 


1260000 
1229000 
893200 
886200 
729400 
729120 
620200 
616000 
482300 
478100 
373240 
372680 
282380 
277200 
129780 
127120 
109200 
104720 
59920 
52500 
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DEFLECTION IN INCHES 
FIGURE A.2 STORY- DEFELECTION DIAGRAM 
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FIGURE A.3 HINGE PATTERN IN THE ANALYTICAL MODEL AT ULTIMATE LOAD 
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APPENDIX B 


DEVELOPMENT OF THE APPROXIMATE 
ANALYSIS 
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B.1 Derviation of Joint Rotation Equation 


The frame system has been forced into the deformed Shape of 
the wall shown in FIGURE 3.6. Applying the Slope-deflection equations 


5 


the moment Mews , at the end B. of the right beam is given by: 


Maui = Spwi 49; * 28ypgy - © R'y) 


WBi 


where Spy = EI/L is the stiffness of the right beam (B.W. ) and Ri. 
is the vertical chord rotation of the right beam connecting the 

column to the wall. Also in the above equation 6. is the rotation 
of the frame joint and OB is the rotation at the wall end of the 


right beam. 


Considering the left beam B.F. and applying the slope 
deflection equation, assuming that the joint rotation at the two 


ends B. and F. are equal, the moment, Moras at the end F. is given 


by: 


where Spr. = EI/L is the stiffness of the left hand beam, B.F.. 


The column BiB. 4, between the ith and the i+] th story 


has a stiffness of So (441) and the chord rotation of this column is 


R The moment, Mo(a41)? at the bottom end B. of the column 


Line 


B. Be 4] is given by: 
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Mo(itt) = Sccist) (40; * 2654, - 6 Riyy) 


where 0544 is the rotation of the frame joint at the (i+1) th level. 


The moment, Mage at the top end, B.. of the column 


B. Bey is given by: 


Mega iSa. (40, + 20. _ 


ci Ci , 6R, ) 


1 


where So7 = EI/L is the stiffness of the column B. B._, and R. 
is the chord rotation of the column. 85] is the rotation of the 


frame joint at the (i-1) th level. 


Equating the sum of the joint moments to zero 


po ere” Met Mote) 


Substituting the values of the end moments from the above 


equations the joint equilibrium equation becomes: 


; (40, + 20.13 — 6 R'.) + 6 Spr 6. + 


Sie iS cg) 


(46. Bees 6 R-41) = 0 
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after simplification: 
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Equation (B.1) can be written in the following simplified 


form: 
pigs, fe Bet Gerad 
j ae a Lae aa i ae Deco 
Where in the numerator, 
Me a eg etn) 
Baia ae hoe 20455) 
deans Ei Gea A Meas 
and D = 0 


And in the demoninator, 


C = 4S. (i+1) 
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Equations (B.1) and (B.2) are valid for the elastic analysis 
only. When a plastic hinge forms at any one of the potential hinge 
locations, the joint rotation equation must be modified. TABLE B.1 
lists the substitutions to be made in Equation (B.2) for the formation 
of hinges in the various members, so that the joint rotation equation 


will conform to the particular hinge pattern considered. 


In TABLE B.1, MP, MP awa are the plastic moment capacities 


Rik 
of the beams B.F. and B.W. » and MP is the plastic moment capacity 


(reduced for axial load) of the column. 
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TABLE B.1 


MODIFICATION OF ELASTIC SLOPE-DEFLECTION 
EQUATIONS FOR JOINT EQUILIBRIUM 


Substitution to be Made 
in Equation (B.2) 


Member Hinge Location 
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TABLE B.1 (contd. ) 


Substitution to be Made 
Hinge Location 
in Equation (B.2) 
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B.2 Nomenclature for FORTRAN IV Program 


BAKA A SUBROUTINE SUBRPOGRAM 
BMOM BOTTOM MOMENT IN A SEGMENT OF THE SHEAR WALL (KIP-IN) 
BMOMF BOTTOM MOMENT IN THE BOTTOM OF THE LOWEST SEGMENT 


IN A STORY (KIP-IN) 

CON CONVERGENCE LIMIT FOR THE DEFLECTION AND ROTATION 

CDEFF FLOOR LEVEL DEFLECTION OF SHEAR WALL AFTER APPLYING 
CONVERGENCE FORMULA 

CROTF FLOOR LEVEL ROTATION OF SHEAR WALL AFTER APPLYING 
CONVERGENCE FORMULA 


DEF DEFLECTION OF WALL AT A SECTION (IN) 
DEFF DEFLECTION OF THE WALL AT THE FLOOR LEVEL (IN) 
DMBF MOMENT AT THE COLUMN END OF THE LEFT HAND BEAM 


COMPUTED IN THE PREVIOUS LOADING CONDITION (KIP-IN) 

DMBW MOMENT AT THE COLUMN END OF THE RIGHT HAND BEAM 
COMPUTED IN THE PREVIOUS LOADING CONDITION (KIP-IN) 

DMCB MOMENT AT THE BOTTOM END OF COLUMN IN A PARTICULAR 
STORY COMPUTED IN THE PREVIOUS LOADING CONDITION 
(KIP-IN) 

DMCT MOMENT AT THE TOP END OF A COLUMN IN A PARTICULAR 
STORY COMPUTED IN THE PREVIOUS LOADING CONDITION 

| (KIP-IN) 

DMFB MOMENT AT THE LEFT END OF THE LEFT HAND BEAM COMPUTED 
IN THE PREVIOUS LOADING CONDITION (KIP-IN) 

DMWB MOMENT AT THE WALL END OF THE RIGHT HAND BEAM COMPUTED 
IN THE PREVIOUS LOADING CONDITION (KIP-IN) 


(WI-9I) JAW MATZ BHT IO THAMOB2 A wt -TWaNOM worTOg 
TwaMod2 T23wo. 3HT 30 MOTTOS 3HT Wi THaMON woTTos ” 

(i404) yar A MI s 

MOITATOA OMA WOITIBIIIG SHT OF TIMTY 3OMagA3vMOD” 7 
MIVIIGA AATAA JUAN SABE 2 MOITOSIAAG J3V3J A009 | 
AJUMAOA soaaRaVHOD 

DMIYI9IA AATAA JAW anahe 70 OLTATOR 43V34 900.3 
AguMA03 3OM3093 M09 

(MI) WOLTI2 A TA JAW 40 nO1 793.1730 ; 

(M1) J3V31 90019 3HT TA JAW SHT 90 WOrTI92330 bY 

MAga GMAH T3ay ant 40 ana WMUI03 3HT TA TH3MOM \ Li 
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Aa 


HLIR 


HS 


HSF 


IK, IX 


IVL 


KB 


KC 


MAX 


MI 


WIDTH OF SHEAR WALL (IN) 

YOUNGS MODULUS OF ELASTICITY OF THE FRAME (KIP/IN¢) 
YOUNGS MODULUS OF ELASTICITY OF THE WALL (KIP/IN‘) 
APPLIED FORCE ON THE WALL FROM FRAME ANALYSIS (KIP) 
HORIZONTAL DESIGN LOAD ACTING AT FLOOR LEVEL (KIP) 

A SUBROUTINE SUBPROGRAM 

FRAME FORCE AT THE FLOOR LEVEL (KIP) 

FINAL FORCE ON WALL AT EACH FLOOR LEVEL (KIP) 

LENGTH OF LEFT HAND BEAM (IN) 

LENGTH OF RIGHT HAND BEAM (IN) 

PERCENTAGE OF ORIGINAL HORIZONTAL LOAD TO BE INCREASED 
IN THE ELASTIC RANGE 

PERCENTAGE OF ORIGINAL HORIZONTAL LOAD TO BE INCREASED 
IN THE INELASTIC RANGE | 

STORY HEIGHT (IN) 

HEIGHT OF FLOOR LEVEL FROM BASE (IN) 

DUMMY CONSTANTS 

DUMMY CONSTANTS FOR STOPPING THE PROGRAM IF THE 
DEFORMATION EXCEEDS CERTAIN SPECIFIED LIMIT 

SPRING CONSTANT AT THE BASE OF THE SHEAR WALL (KIP-IN 
/RAD) 

SPRING CONSTANT AT THE BASE OF THE COLUMNS (KIP-IN/RAD) 
PRODUCT OF NUMBER OF STORY AND NUMBER OF SEGMENTS IN 


A STORY 
MAXIMUM NO. OF CYCLE TO BE PERFORMED FOR ANY ITERATION 


PROCESS 


MOMENT OF INERTIA OF SHEAR WALL (IN‘“) 
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MIC 
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MOMBW 
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MOMENT OF INERTIA OF LEFT HAND BEAM (IN*) 

MOMENT OF INERTIA OF RIGHT HAND BEAM (IN*) 

MOMENT OF INERTIA OF COLUMN (IN?) 

MOMENT OF THE SEGMENT OF THE SHEAR WALL (IN*) 

NO. OF PROBLEMS TO BE SOLVED 

MOMENT AT THE BOTTOM OF A COLUMN IN A STORY (KIP-IN) 
MOMENT AT THE TOP OF A COLUMN IN A STORY (KIP-IN) 
COUNTER FOR ITERATION WITH THE EFFECT OF AXIAL LOAD 
MOMENT AT THE BASE OF THE SHEAR WALL (KIP-IN) 

MOMENT AT THE COLUMN END OF THE LEFT HAND BEAM IN 

THE CYCLE UNDER CONSIDERATION (KIP-IN) 

MOMENT AT THE COLUMN END OF THE RIGHT HAND BEAM IN 
THE CYCLE UNDER CONSIDERATION (KIP-IN) 

MOMENT AT THE LEFT END OF THE LEFT HAND BEAM IN THE 
CYCLE UNDER CONSIDERATION (KIP-IN) 

MOMENT IN A STORY DUE TO AXIAL LOAD (KIP-IN) 

TOTAL MOMENT AT FLOOR LEVEL ON THE SHEAR WALL DUE TO 
END MOMENT AND SHEAR FROM THE RIGHT HAND BEAM (KIP-IN) 
MOMENT AT THE RIGHT END OF THE RIGHT HAND BEAM IN 

THE CYCLE UNDER CONSIDERATION (KIP-IN) 

PLASTIC MOMENT CAPACITY OF COLUMN IN A PARTICULAR 
STORY (KIP-IN) 

PLASTIC MOMENT CAPACITY OF THE LEFT HAND BEAM (KIP-IN) 
PLASTIC MOMENT CAPACITY OF THE WALL SEGMENT (KIP-IN) 
PLASTIC MOMENT CAPACITY OF THE SHEAR WALL (KIP-IN) 
PLASTIC MOMENT CAPACITY OF THE RIGHT HAND BEAM (KIP-IN) 


NUMBER OF DIVISION IN A STORY 
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ROTO 
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SR 

SR 1] 
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SSB 
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NUMBER OF TIMES HORIZONTAL LOAD TO BE INCREMENTED 
LENGTH OF A SEGMENT OF A WALL IN A STORY (IN) 
SEGMENT HEIGHT FROM BASE OF WALL(IN) 

NUMBER OF STORIES 

VERTICAL LOAD ON COLUMN (KIPS) 

VERTICAL LOAD ON WALL (KIPS) 

VERTICAL LOAD ON STORY (KIPS) 

SWAY ROTATION OF RIGHT HAND BEAM 

A SUBROUTINE SUBPROGRAM 

STORY ROTATION OF FRAME 

RATIO OF MOMENT TO PLASTIC MOMENT CAPACITY OF WALL 
RATIO OF CURVATURE TO CURVATURE CORRESPONDING TO RMWP 
ROTATION OF WALL IN A SEGMENT (RAD) 

ROTATION OF WALL AT THE FLOOR LEVEL (RAD) 

JOINT ROTATION OF FRAME (RAD) 

JOINT ROTATION AT THE BASE OF COLUMN (RAD) 
STIFFNESS OF LEFT HAND BEAM (KIP-IN) 

STIFFNESS OF RIGHT HAND BEAM (KIP-IN) 

STIFFNESS OF COLUMN (KIP-IN) 

SHEAR AT THE ENDS OF LEFT HAND BEAM (KIP) 

SHEAR AT THE ENDS OF RIGHT HAND BEAM (KIP) 

SHEAR IN COLUMN IN A STORY (KIP) 

A SUBROUTINE SUBPROGRAM 

A SUBROUTINE SUBPROGRAM 

A SUBROUTINE SUBPROGRAM 

SLOPE OF THE SECOND BRANCH OF SHEAR WALL MOMENT - 


CURVATURE DIAGRAM 
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TOP MOMENT IN THE TOP OF THE TOPMOST SEGMENT IN 

A STORY (KIP-IN) 

TOP MOMENT IN A SEGMENT OF THE SHEAR WALL (KIP-IN) 
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WALL (IN) 
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B.3 FLOW DIAGRAM OF THE COMPUTER PROGRAM 
B.3.1 MAIN PROGRAM 


NO. OF PROBLEM 


NO. OF INCREMENT] ' 


O) 
CALL 
SR 


BAKA 

NO 

YES 
Cikis 12) 
0 


N 
AXIAL LOAD NO 
PRESENT? 


CONVERGENCE TEST 
FOR DEFLECTION WITH 
P-A EFFECT WITHIN 

SPECIFIED LIMIT? 


COMPUTE KORIZONTAL 


YES 
NO 
CALCULATE HORIZONTAL 
FORCES TO BE APPLIED 
@ TO THE SYSTEM INCLUDING 
P-A EFFECT 


FORCES ON WALL 
FOR CHECKING 
DEFORMATIONS 


SET IK 


A eee 
OR IX = 1 
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CALL 
SR 3 


SHEAR WALL YES 
MOMENT > PLASTIC 
MOMENT CAPACITY? 


CALL 
SR 1 


REDUCE THE 
MOMENT OF INERTIA 
OF THE PLASTIFIED 
WALL STRIP 
ACCORDING TO 
M-@ DIAGRAM 


SET eLiGi=..0 
NO INELASTIC YES 
CNL = Jd? > RANGE 
= Au le? 


aq YES 


DECREMENT THE 
HORIZONTAL LOAD 
TO MAKE THE 
FRAME ELASTIC 


NO 


INCREMENT 
HORIZONTAL 
LOAD 


OUTPUT STATEMENTS 
FOR LOAD, 
DEFLECTION AND 
ROTATION FLOOR- 
WISE FOR EACH 
INCREMENT 


SET NL = JJ 
HLI = HLIR 
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B.3.2 SUBROUTINE 'SR' 
B15 


INPUT 
STATEMENTS 
FOR DATA 


COMPUTE 
1. FLOOR HEIGHT FROM THE BASE 
2. STIFFNESS OF MEMBERS 

3. TOTAL VERTICAL LOAD 


OUTPUT STATEMENTS 

FOR DATA OF 
FRAME AND 

SHEAR WALL 


RETURN 
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B.3.3 SUBROUTINE 'BAKA' 


(ws) 


WALL-FRAME 
ITERATION 
BEGINS 


CALCULATE 
1. MOMENT i. AT ALL SEGMENTS 


2. ROTATION OF THE WALL 
3. DEFLECTION ii. AT FLOOR LEVELS 


YES 
YES 


NO 
RETURN 


YES 


Be 
fe 


APPLY CONVERGENCE 
FORMULA FOR ROTATION AND 
DEFLECTION 


CALCULATE 
1. STORY SWAY 
2. BEAM SWAY 


| CALL 
FRAME 


NO 


SET MOMENT (COMPUTED 
IN THE PREVIOUS 
LOADING CONDITION) 
EQUAL TO THE PLASTIC 
MOMENT CAPACITY 


MOMENT (AT 
POTENTIAL HINGE 
LOCATIONS) GREATER 
THAN PLASTIC 
MOMENT CAPACITY? 


NO 


MES 


CALCULATE MOMENTS 
AT OTHER POINTS 
USING MODIFIED 

SLOPE-DEFLECTION 

EQUATION 


CALCULATE MOMENT 
USING SLOPE-DEFLECTION 
EQUATIONS 
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CONVERGENCE 
TEST FOR 
DEFORMATIONS WITHIN 
SPEC TEICAL IML 


RETURN 


CALCULATE FORCE 
FOR THE NEXT 
CYCLE, 


OUTPUT STATEMENT 
FOR INSUFFICIENT 
NUMBER OF 
GYCEES 


RETURN 
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B.3.4 SUBROUTINE ‘FRAME! 


B18 


ITERATION FOR 
JOINT ROTATION BEGINS 


INITIALIZE ALL 
QUANTITIES 
“ CC INcREMENT NO = 12D 


YES 


MOMENT AT 
A POTENTIAL 
HINGE LOCATION 
EXCEEDS THE PLASTIC 

MOMENT CAPACITY 
OF THE MEMBER? 


MODIFY THE PARTICULAR 
JOINT ROTATION 
EQUATION TO TAKE INTO 
ACCOUNT THE FORMATION 
OF PLASTIC HINGE 
ACCORDING TO TABLE A.1 


YES 


CALCULATE ROTATION AT 
ALL THE JOINTS USING 
JOINT ROTATION 
EQUATION (A.1) 


VES 


NO 


CONVERGENCE 
TEST FOR JOINT ROTATION 
BETWEEN SUCCESSIVE CYCLES 
WITHIN A SPECIFIED LIMIT? 


eS 
RETURN 


OUTPUT STATEMENTS 
FOR INSUFFICIENT 
NUMBER OF 
CYCLES 
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B.3.5 SUBROUTINE 'ROFA' 


OUTPUT STATEMENTS FOR 
1. JOINT ROTATIONS 


3. STORY SHEAR AND FORCES 


B.3.6 SUBROUTINE 'SR3' 


OUTPUT STATEMENTS FOR 
hy & SLOPE 

2. ‘DEFLECTION 

3. MOMENT AND SHEAR 
IN THE SHEAR WALL 


RETURN 


2. - MOMENTS AT POTENTIAL HINGE 
LOCATIONS IN FRAME MEMBERS 
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Geen) SUBROUTINE —<SRi* 
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MOMENT AT A 
POTENTIAL HINGE 
LOCATION HAS EXCEEDED THE 
PLASTIC MOMENT CAPACITY 
OF THE MEMBER? 


RETURN 


OUTPUT STATEMENT 
INDICATING LOCATION 
OF PLASTIC HINGE AND 
THE MAGNITUDE OF 
MOMENT 


SET THE MOMENT AT 
HINGE LOCATION .0001 GREATER 
THAN PLASTIC MOMENT CAPACITY 


RETURN 
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B.4 LISTING OF THE PROGRAM 
MAIN B21 


FORTRAN PROGRAM FOR THE ANALYSIS OF SHEAR WALL-FRAME STRUCTURE 


DIMENSION F(99,30),HSF(30) 4, 0W(30) ,F0C( 30)» SHEAR(30) »y ROTF(99, 30), 
LROTFF(99,31),ROTO(99), SHEC{ 30) »FD(30) »FW{ 30) y SHEARW(30) »PS{ 30), 
2SHE ARR (30) »CPHI (300) ,D( 300) 9D1{300) sH(99y 30) yDELTA( 99,30) , 
3FA(30)_PPS( 31) ,CROTF(99430) sCDEFF(99,31), PPD{(30),ZM{300), 

4 DEFF(99,30), HBF (30) »HBW{30),HS( 31), SBF{30),SC (31) ,SBW(30), 
5F11(30)sTMCM(3200) »8MGM( 300) »TMOMF (30) ,BMOMF (30) yDMFB(30), 
6DMBF (30), DMBW( 30) »DMWB(30),CSHEC(30), DMCB{(31),DMCT(31) 

REAL NNM(300) »NNMM{300), MIA(300),MI1(300),MMCT( 30) »MMCB( 30), 
1KByKC, MOMW( 30), MOMWB (30) »MOMBW(30),MOMBF( 30) »MOMFB(30) 
2yMPF(30),MPW(30),MPC(31),MPS{( 300) 

COMMON IUNB 

20 FORMAT(1H1) 
21 FORMAT( 1HK) 
22 FORMAT{1H ) 

READ(5,130) MM 

30 FORMAT(1X,13) 
DO 400 JJ=1,MM 


INPUT STATEMENTS FOR DATA BY SUBROUTINE SR. 


BAL.1 SR (LyNSoNDy EWs KC sKBy FDyNI,HS,DWy PS, MPF ,NNM, 
1 MPWyMPC, SC, SBW,SBF,IVL,PPD, HBF »HBW, 
2 RMWP »RPWP » MAX ,CONygHLI,HLIRsHSFaFeMI Ly MPS aNNMMyMIA PPS, ALPHA) 
NL=JJ-1 

KP=0 


SSB=(RMWP-1.0) /(RPWP-1.90) 
DO 302 NJ=1,NI 


SHEAR WALL ANALYSIS BEGINS. 


tK=0 

IX=0 

MMMM=1 

DO 131 K=1,NS 
H{NJ,K) =F(1,K) 


31 CONTINUE 
32 CALL BAKA (IKy1X~MMNMyLyNS_NDyEWy KC ¥KBy HS y DWy PS yNNMy MPF yMPW, 


LMPC SC» SBW,y SBF yg HBWy MAX »CON,HSF 9 Fy MI Ty NNMMy TMOMy BMOM » MOMWs ROTB, TMOM 
2F, BMOMF ,ROTE, DEFF,ROTFFyRCTO,ITERyNJy DMCByDMCT, OMBFy DMFBy DMBW, 
3NMWByMMCBy MMCT » MOMBWy MOMW Ey FOC» SHEC oC SHEC » SHEARWy NCYCL Ey SHEAR» CDEF 
4F yCROTF,ZMyMIA) 

IF{IX .£Q. 19GO TO 245 

IF( IK .£&Q. 1) GO TO 272 sates 

BS IF(PPS(1) .LE. 0-0) GO TO 2 
SEANMMM 00-1) «AND. NJ «Qe. 1) GO 10 235 
[F(MMMM .EQ. 1) GO TO 256 


NO 234 N=1,NS 
TE CABSL (COCERUMMMM, ND —CDEFF (MYMM=1,N))/COEFF( MMMM» N) J »GE. CON .OR 


1. ABS( (CROTE(MMMMyN)-CROTF( MMMM—1yN))/CROTF(MMMMyN)) 2GE. CON) GO 
270.251 
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| B22 
| 234 CONTINUE 


GO TO 260 


OUTPUT STATEMENTS FCR RESULTS CF FRAME ANALYSIS BY SUBROUTINE ROFA 
(WITHOUT THE EFFECT CF AXIAL LOADS) 


| 235 WRITE(6,240) 

| 240 FORMAT(/45X,"RESULT WITHOUT THE EFFECT OF AXIAL LOAD®//) 

CALL ROFA {(ROTFF, SHEC yFOC »MMCB yMMCT yROTO, NS yNCYCLE, ITER, MMMM) 

DO 242 N=1,NS 

| MOMFB(N)=SBF(N)*6 .O#ROTFF(ITER,N) 

MOMBF{(N)= MCMFB(N) 

SHEARR(N)=2.0*MCOMFB(N)/HBF(N) 
WRITE(69241) NyMOMBW(N) »sMOMWBI(N) »MOMBF(N) »MCMFBI(N), SHEARW(N),SHEAR 
LR(N) 

241 FORMAT (13X_139 1X9 4F190295X%yFl 0025 7X9 F1002/) 

242 CONTINUE 

WRITE(6,21) 

DO 244 K=1,NS 

FL1(K)=F{1,K)-FOC (K) 

F(l,K)=F11(K) 

1244 CONTINUE 

| ix=1 

GO TO 132 

(245 WRITE(6,22) 


OUTPUT STATEMENTS FOR RESULTS OF SHEAR WALL ANALYSIS BY SUBROUTINE 
SR3 (WITHOUT THE EFFECT OF AXIAL LOADS) 


CALL SR3 (BMOMF,TMOMF,ROTB,NSsFIlyFOC,ySHEAR,ROTF,DEFF,CROTF,CDEFF, 
1MMMM) 

IX=0 

WRITE (6,21) 


ADDITIONAL HORIZONTAL LOACLTO SIMULATE P-DELTA EFFECT 


i251 IF( MMMM .EQ. MAX) GO TO 254 
256 PPS(NS+1)=0.0 
HS(NS+#1 )=5000.0 i 
CDEFF(MMMM,NS+1)=CDEFF(MMMM,N 
F(1,1)=HORJ) 1) #AL PH AX (PPS (1) *CDEFF( MMMM, 1)/HS(1)-PPS(2)%(CDEFF( MMM 
1M, 2)-CDEFF(MMMM,1))/HS(2)) 
TANS SEQs, 19460 TO; 253 
DO 252 N=2yNS 
F(1,N)= =H(NJgN) +AL PHA® (PPS (N)#(CDEFF (MMMM, N)-CDEFF (MMMM,N-1) DZ HSIN) 
1-PPS(N+1L) *(CDEFF(MMMMyN4+1)-CDEFF(MMMMyN)) JHSENt+1)) 
252 CONTINUE 
253 MMMM=MMMM4] 
GO TO 132 


OUTPUT STATEMENTS FOR RESULTS OF FRAME ANALYSIS BY SUBROUTINE ROFA 


254 WRITE(6,255) 
>55 FORMAT(10X,*CONVERGENCE (CALCULATION OF THE EFFECT OF AXIAL LOAD) 


IWAS NOT ENOUGH'//) . 


DO 279 N=1,NS 
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TF(MMMM.LEQ.MAX.OR.NCYCLE.EQs. MAX) GO TO 960 
IF(NJ EQ. 1) GO TO 261 
TF(ABS({OMFB(N)) .GE. MPF(N)) GO TO 262 


MOMFB(N )=SPFI(N) #6. 0*ROTFF(ITER,N) 


(271 


273 


‘83 


184 


eo T0 263 

MOMFB(N)=DMFB(N) 

MOMBF (N)=MCMFB(N) 

SHEARR(N)=2 .0*MOMFBIN) /HBFE(N) 


WRITE (69241)Ny MOMBWIN) »MOMWB(N) »MOMBF(N) + MOMFBIN) » SHEARW(N) »SHEARR 
1L{N) 


(270 CONTINUE 


DO 271 K=1,NS 
DELTA(NJ,yK)=CDEFF (MMMM,K) 
F1l(K)=F(1,K)-FOC (K) 
F(1,K)=FL1(K) 

CONTINUE 

IK=1 

ome TO! 132 

WRITE(6,21) 


DUTPUT STATEMENTS FOR RESULTS OF SHEAR WALL ANALYSIS BY 
SUBROUTINE SR3 


CALL SR3 {BMOMF »TMOMF, RCTByNS yF11LyFOC,SHEAR,ROTFy DEFF,CROTF yCDEFF, 
1MMMM) 
WRITE(6,21) 


DETECTION OF HINGES IN SHEAR WALL. 


PRE NECGNESOJJ) GO TO 273 
IF{ABS(DELTA(NJ,NS)) .GT. 
SSB=(RMWP-1.0)/(RPWP-1.0) 
DO 284 J 1yt 
D(J=(TMOM(J)+#BMOM(J))/{(2.0*MPS(J)) 

DiC J)=D(J)*MPS{(J) 

IFCABS(D(J)) .LT. 1.0) GO TO 284 

IF(NUL .NE. JJ) GO TO 274 

CPHI( J) =(ABS{D( J) 1-1.9)/SSB+1.0 

MII(JI=CABS(D(J)) ) *MIA(C I) /CPHI(J) 

KP=K P+] 

IF{(KP .GE. 2) GO TO 281 

WRITE(6,280) 

FORMAT ( 45X,37HDETEC TION OF HINGES IN THE SHEAR WALL///) 

; D(J),D1(J) 

ee ONCE: IN SECTION NOw2I392X_y51HRATIO OF MOMENT IN WALL 
1 TO PLASTIC MOMENT CAPACI TY=,F5.232Xys15HMOMENT IN WALL=3F14.2/) 

J),MIACJ) 

Reyer puceD MOMENT OF INERTIA =", F13.2,3Xye* INITIAL MOMENT 
1 OF INERTIA =!,F13.2//) 

CONTINUE 

WRITE(6,21) 


STOPD) GO TO 960 


— 


DETECTION OF HINGES ON FRAME BY SUBROUTINE SRI. 


(DMF B, OMBF »OMBs DMWB sy DMCT, OMCB yMOMFB y MOMBF » MCMBW,MO 


CALL SR1 


1MWB,yMMCT »yMMCBy NSy M39MPF yMPW MPC) 


IK=0 
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TrANL CEC. JJ) GO TC 285 
IF(KP .GE. 1 .OR. M3 .GE. 1) GO TO 295 


HORIZONTAL LOAD IS INCREMENTEC 


IF{NJ -EFQ. NI) GO TC 960 

WRITE{(6,21) 

WRITE(6,291) 

FORMAT(50X, "HORIZGNTAL LOAD INCREMENTED*®//) 

WRITE(64292) 

FORMAT(10X,*FLOGR NCw*,5Xy*HORIZONTAL LOAD(K)*,5X_"*VERTICAL 


1K)"//) 


DO 294 K=1,NS 

F(1lsK)=H(NJ,K) tHLI*FD(K) 

IF{IVL .LE. 0) GO TC 296 
PPS(KJ=PPS(K)+HLI¥*PPD(K) 
WRITE(6,293)KyF{(l eK) yPPS(K) 
FORMAT (13X913,11XyF8.2516X9F8.2/) 
CONTINUE 

KP=0 

mo 70 302 


HORIZONTAL LOAD IS CECREMENTED 


COMPUTE SECTION HEIGHT FRCM BASE IN INCH 


IF{INJ ~EQ- NI) GO TO 960 
STOPD=90.0*ABS (DELTA(NJ,NS)) 
WRITE(6,300) 


WRITE(6,292) 

DQ 301 K=1,NS 

F(1l,K)=H(NJ,K) —-HLI*FD(K)+HLIR*FD(K) 
IF(IVtL .LE. 0) GO TC 303 
PPS(K)=PPS{(K)-HLI *PPD(K )+HLIR*PPDIK) 
DMFB(K)=0.0 

DMBF(K)=0.0 

DMWB{K)=0.0 

DMBW(K)=0.0 

DMCT(K)=0.0 

DMCB(K)=0.0 

WRITE(65293) KyF{lsK)y PPSIK) 
CONTINUE 

HLI=HLIR 

NL=JJ 

KP=0 

CONTINUE 


OUTPUT THE ROTATION AND DEFLECTION OF EACH HORIZONTAL LOAD 


DO 907 N=1,NS 
WRITE{6,915) N 
FORMAT(48X%> *LOAD-DEFLECTICN DATA FOR FLOOR NO. 'y13//) 


WRITE(6,905) 
Bear cox, * WOR IZONTAL LCAD(KIP)*,5X_, "DEFLECT ION(IN)*//) 


LOAD ( 


300 FORMAT(1LHK,15X,"*HORIZONTAL LOAD DECREMENTED TO MAKE THE FRAME ELAS 
1ITIC SO THAT HORIZONTAL LOAD CAN BE INCREMENTED SLOWLY*///) 
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ITE(65903) H(KyN) ,DELTA(K,N) 
Bere MAT (16X%yF8.2,11X,F£13.6/) 
CONTINUE 

CONTINUE 

) CONTINUE 

peoTorP 

END 
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SUBROUTINE SR READS IN DATA - COMPUTES CO-ORDINATES OF EACH 
FLOCR FROM THE BASE AND THE STIFFNESS OF MEMBERS — ALSO ALL 
THE INPUT QUANTITIES ARE PRINTED OUT 


SUBROUTINE SR {(LyNSyNDy EWyKCyKByFDyNIyHSsDWyPS> MPF yNNM, 

1 MPW,MPC, SCySBW,SBFyIVL»PPD, HBF »HBW, 

im 2 RMWP y»RPWP MAX, CONyHLI gHLIR yHSFa Fy MII 9 MPS pNNMMyMI AyPPS, ALPHA) 

| DIMENSION FD(30),HS (30) ,DW(30 ),PS{(30), SC(30); 

| 1 SBF(30),SBW{30), HBF(30), 
2HBW (30) »HSF(30),F 199,30), PSW( 30), PPS{ 30) » PPD( 30) 

| REAL MI(30),MPF (30) »MPW(30) »MPC( 30) »MPSW( 30) »MIC(30),MIBF{ 30), 

| 1 MIBW(30) »NM(30)_,NNM(300) MIT (300),MPS( 300) yNNMM(300)yMIA (300), 

| » 2KB,KC 

COMMON IUNB 

READ(5,y200)KBy KC, EF yEWgNS gND a MAX »CONYHLIyNI,HLIR 

ago FORMAT (1X,2E18 5, 2F 7.053139 2F 502,13 9F 502) 

READ(5,470) RMWPyRPWP,ALPHA,IVL,IUNB 

470 FORMAT(1X_2F7. 29F 522,213) 

| DO 202 K=1,NS 

READ(5, 204) FD(K),HS(K) »HBFE(K) »HBW(K) »DW(K) sMI(K),MIC(K), 

| IMIBF(K) »MIBW(K) 

04 FORMAT (1LX,F5.2,4F70294Fll 02) 

202 CONTINUE 

BaTiO? K=1.NS 

: READ(5,LO1L)IPS(K) ,PSWI{K) 

101 FORMAT(1X,2F8.2) 

92 CONTINUE 

ie 6205: 5 AK as cli NS 

READ(5,500) MPF(K) »MPW(K),MPCIK) »MPSW(K) 

30 FORMAT ( 1x GAP i4elsh2i.c) 

)5 CONTINUE 


ALL FORCE UNITS ARE IN KIPS AND ALL LENGTH UNITS) ARE 
IN INCHES UNLESS STATEC OTHERWISE. 


L=NS*ND 

NDN=0 

K=1-ND 

CB=0.9 

NO 22 J=1lsNS 

NM( J) =HS(J)/FLOAT(ND) 
HSF(J)=CB4HS{(J) 
CB=HSF{J) 
SBW(J)=ER*MIBW( J) /HBWU J) 
SBF(J)=EF*XMIBF (J) /HBF(JI) 
SC{ J)=EF*MIC(J)/HS(J) 
F(1,J)=FD(J) 
PPS(J)=PS(J)+PSW(J) 
PPD(J)=PPS(J) 

K=K4N) 
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53 
510 
501 

364 


aT 


358 
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NDN=NON4ND 

NO 23 N=K,NON 
NNM(N)=NM(J) 
MIT{NI=MI (I) 
MIA(N)=MI(J) 
MPS(N) = MPSW{ 
CONTINUE 
CONTINUE 
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Sumrure “SECTION HEIGHT FRCM BASE IN INCH 


3C=0.0 

DQ 25 K=1,L 
NNMM(K)=BC+NNM 
BC=NNMM(K) 
CONTINUE 
WRITE(6,324) 
FORMAT(1H1) 


{K) 


Ber erul’ “STATEMENTS” FOR’ YEAVA OF FRAME ANC SHEAR WALL. 


WRITE(6,51) 


FORMAT(50X,29HDATA FOR SHEAR WALL AND FRAME//) 


WRITE(6,351) NS 
FORMAT(1LOX,56HNUMBER OF STORIES = 


Lt 713/) 
WRITE(6,352)ND 


FORMAT(10X,56HNUMBER OF DIVISION TO BE MADE IN A STORY 


1 713/) 


WRITE{6,355) EW 
FORMAT(10X,53HMODULUS OF ELASTICITY FOR WALL = 


i hee AES ae) 
WRITE(6,356) E 


FORMAT(10X,53HMODULUS CF ELASTICITY FOR FRAME = 


1 3F9.2/) 


= 


WRITE(6,354) KC 


FORMAT(10X%,43HSPRING CONSTANT AT BASE OF FRAME = ,€19.5/) 
WRITE(6,353) KB 
FORMAT(10X%,43HSPRING CONSTANT AT BASE OF WALL = ,£&19.5/) 


WRITE(6,510) C 


ON 


FORMAT{10X_*CONVERG ENCE LIMIT *,24Xy_ '=',F20.5/) 
WRITE(65501) ALPHA 


FORMAT(10X," AL 
WRITE(6,5364) 
FORMAT(1HK) 
WRITE(6,357) 


FORMAT(10X,1LOHFLOCR NO.- 
13X,y16HSTORY HEIGHT(IN)» 


PHA", 36X_%=",F20-3//) 


/y2X%y1lOHFORCE(KIP) »3Xy17HMOMENT OF INERTIA, 
3X,14HWALL WIDTHCIN) »3Xy15HAXIAL LOAD(KIP), 


23X91 5HAXIAL LOAD(KIP)) 


WRITE(6,358) 


FORMAT{10X,9HSTORY NO-y18Xe 


20N WALL//) 
DO 10 N=1,NS 
WRITE(6,359)Ny 
FORMAT(13X_913% 
CONTINUE 
WRITE(6,366) 


L3HOF WALL (1N4)943X_,9HON COLUMN,1LOX,7H 
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266 FORMAT(1HK) 
WRITE(6,369) 
360 FORMAT(10X,10OHFLOCR NOe/92Xy14HLENGTH OF REAM,3Xy14HLENGTH CF BEAM 


Ly3X,1L7HMCMENT OF INERTIA, 2Xy17HMOMENT OF INERTIA,3X_17HMOMENT OF I 
2NERTIA) 


WRITE(6,361) 
361 FORMAT(10Xy9HSTORY NO. 24Xy13HCN ROLLERCIN) » 3X, L4HCONNECTED WITHy4X 
Lyl4HOF COLUMN(IN4),5X,17HCF BEAM ON ROLLER»y3Xy17HOF BEAM CONNECTED 
2) 
WRITE (65362) 
362 FORMAT{(42Xy 8HWALL (IN) 532Xy5H(1N4),10X,14HWITH WALL(IN4)//) 
DO 71 J=1,NS 
WRITE (64363) J5y HBF (J ),HBW( JS) pMIC(J),MIBFE(J),MIBW(J) 
363 FORMAT (13X_13 5 8XyF8 02 99XoFBo2 sIXe FLL o299XsF ll o299Xe Flle2/) 
71 CONTINUE 
WRITE(65411) 
411 FORMAT(1HK,30Xy"PLASTIC MCMENT CAPACITY OF BEAMS, COLUMNS AND WALL 
els 7) 
WRITE(6,399) 
399 FORMAT(1CX,10HFLOCR NO./92X;5 
1 LSHPLASTIC MOMENT,2X,15HPLASTIC MOMENT,2X,15HPLASTIC MO 
3MENT, 5X,15HPLASTIC MOMENT) 
WRITE(6,406) 
406 FORMAT(10X,9HSTORY NO. 93X,y 


1 L5HCAPAC ITY OF THE,2X,15HCAPACITY OF THE »,2X,1L5HCAPACITY OF T 
2HE,y 5X, 15HCAPACITY OF THE) 
WRITE(6,407) 


| 407 FORMAT(22X, LSRBEAM 


20N ROLLER,2X,15HWALL SIDE BEAMs2X,15HCOLUMN (KIP-IN)» 5X_15HWALL 
3 {KIP-IN).) 
WRITE{6,408) 


1403 FORMAT(26X; BHIKIP-IN 


1), 7X,8H(KIP-IN)//) 
DO 410 K= 1,NS 


| WRITE(63409) Ky MPFCK) »yMPW(K)eMPC(K) »MPSW(K) 
} 409 FORMAT(13X,13_ 6X, F142 2X9F14e25 3XyF14.252XeF18 
1.8/) 


410 CONTINUE 
WRITE(6,385) 

385 FORMAT{(1LHK,45X,39HSTIFFNESSES (£I/L) OF BEAMS AND COLUMNS//) 
WRITE(6,365) 

365 FORMAT(10X,1OHFLCCR NO./,4Xy1L9HSTIFFNESS OF COLUMN, 4X, 17HSTIFFNESS 
1 OF BEAM,4X,17HSTIFFNESS CF BEAM) 
WRITE (6,367) 

367 FORMAT(10X,9HSTORY NO« yLOX,8H(KIP-IN) »11Xy15HCONNECTED WITHs5X917 
LHON ROLLER(KIP-IN)) 
WRITE(6,368) 

368 FORMAT(49X,12HWALL(KIP-IN)//) 
Ba 9720)=1,NS 
WRITE( 65369 )Jy SC(J) oSBW(J) 2SBF(J) 

369 FORMAT(13%_13y9XeF 15029 7X 9F 156296X9F15.2/) 

72 CONTINUE 
WRITE(65481) 

481 FORMAT(1HK, 40X,53HMCMENT CURVATURE RELATICNSHIP OF THE SHEAR 
1 WALL///) 
WRITE(6,482) 


482 FORMAT(10X,9HPCINT NOes5Xy19HRATIO OF MOMENT IN, 5Xy18HRATIO OF CU 
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WRITE (5,483) 
FORMAT(24Xs1L9HWALL TO THE PLASTIC,5X,1L8HIN WALL TO tHe? 
WRITE(6,484) 
FORMAT(24X,19HMOQMENT CAPACITY OF,5X,18HCURVATURE OF THE) 
WRITE (65485) 
FORMAT (29X,9HTHE WALL.y1LOX,18HWALL AT YIELD PT.///) 
WRITE (6,488) 
FORMAT(13X,_"NO.1 1.00 1.-0C'/) 
WRITE(6,487) RMWP,RPWP 
FORMATLI3Xs*NO.2* pL 5XeFT C291 BXyFT 2/9 
WRITE(6,324) 
RETURN 

END 
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SUBROUTINE BAKA COMPUTES THE CISTRIRUTION OF LATERAL LOAD BETWEEN 
THE FRAME AND SHEAR WALL. MOMENTS AND DEFORMATIONS AT EACH SEGMENT 
QF THE WALL ARE COMPUTED. THE CONVERGENCE FORMULA IS APPLIED 

( EXCEPT THE FIRST CYCLE ) FOR DEFLECTION AND ROTATION. THE 
NEFCRMATICNS COMPUTED BY THE CONVERGENCE FORMULA ARE ENFORCED ON 
THE FRAME SYSTEM. THE JOINT ROTATION CF THE FRAME ARE COMPUTED IN 
SUBROUTINE FRAME, 
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SUBROUTINE BAKA (IKsIX+MMMMsL yNSyNDyEWy KC KB yHSyDWy PSyNNMy MPF yMPWy 

IMPC ySC 4 S3Wy SBF yHBWs MAX ¢ CON ,HSFy Fy MI Ty NNMMy TMOMy BMOMy MOMWy RCTBy TMOM 
2F,BMOMF »RCTF,DEFF,ROTFFsRCTOs ITERyNJy DMCByDMCTy DMBF,DMFB,DMBW, 
3DMWBsMMCB,MMCT y MOMBWsMOMWE», FOC, SHEC») CSHEC y SHEARWy NCYCL Ey SHEAR» CDEF 
| 4F,yCROTF) 
. DIMENSTON BMOM(300 ),BMOMF{( 30) yCSHEC( 30) ,CNEFF( 99430), CROTF(99,30 
| 1) ,DEF(300),DEFF (99, 30) »9MCB{ 30) ,NMCT( 30) »DMBW(30) »DMWB{ 30) ,OW(30) , 
| 2DMFRB( 30), DMBF(30),F(99,30),FOC( 30) yHSF( 30) »HS(30C) sHBW(30) yPS( 30), 
3ROT (300), ROTFI99,30),ROS (30) »yROB{ 30), ROTFF (99,30), ROTO(99),5C(30) 
| 4ySBW(30),S8F (30), SHEC(30) »SHEARW( 30) » SHEAR ( 30) » TMOM( 300) , TMOMF (30) 
| 5,VNFF (30) 

REAL KB,KC,MOMB,MOMW{30) MII (300),MPC( 30),MOMP{ 30) ,MPW( 30), 

IMPF (30) »MMCB(30),MMCT (30) »MOMBW(30) » MOMWB ( 30) »NNMM( 300) >NNM(300) 
| COMMON IUNB 
\gees2 DO 232 M=1, MAX 


' MOMENTS AT VARIOUS SECTIONS (KIP-IN). 


MOMB=0.9 
| NDN=0 
K=1-ND 

DO 135 N=1,NS 
MOMB=MOMB4E(M,N)XHSFI(N) 
K=K4ND 
NDON=NDN4ND 
DO 134 J=KyNDN 
RS=0.9 
DO 133 I=N,NS 

133 RS=RS4+F(M,1I9*{ HSF CI I-NNMM(J)) 
TMOM(J)=RS 

1324 CONTINUE 

135 CONTINUE 
=e SECS IVOR IER PLEQ. 1) GO TG 140 
IF(M .F€Q. 1) GO TO 145 

140 DO 142 N=1,NS 
MOMB=MOMR4+MOMW(N) 
RS=0.0 
NO 141 I=N,NS 
RS=RS#+MOMW( 1) 

141 CONTINUE 
MOMW(N)=RS 

142 CONTINUE 
NON=0 
K=1-ND 
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DO 144! .1=V NS B31 
K=K+4N9 

NON=NON4ND 

NN 143 J=K,NDN 

TMOM( J)=TMOM(J) #MOMWCT) 

RMOM( JS) =TMOM(J) #SHEARCT)*®ANM( J) 
CONTINUE 

CONTINUE 

GO TO 150 

BMOM(1)=MOMB 

mex 0.EO 1) GO To 150 

DO 146 J=2,L 

BMOM(J)=TMOM(J-1) 


GeeLECTLONS AT \VARTCUS SECTIONS (IN) 


A = (TMCM(1L)I+MOMB)FANM (1) /02.0*EWSMITI(1)) 
ROTB=MOMB/KB 

ROT(L)I=ROTB+A 
DEFCLI=ROTB*NNM(L)#A*¥NNM( 1072.9 

Prete ieGGe (1) 160 TO 1s2 

90 151 J=2,t 

B=(TMOM(J1+BMOM(J)) #NNM(JS)/(2.0*EW*MIT( J) ) 
WOT 099 t= "ROT(J-1) + 8B 

NEF (J) =DEF( J-1 9 4#ROT (J-1) *NNM( J) #B*NNM(J1/2.20 
CONTINUE 


MOMENTS, ROTATICNS AND CEFLECTIONS AT EVERY FLOOR LEVEL 


DO 153 N=1,NS 

T=N*ND 

K=1+1-N9) 
TMOMF{N)=TMOMCUT) 
BMOMF(N)=BMCM{K) 
ROTF(M,NI=RCT(T) 
DEFF(MyN)=DEF(T) 
CONTINUE 
DRCEXISEQ@s).1) RETURN 
IFC IK «EQ. 1) RETURN 
PrieJd.6G0 <1) 'GOOTO 155 


DQ 154 N=lyNS 
ROTF(M,N) =ROTE( Ly N) #ROTE(M—1yN) / (ROT (M-1LyN)-ROTF(MyN) 


DEFF(M, N)=DEFF (1,N) #DEFF(M-1,N)/ (DEFF(M-L»N) -DEFFIMyN)} 
CONTINUE 

ROS(1)=DEFF(M,1)/HS (1) 

PP(NS SEQ. 1) GO TC 161 

90 160 T=2,NS 
ROS(I)=(CEFF(M, 1) -DEFF(M, 1-1) )/HS(T) 
CONTINUE 

00 162 N=1,NS 

VDEF(N) =-ROTF (M,N) *OW(N)/ 2.0 
ROBI{N)=VDEF(N) /HBWIN) 

CONTINUE 


COMPUTE MOMENT, SHEAR AND FORCE ON FRAME BY SUBROUTINE FRAME 


CALL FRAME (MAX,ROTFF»ROTO,SC »ROS sKCe SBW, SBF »sROBsROTF yNSy ITER, MMMM 
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bes My DME By NMC Ty OMBWyDMWBy EMRE yDMEB MPC, MPF, MPW, NJ CON) 
PAR.) FSC HELY (90 sTOMRRO 

PRCABSGAOMCA(1)) .GE. MPC(1)) GO TN 181 
MMCBR(L)=-KC*XRNTOCITER) 

GY 210, bar 

MMCB(1L)=0MCB(1) 

MMINGNSEO. 1) 'GO TO 183 

PROABS(OMCT(1)) ~GE. MPC(1)) GO TO 184 

PAWABSTOMORA)) o.GE. MPC(1)) GO TO 185 

MMC T(1)=SC(1)*{(4.0*ROTFFCITER,1)4+2.0*ROTO(ITER )-6.0*ROS(1L)) 
GO TO 192 

MMCT(1)=OMCT(1) 

Go 19 192 

MMC T(1L)=SC(1L)*(3.0*#ROTFF(ITER,1)-3.0*ROS{1))40.5*DMCB(1) 
SHEC(1)=(MMCT(1)4#MMCB(1))/HS(1) 

PRUNGD LEO, Lv 1GG TO)210 

Pa 1205 “= 25:NS 

Svein EQ. 1) GO To 195 

IF{ARBS(OMCB(K)) .GE. MPC(K)) GC TO 196 

IFC ABS(OMCT{K)) .GE. MPC(K)) GO TO 197 
MMCB(K)=SC(K)*(4.0¥*ROTFE( ITER yK-1)4+2.0*ROTFF(ITER,K)I-6.0*ROS(K)) 
GO T0 198 

MMCB(K)=DMCB(K) 

GO 10 198 

MMC 8(K)=SC (K) *(3.0*ROTFF( ITER yK-1)-3.0*ROS(K) )40.5*DMCT(K) 
iFUND £0. 1) GO Te 199 

IF{ABS(CMCT(K)) .GE. MPC{K)) GO TO 200 

IF(ABS{DMCB(K)) ~GE~. MPC{K)) GO TO 201 

MMC T(K)=SC(K)*(4.0*ROTFFCITER ,K)+2.0*ROTFF( ITER »K-1)-6.0*®ROS(K)) 
GO 70) 2.02 

MMCT(K)Y=DMCT(K) 

GOvToO202 

MMCT(K)=SC(K) * (3. O¥ROTFE( ITER ,K)—-3.0*®ROS(K)14+0.5*DMCR(K) 
SHEC(K)=(MMCT(K)4#MMCB(K))/HS(K) 

FOC (K-1 )=SHEC(K)-SHEC(K-1) 

CONTINUE 

FOC (NS) =-SHEC(NS) 

SHEAR(1)=0.0 

NO 219 K=1,NS 

PRINd ~wEQ.) 1) GO TG 211 

IF{ABS(DMBW(K)) .GE. MPW(K)) GC TO 212 

IF{ABS({ DMWB{K)) .GE. MPW(K)) GO TO 213 

MOMBW(K)=SBW(K) *(4.0#ROTFFLIT ERsK)4+220*ROTFIM,K)-6.0#ROBIK)) 
Go TO 214 

MOMBW(K)=DMBW(K) 

GO TO 214 
MOMBW(K)=SBW(K) *( 3. O0#RCTFF(ITERsKI—320*ROBIK) 14+0.5*DMWBI(K) 
Pend? 5 EO. 1) GOTO 215 

IF(ABS( DMWB{K)) GE. MPW(K)) GO TO 216 

IF(ABS(DMBW(K)) .GE. MPW(K)) GO TO 217 

MOMWB{(K)= SBW{K)*{2 OXROTFF(ITERsK) +4.0¥*ROTF(M,K) -6.0*ROBIK)) 


BO 10 218 
MOMWB(K)=DMWB(K) 


60: T0 218 
MOMWB(K)=SBW(K) *( 3 O*ROTF (MyK) -3-0*ROB(K) ) +0.5 *OM BW IK ) 


SHE ARW(K)=( MOMBW(K) +MOMWB (K)) /FBW(K) 
MOMW(K) =-MCMWB (K) -SFEARW(K) *DW(K)/2.0 
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B33 
F(M+1,K)=-FCC(K) 


| SHFAR(1)=SHFAR (1) 4F(M4#1,K) 
| 219 CONTINUE 
| PFINSSSEOL! 19 GO TD. 227 
| D9 226 K=2,NS 
| SHEAR (K )=SHEAR (K-1)-F(M4#1,K-1) 
| 226 CONTINUE 
Mer iF (M .EG. 1) GO Ta 232 
IF(M .EQ. MAX) GO TC 221 
DO 220 N=l1,NS 
TF (ABS{ (OEFFIM,N)-NEFF(M—-1,N) )/DEFF(M,N)D) «GE. CON «OR. ABS((ROTF( 
1M,N)-ROTE(M-1,N))/RCTF(M,N)) «GE. CON) GO TO 232 
220 CONTINUE 


6710) 223 

| 221 WRITE{6,222) MMMM 

222 FORMAT(10Xe"CYCLE NO.w=",13,7X_"*CONVERGENCE (CALCULATION ON WALL) W 

LAS NOT ENGUGH'/) 

| 223 NCYCLE=™ 

SHEAR(1)=0.0 

DO 224 N=1,NS 

| CDFFF(MMMM,N)=DEFF(NCYCLE WN) 
CROTFE( MMMM, N)=ROTF(NCYCLE »N) 
IF (MMMM .EQ. 1) CSHEC(N)=SHECI(N) 
SHEAR (1 )=SHEAR(1) +F(1,N) 

224 CONTINUE 

IF(INS .EQ. 1) GC TO 230 

DO 225 K=2,NS 

SHEAR (K )=SHEAR (K-1)-F(1,K-1) 

CONTINUE 

230 DO 231 K=1,NS 

SHEAR(K )=SHEAR (K) +SHEC(K) 

231 CONTINUE 
RETURN 

232 CONTINUE 
RETURN 
END 
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| IN SUBROUTINE FRAME, THE JOINT ROTATIONS CF THE FRAME FOR A 

| SWAYED POSITION (ENFCRCED BY THE WALL ) ARE COMPUTED. 

THIS IS PERFORMED BY GAUSS-SEICDEL ITERATICN METHOD. IF A HINGE 
FORMS IN THE STRUCTURE THE JOINT ROTATION EQUATION IS MODIFIED. 


SUBROUTINE FRAME (MAX, ROTFFyROTO,SC yROS yKCy SBWy SBFy ROByROTE yNS, ITE 
1RyMMMM,M,DMCBy DMCT,OMBW,DMWBy COMBE, DMFBy MPC, MPF, MPWy NJ yCON) 
DIMENSION ROTFF(99,31),ROTO(99),SC(31),ROS(31) » SBW{ 30) 2 SBF( 30) ,ROB 
1(30),ROTF (99,30), DMCB(31) »DMCT(31)»DMBW( 30) »DMWB( 30)» DMBF( 30), 
| 2DMFB(30) yAA(30) BB 30) CC (30) »DD( 30) gEE(30) » FF (30) »GG( 30) y HH( 30148 
| 3A(30) 
| REAL KCyMPC(31) »MPF(30),MPW( 30) 
COMMON IUNS 
| 90 174 T=1,MAX 
| fet 2671 1). GO T0164 
| DO 163 K=1,NS 
| ROTFF(1,K)=0.0 
| 163 CONTINUE 
(164 N=I 
PRINT EO.) 1) GO TO 165 
N=N-1 
165 SC(NS+1)=0.0 
ROS {NS#1)=ROS(UNS) 
| ROTFF(N,NS#1L)=ROTFFIN NS) 
| MPC({NS#+1)=500.0 
DMCB(NS#1)=0.0 
DMCTINS#1)=0.0 
IF(NS .EQ. 1) GO TO 10 
IF(ABS{OMCB(1)) .GE. MPC{1)) GO TO 20 
IF(ABS(DMCT(1)) .GE. MPC(1)) GO TO ll 
10 ROTOC(I) =(6. O¥SCUL )*ROS(1)-2.0*SC( 1) ROTFF UN, 19/04. 0*SCC1) #KC) 
GO TO 20 
11 ROTO(I)=(3.0*SC (1) *ROS(1)-0.5¥*DMCT(1))/(3 -C#SCULIFKC) 
20 90 170 J=1,NS 
PRUNU?2EQ.1 d45G0 TO 21 
IF({ABS(DMBW(J)) «GE. MPW(J)) GO TO 30 
IF(ABSCDMWB(J)) «GE. MPW(J)) GO TO 40 
21 AAC JI=SBW(J)*( 6.0*ROB(J)-2.0*ROTFIMSJI)) 
EE( J)=4.0*SBW( J) 
GO TO 50 
3C AA(J)=-DMBW(J) 
EE(J)=0.0 
GO TO 50 
40 AA( J)=3~.0*SBW(J)*2CB(J)-0.5*DMWBE J) 
EE(J)=3 .0*SBWI J) aE 
BO TPUNJ .£Q. 1) GO T 
IF (ABS{DMBF(J)) »GE. MPF(J)) GO TO 60 
LFE(ABS(DMEB(J)) «GE. MPF(J)) GO TO 70 
51 BBiJ)=0.0 
FF(J)=6.0*SBF(J) 
GO TO 80 
60 BB( J)=-DMBF (J) 
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B35 
FF(J)=9.9 


| fo 10 89 

| 70 BB(J)=-9.5*DMFB(J) 

| FF( J)=3.0*SBF( J) 

le fF(NJ .£O. 1) GO To 81 

MEP UASSIDMCT(J)) «GE» MPC(J)) GO TO 90 
IF (ABS(DOMCB(J)) .GE. MPC(J)) CO TO 100 

mreir (iJ) EQ. 1) GO TO 82 

COU SI=SCHA J941 68.0*ROSAd) —2. OFROTER(I J—1)) 

GO TO 83 

. 


| 82 CC(J)=SC( J) *(6.0*R0S(5)-2.0*ROTO(T)) 
83 GG(J)=4.0*SC(J) 
GO TQ 110 
| 90 CC{J)=-DMCT(N) 
GG(J)=0.0 
Sa) To10 
1090 CC(J)=3.0*SC(I)#ROS{(JS)—-0.5*DMCBIJ) 
GG(J)=3.0*SC(J) 
PO TFUNUWEQ.))1). GO TO 221 
IF(ABS(DMCB(J4+1)) .GE. MPC{J+1)) GO TN 120 
IF(ASS(CMCT(J+1)) »GE. MPC(J+1)) GO TO 130 
111 DD( J)=SClJ+1)* (6. O¥ROS (541) -2 .OXROTFFIN  Jt1)) 
HH( J) =4.0*SC(J41) 
GO TO 1490 
120 DD(J)=-DMCB{(J+1) 
HH{ J)=0.0 
GO T9 140 
130 DD(J)=3 .0¥SC(JS4+1) #ROS(J+1)-0. 5*OMCT (S41) 
HH( J)=3.0*SC(J41) 
140 BA( J)=FE( J) 4FF (5) 4+GG(J)4HH(S) 
IF (ABS(BA(J)) .LT. 0.0001) GO TO 170 
ROTFFE(I,J)=(AA( J) BB (J) 4+CC( J) 4001) /BA(J) 
170 CONTINUE 
171 ITER=! 
PEVLIGR (EQ. 29GO TO174 
) IF UABS(OMCB(1L)) »-GE. MPC(1)) GO TO 176 
) IF( IUNB.EQ.1)GO TO 176 | 
j LF (ABS ((ROTO(L)-ROTCLI-1L) )/ROTO(1)) .GT. CON ) GO TO 173 
eee one Bins 0001) GO To. 172 
IF (A IACKO)) «LE. C- Gg Te 
ee ion, c-ROTFR CIEL, 4) 17901 EPP ROD GT. CON 160 TO 173 
172 CONTINUE 
RETURN 
mis TF RL Seas (MAX GO TO 175 
174 CONTINUE 
175 WRITE(6,181) M,yMNMM 
181 FORMAT(10X,*CYCLE NOo="e1I322Xo'T 
1 ON FRAME) WAS NOT ENOUGH S/) 
RETURN 
END 


O%,1397Xy "CONVERGENCE (CALCULATION 


boy de ~ 


of AT ORs widest IMG 
OOF (OT Bo: CL 20m enya ran aC 
; S@ OF 0d tf .Os, 

(tite Dbiinic hentia cet. 


ae Cue id 
tl PORIARO SS OU 


Lpaamaya .0- (LI 2088 (L)9; 
(tp | 
112 (01 09 41 if 

OS1 OF O© C(f4L IIWM 23D. 41 eb pas 
O£1 OT OO (ULL ID9IM 13d. VLEs 
(C1+L WD I9TOROO. SELL 20880 ca) mELeL dE 
(14 UL) 32*Q. det 
A r ye < 


7 
trend err 


ITG ; 


; CLG Ts O-( 1 +L 1208" ( ve 
(1#uyyes 
ore tteehaests 
OT1 OF OO (2000.0 oTJde LE 
(LIUAGNELL 1OO* (L192 (07864 (LAA = Chee 
ws TTC 
aa eo 4: ed [2 
err Ot, OCT 03s 
ovr OT MD CHLISIM 239, CE LPHSR 

ats OT Odtt. 0 
EYE OT OD 4 WOD .TO. LITIOTORNG ET~ Reda ik bl 
oa 
STE OF 08 12O00.9> ease Cone ¥ 
VP OT BOC 409. Fo. ttm tT DIATOR NEM ET AOS SFA 6 1 1997 et 
pRRBE Mh, Hoe fee. HH 
Wits, 
ers 07 92. (XAM 


cil | ee ae 
TAIUDIAD) FIVIDAIVUND® .XT CEL, "OT eKSyETe* =s0M 3IIVI* KORY 
BRT POUNEA TAM’ 2AM 

4 v4 i wh, cite © 

; ‘ ; f ) 

can 

Q: 
# de! 


i Ba 9 I se: 


Jae St al ae lem SS 


'ROFA' B36 


OUTPUT STATEMENTS FOR FRAME FORCES AND DEFORMATIONS 


SUBROUTINE ROFA (ROTFFs, SHEC sFOC sDMCByDOMCT, ROTO,NSyNCYCLEy ITERy MMMM 
1) 


DIMENSION ROTFF(99, 30) ySHEC(30) ,FOC (30) »ROTC(30) »DMCB(30), CMCT (30) 
COMMON IUNB 
WRITE(6,371) 
371 FORMAT(50X,28HRESULTS OF FRAME ANALYSIS///) 
WRITE(6,372) MMMM,NCYCLE 
372 FORMAT(LOX, "CYCLE NO» MMMM=",13,2Xy"NCYCLE=",13///) 
WRITE(6,200) 
2CO FORMAT(10X,7HCOLUMNS//) 
WRITE(6,373) 
373. FORMAT(10X,10HFLOOR NO./,2Xy,1L9HJCINT ROTATION(RAD) »3X,L8HBCTTCM MO 
LMENT(KIN),3Xy16HTOP MOMENT (KIN) » 3X» LOHSHEARCKIP) » 3X, LOHFORCE(KIP) 
2) 
WRITE(6,374) 
374 FORMAT(10X,9HSTORY NO.//) 
WRITE(65375)ROTOC(ITER) 
375 FORMAT(12X,4HBASE ,10X,F13.6/) 
DO 377 K =1,NS 
WRITE(6,378)Ky ROTFF LITER 4K) yDMCB(K) yDMCT{(K)y SHEC(K) y FOC(K) 
378 FORMAT(13X_13y10XyE13 0698X9 F140 29 OX pF 14029 4X9 FL 0029 3X9F 10. 2/) 
377 CONTINUE 
814 WRITE(6,386) 
386 FORMAT(LHK) 
WRITE(6,201) 
201 FORMAT(//10X,5HBEAMS//) 
WRITE(6, 380) 
380 FORMAT(10X,9HFLOOR NOw,»3Xy16HMOMENT IN WALLy3Xe16HMOMENT IN WA 
LLL» 3X») 16HMOMENT AT COLUMN,3X~,1L6HMOMENT AT ROLLER,3Xy14HSHEAR AT 
2THEs3Xy14HSHEAR AT THE) 


WRITE(6,381) 

381 FORMAT(22X,16HSIDE BEAM AT THE,3X,16HSIDE BEAM AT THEy3Xy16HEND OF 
1 THE BEAM,3X,16HEND OF THE BEAMy3X,14HENDS OF WALL,3X,14HENDS 
2 OF BEAM) 

WRITE(6,382) 

382 FORMAT(22X,16HCOLUMN END(K-IN),3X,1L6HWALL END (K-IN),3Xy16HON RD 
ILLER(K-IN)» 3X, L6HON RCLLER(K-IN) »3X,14HSIDE BEAM(KIP),3X_,14HON RO 
2LLER(KIP)//) 

RETURN 


END 
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QUTPUT STATEMENTS FOR SHEAR WALL FORCES AND DEFORMATIONS 


SUBROUTINE SR3 {BMOMF,TMOMF,ROTByNS,F11L sFCC, SHEARyROTF,DEFF yCROTF, 
ICDEFF MMMM) 

DIMENSION F11(30),F0C(30), SHEAR(30) ,ROTF( 99,30) ,.DEFF(99,30) y 
1CROTF (99, 30),CDEFF(99,30) »,BMOMF (30) » TMOMF (30) 

COMMON TUNB 

WRITE(6,55) 

FORMAT (39X,52HSHEAR WALL ANALYSIS AND FINAL SLOPES AND DEFLECTIONS 
ry sy} 

WRITE(6,312) 

FORMAT(3Xy,lLOHFLOOR NO/s3X,15HWALL FORCE(KIP) » 3Xy LOGHFRAME FORCE(KI 
1P)53Xy,yL5HWALL SHEAR(KIP) 9X sl 7HWALL MOMENT(K-IN)»99Xy1LOHSLOPE(RAD) , 
23X%, 14HDEFLECTION(IN)) 

WRITE(6,401) 

FORMAT(3X,9HSTORY NO. y62X p6HBCTTOM, 1 3X» 3HTOP/S/) 

WRITE (6,20)ROTB 

FORMAT (5X s4HBASE,93X,E13.6/) 

N90 400 K=1,NS 

WRITE (6,313) KyFL1(K), FOC (K) » SHEAR(K) » BMOMF (K) y TMOMF(K),ROTF( 14K), 
IDEFF(1,«) 

FORMAT( 6Xo 139 7X 9F SB e2elL LX F802 IX gFlLOe2s SXgFl4e2 a 3Xe Fl ae 2 4XzsEl3eby 
13X,F13.6/) 

CONTINUE 

WRITE(6,404) 

FORMAT(1HK,50X y3LHCHECK ON SLOPES AND DEFLECTIONS//) 


CBAMKPCEREOHEL COR NO. ,5Xy LOHSLOPE (RAD) » 4X, 1L4HDEFLECTION(IN)//) 
NO 402 N=1,NS 
WRITE(6,403) N, CROTF (MMMM »N),CDEFF(MMMM,N) 
FORMAT (8X 9139 7X9E13 069 4X9 E130 6/) 
CONTINUE 
RETURN 
END 
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IN SUBROUTINE SR1 FORMATICN OF PLASTIC HINGE IN A MEMBER” CF 
THE FRAME [IS DETECTED 


CVO OVE Oye? 


SUBROUTINE SR1 (DMF PE, OMBF y0MB Ws DMWB yDMCT»DMCByMOMFB yMOMBF , MOMBW y MO 
DIMENSION DMFB (30), CMBF (30) ,DMBW{ 30) ,DMWB{ 30) »DMCT{ 30) »DMCB(30) 
REAL MOMFB (30), MOMBF (30), MOMBW( 30),MOMWB{ 30) ,MMCT(30), 
| 1MMCB(30),MPF(30),MPW(30),MPC(30) 

COMMON TUNB 

M3=0 

DQ 403 K=1,NS 

DMFB(K)=MOMFB(K) 

DMBF (K)=MCMBF(K) 

DMBW(K)=MOMBW(K) 

DMWB{K)=MOMWB(K ) 

DMCT(K)=MMCT(K) 

DMCB(K)=MMCB(K) 
403 CONTINUE 


DO 619 K=1,NS 
WRITE(6 7633) 
633 FORMAT(1HK) 
IF ((ABS(DMCB{K))).LT.MPC(K))GO TO 616 
M3=M3+1 
WRITE(6,617)Ks DMCBI{K) 
617 FORMAT(LOX,SOHHINGE AT BOTTOM POINT OF COLUMN IN STORY NOwyI 
13,5X,4O0HMOMENT AT BOTTOM POINT OF COLUMN =,F14.2,5H K-IN/) 
LE{DMCB(K) .sLT.0%02 GO TC 10 
DMCB{K)=MPC(K)+0.0001 
: GO T9 616 
: 1C DMCB(K)=-MPC{K)-0.0001 
616 IF((ABS(DMCT(K))).LT.MPC{KIIGO TO 602 
M3=M341 Pe 
8)K,DMCT(K 
618 BF dar iGn, SONHANGE AT TOP POINT OF COLUMN IN STORY NOwyI 
13,5X,4O0HMCMENT AT TOP POINT OF COLUMN =,F14.2,5H K-IN/) 
IF(DMCT(K) .LT. 0.0) GO TO 20 
DMCT(K)=MPC(K)+0.0001 
GO TO 602 
20 NMCT(K)=-MPC{K)-0.0001 
602 IF({ABS{DMFB(K))).LT.MPF(KIIGO TO 604 
M3=M3+] “a 
F MFR 
603 Pies criiay Sonnince AT ROLLER END OF BEAM CN ROLLER OF FLOOR NOwyI 
1345X,40HMOMENT AT ROLLER END OF BEAM CN ROLLER =yF14.2,5H K-IN/) 
IF(DMFB(K) «LT. 0.0) GO TC 30 
DOMFB(K)=MPF(K)+0.0001 
GO TO 604 ee 
30 DMFB(K)=-MPF(K)-0.- 
604 IF((ABS(DMBF(K))).LT»MPF(K)IGC TO 606 
M3=M341 
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695 FORMAT(10X,50HHINGE AT COLUMN END OF BEAM CN ROLLER OF FLOOR NO«yI 
13,5X,40HMOMENT AT COLUMN END SF BEAM ON ROLLER =5F14.2,5H K-IN/) 
IF(DMBFI(K) .LT. 0.0) GO TC 40 
DMBF(K)=MPF(K)+0.0001 
GO TO 696 
40 OMBFE(K)=-MPF(K)-0.0001 
606 IF((ABS(DMBW(K))).LT.MPW(KI)GO TO 608 
M3=M34] 
WRITE (6,607)K,DMBW(K) 
| 607 FORMATL1CX,50HHINGE AT COLUMN END OF WALL SIDE BEAM OF FLOOR NOey! 
1395X4OHMCMENT AT CCLUMN END OF WALL SIDE BEAM =,F14.2,5H K-IN/) 
IF(DMBW(K) .LT. 0.0) GO TO 50 
| DMBW{K)=MPW(K)+0.00C1 
GO TO 608 
50 NMBW(K)=-MPW(K)-0.0001 
608 IF((A3S(DMWB(K))).LT.MPW(K))GO TO 610 
M3=M3+1 
WRITE (6,609)K,DMWBIK) 
609 FORMAT(1CX,50HHINGE AT WALL END OF WALL SIDE BEAM OF FLOOR NOoyI 
13,5X%,4OHMOMENT AT WALL END OF WALL SIDE BEAM = yF14-2,5H K-IN/) 
IF(DMWB(K) .LT. 0.0) GO TO 60 
DOMWB(K)=MPW(K)+0.0001 
| 60 TO 61C 
60 DMWB(K)=-MPW(K)-0.0001 
610 CONTINUE 
RETURN 
END 
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APPENDIX C 


DERIVATION OF MOMENT-THRUST-CURVATURE RELATIONSHIP 


FOR REINFORCED CONCRETE CROSS-SECTION 
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C.1 Cross-Section 


The cross-section considered is shown in FIGURE C.1. The 
nomenclature used to define the cross-section is also shown. The 
subscript n refers to the total number of layers of steel. The steel 
layers are numbered starting from the tension face. Only rectangular 
cross-sections are considered. The stress-strain curves of the 


concrete and steel are defined in section 4.2.1 and 4.2.2, respectively. 


The corresponding areas of steel, ae and percentages, 


> 


bare sn : 
Ben a ahi of steel in each layer are given by Ay: Aco» Rogers 


ee and Pry Pros P 3 Sy eee R ae respectively. The total area and 


percentage of steel are given by equations (C.1) and (C.2). 


st 7 S] 52 53 


The strain in various layers of steel, starting from tension 
face ere qrven* by e(1), €(2)..is.... , e(n) whereas the strains at the 
extreme fibre of the section are given by Ey and ey at compressive 
and tensile (or least compressed) face respectively. The positive 
value of Ey shows compression on the tensile face whereas negative 
value shows tension on the tensile face. FIGURE (C.1) also shows a 


typical strain configuration. For any strain configuration the following 


hold true. 
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given by: 
and f is 
Piety 
by: 
or 


C3 


The curvature, o, is given by: 


ee? 
o= 4 ] 
t 
or ot = Spa Cs 3) 
The strain, «, at any depth, h, from compression face is 
e = Eq-oh (C.4) 
The strain in any layer of steel is given by: 
Aig loye este oo (C.5) 
4 t ; 
The stress, Lng in any layer of steel is given by: 
Tes S(O) Xx Es i) 1Ce6) 


Where E. is the modulus of elasticity of reinforcing steel 


the yield strength. 


The total axial force in the reinforcing steel will be given 
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The moment contribution of these steel layers will be 


given by: 
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The derivations of loads and moments relating curvature and 
fibre strains for different possible strain domains are described 
below. The reduction in the gross area of concrete due to presence of 


steel has been neglected. 


C.1.1 Sections under Compression Strain Only 


This strain configuration is shown in FIGURE (C.1). 
Considering an element of depth dh, at a depth h from top, 


where the concrete stress is Pes the axial load and moment will 


be given by: 
Pp 1 ae 
pT = A + Bt bf dh (C.9) 
and + = B + s J ; bf (-h) dh cs Gene) 
bt ie bt i 0 


Substituting the value of f, from equation (4.1) and integrating 


we have: 
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C.1.2 Pure Axial Load Capacity 


In the case of uniform strain over the cross-section the 
curvature becomes zero and equations (C.11) and (C.12) do not 
hold true, since M will be zero in this case. Therefore the 
value of P in this case can be found by integrating the uniform 


stress over the section which will be given by: 


2€ € 1.67 
wr = Att - 1.03 (=) Pee) 
c me 0 
M = 0 Ree nS 
and od = O os GS ep 


For the stress-strain curve presented in section 4.2.1, 
this will lead to a maximum value of P when Eq = 1.252, , eas 
the pure axial load capacity of, the section allowing for the 


reduction in gross area of concrete due to presence of steel, 


will be given by: 
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btt = A + 1.002 (1-p,) rey lpi Ks.) 


These values are essentially similar to those presented 


in reference (Al). 


C.1.3 Sections Partly in Tension but Uncracked 


The strain configuration is shown in FIGURE (C.2). 
The depth of neutral axis DN, is given by: 
DN = Eq /o 


The axial load and moment for this strain configuration 


will be given by: 


P 1 ‘DN t 
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btf. btt. ; c DN t 
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Diets ber 0 DN 
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Substituting the values of - and f, from equations (4.1) 


and (4.2) and integrating, we have 
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| C.1.4 Section Partly In Tension and Cracked 


The strain configuration is shown in FIGURE (C.3). 


The depth of neutral axis, DN, and the depth below which 


the section is cracked, dcr, are given by: 


and dcr = 


The axial load and moment for this strain configuration 


will be given by: 


Be fie er { ie bf cdh + Set be Anh, Wo ce 
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ber peer i we 
+ ge bron} .. 2. (6.22) 
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Substituting the values of i and f, from equations 


(4.1) and (4.2) we have 
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COMPRESSION <4 
FACE . 


TENSION 
FACE 
FIGURE C.1 
aa 
ies 
€, ae 
FIGURE - C..2 FIGURE C.3 


UazZ Program Nomenclature 


AS(N) 


FC] 
FC2 


FT2 


AXIAL LOAD CONTRIBUTION DUE TO REINFORCEMENT IN THE FORM 
OF P/bt IN bs. AND INCH UNIT 

GIVEN AXIAL LOAD ON SECTION IN 1bs. 

AREA OF REINFORCEMENT IN Nth LAYER IN IN¢ 

WIDTH OF SECTION (IN MAIN PROGRAM) IN INCHES 

MOMENT CONTRIBUTION DUE TO REINFORCEMENT IN THE FORM OF 
(M52), IN SUBROUTINE PM. 

DISTANCE OF CENTROID OF STEEL REINFORCEMENT IN Nth LAYER 
FROM BOTTOM FACE IN INCHES. 

THE RATIO C(N)/T 

STRAIN AT MID-DEPTH OF SECTION 

MODULUS OF ELASTICITY OF CONCRETE IN 1bs/IN@ 
ECCENTRICITY OF AXIAL LOAD IN INCHES 

c, (AS DEFINED IN EQUATION 4.1) 

MODULUS OF ELASTICITY OF STEEL IN 1bs/IN¢ 

LIMIT ON COMPRESSIVE FIBRE STRAIN, e, 

CRACKING STRAIN, €, 7, 

STRAIN AT BOTTOM FACE OF SECTION 

STRAIN AT TOP FACE OF SECTION 

CONCRETE STRENGTH IN DIRECT COMPRESSION, f1, in psi 
CONCRETE STRENGTH IN FLEXURE, 0.85f!, in psi 

STRESS IN Nth LAYER OF STEEL IN PSI 

YIELD STRENGTH OF REINFORCEMENT IN PSI 

CONCRETE STRENGTH IN TENSION IN PSI 


MAOF SHT MT: THANE oT aud norTua ATH | 
rity How OHA 2df salient 

L201 HT MOTOR WO BAD. AAT, WavTA 

Syn wr ava Aan Mt THEM DAOIIIA 70 ‘AaSA 

23H3M1 M1 (MARQORG MIAM wt) wortoae 90 ‘HTOIW 

* * 

#0. MAO? 3HT HI TH3MGOROANT 3A oT aud norTussrTaoD TMaMOM 
| ns SMC TUORRUE, HI A Sgacy: 


AaVAJ At VI THAMAIAOAMTRA 4133T2 30 GTORTHAD 30: JOMAT2IO 


2349K ME SOAR MOTTOS MOST. 

| T\(A)S OTTAA 3HT 

WOLTII2 TO NTGIG-OIM TA WEART 

Sy1\adf AI aT290N09 30 YTISITAAII 30 Su JUOM 
23HM WI GAOd whe de | 


; yt 30 30890 TA 

. tag nt.) <o1z2399"09 Toaan@ wr HTaAaATe 3TARONOD 
 Pequaty «5388.0 SAR ATRIA T3AQMOD 

129 WI one 5 at a) 


Cll 


M MOMENT IN THE FORM OF M/bt?F" 
NL NUMBER OF REINFORCING LAYERS 
NP NUMBER OF PROBLEMS TO BE ANALYZED 
| P AXIAL LOAD IN THE FORM OF P/bt ft. 
| PG GIVEN AXIAL LOAD IN THE FORM OF P/btt. 
PHIT CURVATURE TIMES THE DEPTH OF SECTION 
| PO ULTIMATE AXIAL LOAD IN THE FORM OF P/btf- 
PLL TOTAL PERCENTAGE OF STEEL IN THE FORM OF A. ,/bt 
| PT (N) PERCENTAGE OF STEEL IN Nth LAYER OF REINFORCEMENT IN 


THE FORM OF A, ,/bt 
| it TOTAL DEPTH OF SECTION IN INCHES 
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C.3 FLOW DIAGRAM OF THE COMPUTER PROGRAM 
C.3.1 MAIN PROGRAM 


READ NO. OF PROBLEMS 
TO BE SOLVED 


cs SET INDEX FOR THE 
PROBLEM 


READ AND WRITE DATA 
FOR THE PROBLEM IN 
PROCESS 


COMPUTE ULTIMATE 
AXIAL LOAD, Pa 


ES 
GIVEN AXIAL LOAD 
YES 


> 
ULTIMATE AXIAL LOAD 


(ar) COMPUTE e€ AND 
1 
CALL SUBROUTINE 'PM' 
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Cis 


IS GIVEN AXIAL LOAD} 
Pe. > 0:05 


IS 
ABS(PG-P) 
mire 


eS 


INCREASE 
ot BY 0.001 
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C.3.2 SUBROUTINE 'PM' 


COMPUTE STRESS AND STRAIN IN ALL 
STEEL LAYERS. SET THIS STRESS 
EQUAL TO f. INTHE LAYER OF STEEL 
WHICH HAS ¥ YI®LDED 


COMPUTE A AND B 


YES 
IS ot = 0.0 


NO 


COMPUTE P AND M 
FOR COMPRESSION 
FAILURE 


COMPUTE P AND M FOR 
UNCRACKED TENSILE 
FAILURE 


“COMPUTE P AND M FOR 
CRACKED TENSILE FAILURE 


COMPUTE P FOR M AND ¢ COMPUTE ECCENTRICITY 
EQUAL TO ZERO 
RETURN 
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C.4 LISTING OF THE PROGRAM 


MAIN C15 


DIMENSTON DT(10)ePT(10),FS(10),E110),AS{(1LO) ,C(10) 
REAL M 

READ(5,105)NP 

FORMAT(1X,13) 

NC 20 JJ3J=1,NP 

WRITE(6,510)JJ 

FORMAT(LHK/54X_ "PROBLEM NC.*,I12///) 
READ(5y,1LOO)ESyFCLyFY NL y Bey T AL eg EUL 
FORMAT(1XsE9.1 sF7 a1 oF 8.191292 F602 F901, F906) 
REAC(5,101){AS(K) ,K=1LyNL) 
FORMAT( (1X ,10F7.4)) 

Renu os 10294C( 0) L=1,NL) 

FORMAT ((1X,10F7.3)) 

Woe 5" T=15'NL 

PT(1)=(AS(1)7(B*T)) 

NTCVI=HCCC1I/T) 

CONTINUE 

FC2=0.85*FC1 

EC=5740C.0*{(FC1 **0.5) 
EQ=((2.9*FC2)/EC) 

EY=FY/ES 

FI2=7 .0*(FCL*¥*0.5) 

EU=0.09040 

WRITE(6,516) 

FORMAT(50X,"DATA FOR THE ANALYSIS!///) 


WRITE(6,517) 
FORMAT{5X,"E FOR STEEL ',3Xy*CONC. STRENGTH',3Xy"STEEL STRENGTH!® ,3X 


L,*NO. OF STEEL LAYERS" ,3X,"GIVEN AXIAL LOAD",3X,*CRACKING STRAIN*/ 
2) 


WRITE(6,519)ES,FC1l,FYyNL,AL,EUL 

FORMAT( 6X9E9.357TXoF Tol pIX oF Be ly lSOXeg L291 5X eF9.1,10XeF9L6///) 
PTL=0.0 

WRITE(6,103) 

FORMAT (5X,°% OF STEEL’ »3X,*COVER RATIO'/) 

Deen.) = 1» NL 

WRITE(6,104)PT(J),0T{J) 

FORMAT( 7X9F 7049 7XyF623/) 

PV ES=PTL+PT( 3) 

CONTINUE 

FIND ULTIMATE AXIAL LOAD 
PO=(1.002%(1.0-PTL))4#( (PTL*FY )/FC2) 

CALCULATICNS OF MOMENTS FOR GIVEN AXIAL LOAD AND CURVATURE 


WRITE( 6,514) 
FORMAT (34X, "RESULTS FOR THE MOMENT CURVATURE RELATIONSHIP FOR GIVE 


INFAXTAL LOAD*/ 7) 


WRITE(6,513) 
FORMAT(5X,"TOP STRAIN's3Xy"BOTTOM STRAIN® y3Xy "AVERAGE STRAIN", 5Xy 


LCURVATURE" 98Xy "ECC. 'y10X, AXIAL LOAD* 8X9 "MOMENT! y3Xe°GIVEN AXIAL 
2LOAD"/) 


PG={AL/(B*T*#FC2)) 

IF(PG.GE.PC)GO TO 20 

E4=0.090 

PHI T=0.0 

E1=E4—-PHIT 

ae PM{E4, Ely PHITsPsMy FEO ,EULy FC2yFSyPT 2 ESaFY¥eNT,E,FT2,ECC gNL) 
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GO TO 47 
48 L=ABS(PG-P)/PG 
IF(Z.LE.0.010)G0 TO 36 
47 IF(E4.GT.0.0039995)G0 TO 36 
39 E4=F44+0.0000010 
SHI TU 37 
36 EAV=(E4+4+E1)/2.0 
WRITE(6,520 )E4,EL,EAVy PHI Ty ECC yPyMy PG 
a20 FORMAT (6X 9F 8.6 96X 9F 8069 9X 9 FB 6 69 BX oF Ba Ge TX 9 FB 04s BX 9 F904 9 OX pF G0 Ge Xe 
1F9.4/) 
IF(E4 .GT.0.0039995)G0 TO 20 
IF(PHIT.GT.0.00295)GO0 TO 45 
PHIT=PHIT+0.0003 


GO TQ 39 

45 PHIT=PHIT+0.001 
GO TO 39 

20 CONTINUE 
STOP 


END 
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‘PM! 
meee TNE PMCE4,E1L,PHIT, PyM, EQ,EUL s FC2,FSy PT sESyFY ys DT Ey FT2yECCyN 
DIMENSION FS(10),0T(10),PT(10) ,E(10) 

REAL ™ 

A=0.0 

B=0.0 

DO 53 J=1,NL 

E(J)=E4-(PHIT*#(1.0-CT(5))) 

FS(J)=E(J)*ES 

IF(E(J).LT.0.0)G0 TO 54 

IF(FS(J).GT.FY)FS(J)=FY 

GQ 10 55 

IFC CABSC(FS(J)) ) GTsFYIFS(J)=-FY 

A=A+(PT(JS)*FS(J)) 

B=B—( PT (J) *(1.0-(2.0*DT(JS)))* FSC) 

CONTINUF 

IF(PHIT.EQ.0.0)G0 TO 25 

Peiri.1 3 .0.0)G0 TO 20 

COMPRESSION FAILURE 

P=(A/(1.0*FC2) )+((E44E1 )/E0)-( (0.386% ( (F4**2.67)-(E1¥*2.67) IS (EO 
1**1.67) *PHIT)) 

M=(B/(2.0*FC2))+(PHIT/(6.C*EO)) 
1-( (0.193% ( (E4**2.67)4( EL®*2.67))9/¢ (EO**1L £67) XPHITI) +#( (0. 105*( (E4% 
2*3267)- (E1**3.67))) SCL EDEL 267)*{ PHIT#*220))) 

ee eee | 

IF(EL.LT.EULIGO TO 26 

P=(A/(1 .O*FC2) 1 4( (E4**2 00 )/( EQ*PHIT))-( (0. 386*( E4**2. 

167))/¢ ( EO#*1.67)*PHIT)) #C (FT 2/FC2 08 (( (CE L**2) / (EUL*PHIT) )-( EL**®4 
2)/(4.0*(EUL **3 ) *PHIT))) 

M=(B/(2.0*FC2))4((E4%*2.0)/(2.0*EO*PH 

LIT) )-€( £4**3.0)/( 3 -OXEO*( PHIT#®2.0) )I-€10. 193% (E4**2.67))/((EO**1. 
267) ¥PHIT) )4+( (0. 105% (E4**3 067) )/(( EO#*L 667) *{ PHI T#*2 20))9-( (FT 2/FC2 
3)y*(( (E1L**29/(2.O0XPHIT *EUL )) + ((E1**39/(3.0*( PHIT#*2.0) XEUL) -(( 
4ELER4)/ (8. O*PH I TX (EUL ®*3)))-( (E1**5)/(20. 0¥ (PHI T#*2 0) *( EUL* 


5*3)))7) 
GU fC "21 


CRACKED TENSILE FAILURE 
P=(A/(1.0*FC2) )4+((E4**2.0)/( EQ*PHIT )I-( (0.386 *(E4**2. 


167))/( (EO**1 067 )*PHIT)) #( (FT 2*0. 75*EUL)/( FC2*PHIT)) 
M=(B/{2.0¥*FC2) 4+ (E4**2.0)/(2.0%E0*PH 

LIT) )—( CE4E*3 00) /(3 O¥EO#( PHIT#®2 0) ))-(106193%(F4**2067))/( (EO**1. 

267)*PHIT) 4 ( (00 LO5*(E4**3 267) )/( (EO##L 667) *( PHIT##2 ODD HK CFT 2/FC2 
3)*( (0 (3 .O*EUL) / (4 .O#PHIT) )*(0. 5-(E4/PHIT) #(EUL/PHIT)) )-(17 .0*( EUL* 


4%*2)/(60.0*{PHIT**2.0))))) 

G70 -2.1 

P=(A/FC2)4+(2.0 *E4/EQ)-(1.03¥( (E4/EO)**1.67)) 
M=0.0 

ECC=0.0 

IF(P.FEQ.9.0)G0 TO 22 

ECC=M/P 

RETURN 

END 
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D.1 Curvature and Fibre Strains from Dial Gauge Readings for Columns 


and Beams 


FIGURE D.1 shows two curvature meter arms mounted on a section 
of depth (h+g). The dial gauges are located at distances b and c from 
the centre line of the section as shown. Let E, and €, are the extreme 
fibre strains as shown. M and (M+1) are the dial gauge number on the 
station number, N which is numbered (M+1)/2 where M is an odd number. 
DR(M+1,K) and DR(M,K) indicates the dial readings in Mth and (M+1)th 
gauge for load increment number K. The dotted line shows the bent 
position of the section. Therefore, if the strain profile over the 
cross-section is extended then the average strain over the length, a, 


at a distance b and c, as shown, from the centre line will be given by: 


G,/a {DR(M+1,K) - DR(M+1,1) } /a Rr 48) 


{DR(M,K) - DR(M,1)} /a Mere, (ieee: 


G,/a 


Therefore the curvature will be given by: 


(G, - G) 
? 4 1 (0.3) 
0) 3 —a(bte) hs 
and from similar triangles the fibre strains will be given by: 
EA = — g(btg) + Gy/a set ie 
Ey = o(b-h) + G,/a eee 


The above derived equations hold true for all possible strain 


configuration. 
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D.2 Curvature and Fibre Strains for Wall 


FIGURE D.2 shows the placing of Whittemore or Demec points at 
a particular section. The reading on the gauges will give an average 


strain over the length, a, if divided by the gauge length. Let Y12¥5 
1? Xo and Xo from the 


inner’ face of the wall. The strains Yq: V5 and Y3 are given by the 


and Y3 represent the strains at a distance X 


difference in reading at some load increment K and the initial readings 
divided by the gauge length for the corresponding gauge line. At each 
load increment a best fit straight line was fitted to the strain 
readings using least squares. This was done assuming that a straight 
line distribution of strains existed at every section and the equation 
of straight line is given by: 


ms a, x fe OL) 


The inner face was taken as the origin. The constant a, and ay in the 


above equation can be derived (S1) as follows: 


IX = x + % 4 ; ee eye) 


= = X ee ee D.8 

7 peace ls Yo t 5 an + a, 2 (D.8) 

eee a 2 2 - 

SRP Paks Ae (D.9) 
= 2 


where n is the number of known points, in this case, three. Solving 


equations (D.8) and (D.10) gives: 
. rY.=uX - nek -..(D.11) 
| (xx)é - n©x 
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The corrected fibre strains can be computed using equation 
(D.6) by substituting the value of a, and a, and proper value of X. 
To find Ey the value of X will be zero whereas E, will be given by 


Substituting t for X. 
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D.3 Program Nomenclature 


A(J) GAUGE LENGTH OF Jth STATION, OVER WHICH THE MEASUREMENT 
HAS BEEN TAKEN, IN INCHES. 

AA A TEMPORARY LOCATION USED IN SUBROUTINE 'BMKP'. 

AS(N,L) AREA OF REINFORCING STEEL IN THE SECTION AT Nth 
STATION AND Lth LAYER 


AO A CONSTANT IN THE LEAST SQUARE LINE EQUATION 
Al A CONSTANT IN THE LEAST SQUARE LINE EQUATION 
B(J) DISTANCE OF DIAL GAUGE, IN CURVATURE METER, FROM THE 


CENTRE LINE OF THE SECTION, AT Jth STATION, TOWARDS 
OUTER FACE OF COLUMNS OR BOTTOM FACE OF BEAMS, IN 


INCHES. 
BB A TEMPORARY LOCATION USED IN SUBROUTINE ‘BMKP' 
BM(N,K) MOMENT IN KIP-INCH AT Nth STATION FOR Kth LOAD 
BMG(1,N) MOMENT, IN THE GIVEN M-@ RELATIONSHIP, CORRESPONDING 


TO Nth POINT IN Ith CURVE, IN THE FORM OF M/bt ert : 
C(J) DISTANCE OF DIAL GAUGE, IN CURVATURE METER, FROM THE 

CENTRE LINE OF THE SECTION, AT Jth STATION, TOWARDS 

INNER FACE OF COLUMNS OR TOP FACE OF BEAMS, IN INCHES 


CA TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 
CB TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 
CC TEMPORARY LOCATION USED IN SUBROUTINE 'BMKP' 
CCA TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 
CCB TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 
CO(N,L) DISTANCE OF CENTROID OF STEEL IN THE SECTION, In Lth 


LAYER AND Nth STATION, FROM THE TENSION FACE, IN INCHES 
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CURV(I,N) 


G4 


Gl 


LS(N) 
NC 


CURVATURE, IN GIVEN M- RELATIONSHIP, CORRESPONDING TO 
Nth POINT IN Ith CURVE, IN THE FORM OF ot. 

TEMPORARY LOCATION USED IN MAIN PROGRAM 

TEMPORARY LOCATION USED IN MAIN PROGRAM 

TEMPORARY LOCATION USED IN MAIN PROGRAM 

READING OF Ith DIAL GAUGE FOR Jth LOAD 

DISTANCE OF Jth STATION FROM THE BOTTOM OR LEFT END OF 
THE MEMBER, ON WHICH THE Jth STATION LIES, IN INCHES 
MODULUS OF ELASTICITY OF STEEL IN Jth MEMBER IN PSI 
STRAIN AT THE OUTER FACE OF COLUMN OR BOTTOM FACE OF 


BEAM OR INNER FACE OF WALL, AT Nth STATION FOR Kth LOAD 


STRAIN AT THE INNER FACE OF COLUMN OR TOP FACE OF BEAM 
OR OUTER FACE OF WALL AT Nth STATION FOR Kth LOAD 
CONCRETE STRENGTH OF Jth MEMBER IN PSI 

YIELD STRENGTH OF STEEL IN Jth MEMBER IN PSI 

DISTANCE OF COMPRESSION FACE OF THE SECTION, AT Jth 
STATION, FROM CENTRE LINE IN INCHES 

DIFFERENCE IN DIAL GAUGE READING AT ANY LOAD, FROM THE 
INITIAL READING ON THE INNER FACE OF COLUMN, OR TOP 
FACE OF BEAMS OR OUTER FACE OF WALL 

DIFFERENCE IN DIAL GAUGE READING AT ANY LOAD, FROM THE 
INITIAL READING ON THE OUTER FACE OF COLUMN OR BOTTOM 
FACE OF BEAMS OR INNER FACE OF WALL 

DISTANCE OF TENSION FACE OF THE SECTION, AT Jth STATION, 
FROM CENTRE LINE IN INCHES 

NUMBER OF LAYER OF STEEL IN THE SECTION AT Nth STATION 


TOTAL NUMBER OF CURVATURE METER 
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ND 
NGC 
NL 
NM 
NPC (N) 


NS 

NSM(J) 
NWD 

NWL 

PHI 
PHIT(N,K) 


TOTAL NUMBER OF 
TOTAL NUMBER OF 
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CURVATURE CURVE 
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DIAL GAUGES IN CURVATURE METERS 
SUPPLIED MOMENT-CURVATURE CURVE 
LOAD INCREMENT 

MEMBERS IN THE TEST FRAME 

GIVEN POINTS TO DEFINE Nth MOMENT- 


MEASURING STATIONS 
MEASURING STATIONS IN Jth MEMBER 
WHITTEMORE AND DEMEC GAUGES 


NUMBER OF WHITTEMORE OR DEMEC LINES ON WALL 


CURVATURE. AT ANY STATION 


CURVATURE TIMES 
Kth LOAD 


DEPTH OF SECTION AT Nth STATION FOR 


TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 


TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 


TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 


TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 


TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 


TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 


WIDTH OF TEST FRAME IN INCHES 


DISTANCE OF Jth 
FROM INNER FACE 


LINE FOR WHITTEMORE OR DEMEC GAUGE 


OF WALL 


TEMPORARY LOCATION USED IN SUBROUTINE 'FSAC' 
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D.4 FLOW DIAGRAM OF THE COMPUTER PROGRAM 
| D.4.1 MAIN PROGRAM 


READ ALL THE DATA REQUIRED 

FOR THE ANALYSIS 

WRITE OUT ALL THE DATA 
SUPPLIED 


CORRECT ALL THE INITIAL 
GAUGE READINGS AND WRITE 
OUT 


CALL 'BMKP ' 
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D.4.2 SUBROUTINE 'FSAC' 
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COMPUTE FIBRE STRAINS AND 
CURVATURE IN BEAM AND COLUMN 
SECTIONS FOR ALL LOADS 


COMPUTE FIBRE STRAINS AND 
CURVATURES IN WALL SECTIONS 
FOR ALL LOADS 


WRITE OUT FIBRE STRAINS AND 
CURVATURES AT ALL SECTIONS 
FOR ALL LOADS 


RETURN 
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D.4.3 SUBROUTINE 'BMKP'! 


SET INDEX FOR STATIONS 
AND LOADS 


SELECT PROPER M - ¢ CURVE 
TO COMPUTE MOMENTS 


WRITE ‘CURVATURE 
BEYOND THE M - CURVE 
GIVING STATION AND 

LOAD NUMBER 


IS CURVATURE 
OUT OF RANGE ? 


INTERPOLATE MOMENTS 


ARE ALL MOMENTS 
INTERPOLATED ? 


WRITE MOMENTS AT ALL 
STATIONS FOR ALL LOADS 
RETURN 
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D.5 LISTING OF THE PROGRAM 
Sa ee aay ee ae oS ee D12 


MAIN 
REBDYANDYARITE ALL. POSSIBLE  BDATAS OF THE TEST FRAME 
DOUBLE PRECISION DRyWR,E4,E1,PHIT »C1,C2,C3 
DIMENSION A150) ,8(25),0(25)9G6(50),H(50),0R(50915),E4( 50515) ,T(5C), 
1E1(50515),PHIT (50,15) »AS(50 910) -C0(50,10) »LS150) eWRI75,15) »XT10),N 
2PC(35) »CURV(35 4 70),NSM(1L5 )ySM(15) »0S(50),BML(15915) ,BMR(15 415) 9 BMG 
Bis 5401s OeM(S0, 15) ,FCL(15) ,FY(15),E€S8 (15) 
RFACD(5,1L01LIND NWO, NWL a NL NM aNCCoW 
101 FORMAT(1X,614,F6.2) 
NC =ND/2 
NS=NC#NWD/NWL 
READ( 5,503) (NSM(J)_,J=1_9NM) 
503 FORMAT((1X,2613)) 
READ( 5,799) (SM{J) 9 J=19NM) 
799 FORMAT((1X,13F6.2)) 
READ(5,350)(FC1(J),J=1l,NM) 
850 FORMAT((1X,1l1LF7.1)) 
READ(5, 851) (FY(J) »J=1l_NM) 
851 FORMAT((1X,9F8.1)) 
READN(5,30CI(ES( J) ,J=1L9NM) 
3090 FORMAT((1X,3E9.1)) 
READ(5,503)(LS(J)eJ=19NS) 
NO 30 N=1,NS 
T=LSUN) 
READ(S, 102) ((AS(NoL) yCO(NL)) pL=1l5T1) 
390 CONTINUE 
READ(5,102)(A(J)_J=1,NS) 
REAC(5,102) (B(J)5 J=1sNC) 
READ(5,102)(C{ J) 9 J=lyNC) 
READ(5,102)(G(J),J=1,NS) 
READ(5,102)(H(J)_ J=lsNS) 
READ(5,102)(D0S( J) J=1leNS) 
READ(5,102)(X( J), J=lsNWL) 
102 FORMAT((1X,9F8~.4)) 
READ(5, 107) (( ORCL J) 9 d=LyNL)y T=1yND) 
READ(5,1O7)C(WRIT sd) eJ=LoNL)y L=1lyNWD) 
107 FORMAT((1X,9F8.5)) 
READ(5,503) (NPC(N) yN=lyNGC) 
N0 600 [=1,NGC 
JJ=NPC(T) 
READ(5, 502) ({CURV(I,N) sN=leJJ) 
REALC(5,502)(BMG(IyN)N=1s JJ) 
502 FORMAT( (1X,11F7-4)) 
600 CONTINUE 
WRITE(6,359) 
350 FORMAT(50X,'VALUES CF CCNSTANTS®//) 
WRITE(6700)NMyNL pNOaNWD a NC yNSaNWL sy NGC IW 
TOO FORMAT(1H0,5X_"NO- CF MEMBERS IN THE FRAME =",13/6X_,'NO. OF LOAD I 
INCREMENT =¥,13/6X,"NO- OF DIAL GUAGES =',13/6X%,"NO. OF WHI TTEMORE 
DIAL GUAGES =*,13/6X_,"NG. OF CURVATURE METERS =',513/6Xy'NOe OF POS 
3ITIONS ="',13/6X%y'NO- OF WHITTEMORE LINES =9,13/6X,'*NO. OF GIVEN M- 


6 852) 
WRITE(6y eMEMBER NOo' 93Xe"CONC. STRENGTH(PSI) '93Xy STEEL STR 


852 FORMAT(1HO, 10X 
7 ENGTH( PSI)" 3X 9 !SPANCIN) #9 3X9 'NO-OF STATIONS'y3Xy*E FOR STEEL'//) 


DO 853 N=1lyNM 
WRITE (6 4854 )Me ECL (MDs FY (MD 6 SM(M) sNSMCM) 2 ESM) 


854 Beat Ahy 13 2 UBh ob te Lele he PO el prOk hea Beeler hath Tats 


(ee ibe 


i ae 


‘ iy ay: in 4 ed Mi 
AIAN if. 


A 


MA9A T23T IMT AB Zara | sasnzece, jak aria 
| 3s 1SIeh De TYH% 9 [3483 MWe AO wot213399 
NT, HOP p OP HI, 171, OF FRO, (OOP, £02) De (ESI Dy (2S )Gy LOFEA. Wi 3 
LM et OL. ey) As (O#] ‘tee eee IDs {04,08 2H 412108) TINO 2 f Oe 


(Efe FEV AMA. (OTe FI) ING ¢ (OB VEO, CAL) Meee MeM OT. BEIVAUD, (28 
CEP 2F AREY, (REND LOLs ee 
rh Vg INV OWihs OVO Poe: 
J) Vee +81 ATO KET 
x \oy : 
7 su OW OA | 
(Mit. [= =LethI M2) (£088 
| ae | © “Cle Tas eriits 104 €¢ 
wr MA, bebe (LYM (OT 92 10A35 


me MOS. 0 FEL 028 (oe 
MMe Laue (LITIINO@B ge 
((f. TRL e x00 TT AMR 
iM, d= Le (LIYA) (fee ee 4 
(tl. BAC XLVI ae aC 
(MVie f= Ay (LEZIN(O0Es 
i “W4Ee 038k L1) ry 
(2m et=be(L) 24) (£0Re OASs 
enw t= € 
| ; - a 2 I= 
(Pe fade (Se K) 024 (Je) 2407 (S01 9290s 


ta _ URI TAG 
(2VelHbe (LIAI(SOL 2004 
(DW. i=l (LI0)( S01 42 )GASS 

(OKs I=L (LIDICSOL gE fhe 

row faUetL O(SOL 52 OA 

ith (2Ueleb. (LUHD(SOL, oe 
ai t2 Mel =be (L120) (50 48 ORs 
AA) I ASO des OA 
€4d, 890,XI)) c 
ton. anaead Pls (hy TIAOVINTOL Ge) OA 
(OHM gL=T eCity LOSE TIAWI ET 12109 
St | (2. 8 IO eXL) IT MA05 
‘Son ronctmSa toe es 
: “OC 

| 41) 994st) 
“Utes ainsi ae Xs 35 


ae 1 eee eee yanan 
(ih ye Se (OT ALLL) ; 
iy \ ~ i nh tay ve ‘ ae wie 
ns 2 THATeMDD. 90. 230M 
ae 


Ws DOW AW Ze Ie OWN ¢ Dy JM e MI 
RO. OM! eXO\EL = IMART ABT AL 2n3eMay 49 LOK" ’ 
ATH! 30 LOW XoVeT,y "= 230Au0 Secrets 
ete NONE *s 2827 3M ARUTAVAUD 30 | 


fi ja 40 OW", fly*= 239012 3AQMa’ he: 
2 a ‘ANAS: oR ete — se nto a eine 
‘gel OTe Uae 
are “aXe, C129) H Ta WIAT2 ng exe, 146 ' 838 nage i) fs aa st . 
| 2 909 i Re Pd CAT IMARE® exe *4129)HTO 
| mes rar 


: be ; th ee ed Ny j My tes ree 
be oeyzoan | eee en He 9M YAGEMY OFM AE y 
tf 1%, KO cE eR eS of : a . . ibn 


D13 
852 CONTINUE 
WRITE(6,5353) 
353 FORMAT(1H0,45X, "REINFORCEMENT AREAS AND COVER'//) 
WRITE(6,354) 
354 FORMAT(LX,*POSITICN NO. *) 
NO 31 N = 1,NS 
R=isS tN} 
WRITE(65355)Ny(LyL=1l,1) 
355 FORMAT(8X¥1394Xe"STEEL LAYER NO.%94X_13,9ITL0/ (34X—913,9110) ) 
WRITE(6,356)(AS(NoL),L=1,17) 
356 FORMAT(15X,"STEEL AREAS',4X,10F10.4/(30X%, 10F10.4) ) 
WRITE( 6,357) (CO(NyL I yL=1,1) 
357 FORMAT(15X, "COVER", 10X,y10F10.4/(30X »L0F10.4)) 
31 CONTINUE 
WRITE(6,704) 
704 FORMAT(1H0,30X,"*VALUES OF CONSTANTS FOR MEASURING STATIONS") 
WRITE(64150) 
150 FORMAT(1HO,10Xy"*POSITICN NOo%s7Xo "G%y LOX 9 "Hy IX y*DS8yLOXy "A*y LOX, ® 
1B*,10X,"C'//) 
D0, 50, J=1,NC 
WRITE(6 5151)39,G(J),H{J) sDS(J) pA( JS), BUS) C(I) 
151 FORMAT(13X_,13,6X,6F11.4) 
50 CONTINUE 
MM=NC 41 
DO 52 J=MM,NS 
WRITE(6 9145) Jy G(J),H(JS) pDS(S) pA J) 
145 FORMAT(13X_,13,6X,)4F11.4) 
52 CONTINUE 
JM=NWL-1 
WRITE (69116) (X(N) »N=1_ JM) 
116 FORMAT(1H0O,10X,"WHITTEMORE LINES ARE AT A DISTANCE OF ',(10(F6.2,y°, 
L*)/)) 
WRITE(6,117)X(NWL) 
117 FORMAT(11X,*AND',F6.2," INCH FROM THE OUTER FACE OF WALL*) 
WRITE(6,152) 
152 FORMAT(1H0O,50X,"DIAL GAUGE READINGS*//) 
WRITE (69153)(KyK=1,NL) 
153 FORMAT(6X_,"DIAL NCe/LCAD NOw* 9 3Xe1359T10/ ( 26X_13,9T10)) 
DO S51 N =1,ND 
WRITE(65154)Ny (DRINGK) »K=1—NL) 
154 FORMAT(9X,13310X,10F10.5/(22X,10F10.5)) 
51 CONTINUE 
WRITE(6,115) 
115 FORMAT {1H0»50Xy ‘WHI TTEMORE GAUGE READINGS'//) 
WRITE(6 9153 )(KaK=1 5NL) 
DO 53 N =1,NWD 
WRITE(6,154)Ny (WRIN 9K) »K=1,NL) 
53 CONTINUE 
DO 601 I=1,NGC 
WRITE(6,701)1 
701 FORMAT(1HO,20X%9'M-PHI CURVE NO.',13/) 
2 A ee eee NO." ,3X,"CURVATURE X T"y3Xy*MOMENT/BITFC2%/) 
JJ=NPC(T) 
0 602 N=1,JJ 
Mee ik 303) N; CUR VCH UO 
7O3 FORMAT(14X_1398X9FT7 042 8X9 F7 4) 
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692 CONTINUE 
601 CONTINUE 
99 75 N=1,ND 
C 3=DR(N,3)-DRIN,2) 
IF(0ABS(C3).EQ.0.0)GO TO 90 
02=44.44/C3 
C1=5.77-C2*NR(Ny2) 
DR(N,1)={-C1/C2) 
GO> fr 17's 
90 DRI{N,g1)=DRIN,g2) 
75 CONTINUE 
DO 76 N=lyNWD 
C 3=WR(Ny3)-WRUNY2) 
IF(DABS(C3).EQ.0.01GO T9 91 
C2=44 .44/C3 
C1=5.77-C2*WRIN 2) 
WR(N,L)=(-C1/C2) 
GO TO 76 
91 WR(N,1l)=WRIN,2) 
76 CONTINUE 
WRI TE(6,400) 
400 FORMAT(1HO,10X%,'CGRRECTED DIALS AND WHITTEMORE INITIAL READINGS®// 
1) 
DON7T7 "N=1,ND 
WRITE(6,40LIN,DRINy LI) WRONG L) 
401 FORMAT(10X,13,2F190.5) 
77 CONTINUE 
TI=ND+1 
DO 78 N=II,NWD 
WRITE(6,402)N,WR(Ny 1) 
402 FORMAT(10X,13,10XyF 10.5) 
78 CONTINUE 
CALL FSAC{NCyNDgNWDyNWLyNLyNS oA eB oC 9GyHs£4y Ely PHITs DReWRy X) 
CALE BMKP (PHIT »NSeNLoFCly We Ty NGCyNMaNPCyCURV? BMG aNSMy BM) 
STOP 
END 
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SUBROUTINE FSAC UNC s ND a NWN yNWL sNLyNS 9 A9B Cy Ge Hy 4y Ely PHI Ty DRyWR, X) 
PROGRAMME FOR CONVERSION CF OTALS AND WHITTEMCRE READINGS TO 
STRAINS AND CURVATURE 

DOUBLE PRECISION DReWRrE4eF1lyPHIT +G4yGl yPHIy Vy SY¥ySXYyAL,AO 
DIMENSTON A(50)5B(25)4C(2£) 46150) H(50) 4E4(50915),E1( 50915) »PHIT(S 
19915)5DR(50515) WR( 75415) sX(10),Y(10) 


an 
ie 


58 


81 


61 


82 


33 


DO 22 M=1,ND,2 
N=(M4+1)/2 

DO 23 K=2,NL 
S4=DR(M41,K)-DRI(M4+1,1) 
G1=DR{M,K)-DRI(M,1) 
PHI=(G4-G1L)/{A{N)*(B(N)}4C(N))) 
E4(NyK)=PHI*(BIN)4+G(N) )4#G1/A(N) 
EL(N,sKI=PHI*{B(N)—-H(N)DVFtGI1/ACN) 
PHIT(N,K)=E4(N,K)-EL(N,K) 
CONTINUE 

CONTINUE 

SX=0.0 

SX2=0.90 

DO 58 N=1l,NWL 

SX=SX+#X(N) 
SX2=SX2+( X(N) **2) 

CONTINUE 

KK=NWO/NWL 

CA=FLOATINWL) 
CB=(SX**2)-(CA*SX2) 

DO 59 J=1,KK 

N=NC+J 

DC 6C K=2,NL 

SY=0.0 

SXY=0.0 

TJ=0 

SCX=SX 

SC X2=SX2 

DO 61 £ =1,NWL 
M=J4+{((I-1)*KK) 
IF(WR(M,K).EQ.0-0)GC TC 80 
YC IT =CWRUM,KI-WROUMs 1D) SAN) 
GO TO 81 

Y{(1I)=0.0 

[J=1IJ+1 

SC X=SCX-X(1) 

SCX 2=SC X2—(UX(1 ) #*2) 
SY=SY+Y(1) 

SX Y=SXY#X(1I)*Y (1) 

CONTINUE 

IF(IJ.GT.0)GO TO 82 
Al=(SY*SX-CA*SXY)/CB 
AO=(SY-A1L¥*SX)/CA 

GG TD 83 

CCA=FLOATINWL-IJ ) 

CCB=( SCX**2)-{ CCA*SCX2) 
Al=(SY*SCX-CCA*SXY)/CCB 
AQ=(SY-AL*SCX)/CCA 
E4(N,K)=AO 
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EL(N,K)=A0+A1L*(G(N) #H(N)) 


PHITIN,K)=E4(N,K) —-EL(N,K) 
60 CONTINUE 
59 CONTINUE 
WRITE(6,200) 
200 FORMAT(1H0,25X,*I NNER FACE STRAIN IN COL. OR TOP STRAIN IN BEAMS 0 
1R OUTER FACE STRAIN IN WALL*//) 
WRITE(6,201)(K,K=2,5NL) 
201 FORMAT(2Xs*POSITION NO./LCAD NOW *, 3X_,13,9T10/(26X%_,13,9110)) 
NO 25 N=1,NS 
WRITE(69202)Ny (E4(N eK) »K=2,NL) 
202 FORMAT(9X,13,10Xy L0F10.6/ (22X y10F10.6)) 
25 CONTINUE 
WRITE(6,203) 
203 FORMAT(1HO,25X,*OUTER FACE STRAIN IN COL. GR BOTTOM STRAIN IN BEAM 
1S OR INNER FACE STRAIN IN WALL'//) 
WRITE (65201)(K,K=2,NL) 
DO 26 N=1,NS 
WRITE(69202)Ny (ELIN,K) »K=29NL) 
26 CONTINUE 
WRITE(6,206) 
206 FORMAT(1H0O,50X,"CURVATURES AT VARIOUS POSITIONS'//) 
WRITE(6,201)(K,K=2, NL) 
DO 27 N=1,NS 
WRITE(6,202)Ny ({PHITIN,K) »K=29NL) 
27 CONTINUE 
RETURN 
END 
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SUBROUTINE BMKP(PHITyNSyNL oFC 1 yWe Tp NGCyNMeNPC 
PROGRAMME FOR FINDING MOMENTS FROM KNCWN AEA a TNY a 
DOUBLE PRECISION PHIT 
DIMENSTCN NPC(35) sCURV(35 970) ¥BMG(35y 70) y PHIT (50415) ¢ BMC 50915) »NSM 
1(15),T(50),FC1(15) 
90 605 M =1,NM 
TA=NSM(M) 
IFAM.SEQ.1)G0 TO 626 
KK=NSM(M-1) +41 
GQ TC 627 
626 KK=1 
627 DO 603 N =KKyIA 
DO 604 K=2,NL 
IF(DABS(PHIT(N,K)).LT.10-9)GO TO 625 
IF(M.LE.4.0R.M.GE.9)GO TO 621 
I=M 
IF(INEQ.19)1=2 
IF(N.EQ.21)1=3 
IF(N2EQ.23)1=4 
GC TO 610 
621 IF(K.EQ.2)G0 TO 608 
IF{K.EQ.NL)GO TO 609 
IF(M.LE.4) 1=M+8 
IF(M.GE.9)IT=M+12 
GO TO 610 
608 IF(M.LE.4)1=M4+12 
IF(M.GE.9)I=M+16 
GO TO 610 
609 IF(M.LE.4)1=M+16 
IF(M.GE.9) L=M+20 
61C JJ=NPC(T) 
Domolle) =2,3J 
A=(DABS(PHIT(N,K)))-CURVUIIsJ) 
IF{A.GT.0.0)GO TO 612 
IF(A.EQ.0.0)GO0 TO 613 
IF(A.LT.0.0)GO TO 614 
612 IF(J.EQ.JJIGO TO 615 
GO Td 611 
613 BM(N,K)=BNG(I, J) 
GO TO 616 
614 AA=(DABS{PHIT(IN,K)))-CURVII yJ-1) 
BB=BMG(1,J)-BMG(I,J-1) 
CC=CURV(I,J)-CURVII,J-1) 
BM(N,yK)=BMG (I, J—-1)+AA*BB/CC 
GO TO 616 
611 CONTINUE 
Kyl 
ese ORME ie lax, | CURT ATURE AT PCSITION NO.*,I3," AND LOAD NO.',13,° 
11S BEYOND THE GIVEN CURVE NO.',13) 


GO TO 604 
616 BM(N,K)=(0.85*FCL (M) *W*(T(N)**2)*BM(N,K))/1000.0 


IF{PHIT(NsK)-LT~0-0)GO TO 630 


GC TO 604 
625 BM(N,K)=0.9 
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GD) T2 694 
BM{(N,K)=(-BA(N,K) ) 
CONTINUE 


2 CONTINUE 
5 CONTINUE 


WRITE(6,707) 
FORMAT(CLHO,50X,"*MOMENTS FCR KNOWN P IN KIP INCH*//) 
WRITE(6 ,708)(KyK=2yNL) 


3 FORMAT(2X,*POSITION NO./LCAD NO’ *y3X_13,9T10/( 26X,13,9110) ) 


DO 617 N=1,NS 
WRITE(6,709 Ny ( BMIN GK) »K=2,NL) 
FORMAT(9X,_,13,10X,10F10.2/(22X,10F10.2)) 
CONTINUE 

RETURN 

END 
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